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Erganzung zu dem Aufsatz von N. Scholz, Beitrage zur Theorie 
der tragenden Flache. 


Von P. F. Byrd. 


N. Scholz* hat gezeigt, daB die Streifenmethode in Verbindung mit dem 31/4-Theorem fiir eine 
ebene tragende Flache bei einer Aufteilung der Flache bis zu vier Streifen die exakten Werte des 
Gesamtauftriebs und des Gesamtmoments gema der Birnbaumschen Theorie der tragenden 
Wirbelflache liefert. Bei einem Kreisprofil dagegen erhalt man nach der Streifenmethode zwar 
den Gesamtauftrieb exakt, das Gesamtmoment aber nur noch in einer (allerdings ausreichenden) 
Naherung. Im folgenden werden diese Ergebnisse auf den allgemeinen Fall ausgedehnt, daB die 
tragende Flache in beliebig viele (m) Streifen aufgeteilt wird. Damit werden zugleich explizite 
Formeln [Gleichungen (16) und (18)] fiir die Teilzirkulationen gewonnen, die bei numerischen 
Rechnungen mit Vorteil verwendet werden. 

a) Die ebene Platte. Teilt man die Tragflache in n Streifen gleicher Tiefe auf, so lautet 
die Strémungsbedingung? im 31/4-Punkt des j-ten Streifens 


n c ms l ; 
pad ere mer G=1,2,...,n) (1) 
mit 

i SS ithe i; =2j-1. (2) 
Dabei ist G; die dimensionslos gemachte Teilzirkulation des i-ten Streifens. Sollen die durch (1) 


bestimmten Zirkulationen G; die exakten Werte des Gesamtauftriebs und des Gesamtmoments 
liefern, so miissen sie zusadtzlich noch den Bedingungen ® 


Ca Str =a ae es 
is ee :) 
und. 
Cm Str 1 = 
ie an ee Ce (4) 
genugen. 


Im folgenden wird bewiesen, daf die Lésungen G; von (1) diese Kigenschaft besitzen. Hierzu 
lésen wir zuerst die Gleichungen (1) nach den G; auf und verifizieren dann die Gleichungen (3) 
und (4). 

Zur Vorbereitung betrachten wir die Polynome 


f (@) =@—n) (2m) -+-(%—m), | 


| (5) 
& (2) =(% —t) (2 —ty) .. . (@ —th) 
und bilden damit die rationalen Funktionen 

Px (2) ey Ae fir ie O 1 2 (6) 


f (2) 


mit den sofort anzugebenden Partialbruchzerlegungen 


n : 1 
nait See 0 
t=1 


1 N. Scholz, Ing.-Arch. 18 (1950), S.84, insbesondere 5. 88f. 
2 Dies ist die Verallgemeinerung der Gleichungen (4.4) von N. Scholz. ; 
3 Gleichungen (4.5) und (4.6) bei N. Scholz. Die Gleichung (4.6) enthalt einen Druckfehler: die rechte 
Seite lautet n statt 1. 2 
99 
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<7 (7;) < r8(7) 1 ; 8 
n= —a+ Sth 4 Dae (8) 


MAY =FAQ=P 1% 2 pat 2. Pa) ) ae Fe) An 


| : &(vi) | S ; 8 (Ti) “gis (9) 


die sich noch vereinfachen lassen. Bezeichnen wir nimlich das Residuum! einer Funktion y (2) 
an einer Stelle a kurz mit Res(y; a), so ist nach bekannten Satzen der Funktionentheorie und 


wegen (2) 


788) ae eae : 
a Ga) 2 Res (% 3 r;) = — Res (3 ©) “, (10) 
und ebenso 
ay = — Res (9; 0) = — = (n +3), | 
a | (11) 
saves ; = ye el l 
Fn) = — Res(9y; 0) = — 5 (n*+3n-+4), 
wie man unschwer feststellt. Wegen ¢;(t;) =0 folgen damit aus (7), (8) und (9) die Beziehungen 
tte) iE “ 
cont f’ (Ti) Ciel ? 
TiS (0; eee ee 3 
relia) CE eae (13) 
r} 8 (Ti) pitohye fone (n+ 3) Gane eo ie ON (14) 
al if Gy a J Lowe oe, 


Nach diesen Vorbereitungen bemerken wir, das die Determinante von (1) 


A ={a;;\ ni, po = d 


a tr 


nicht verschwindet. Multipliziert man namlich die Zeilen von A nacheinander mit f(t,), f(t), .-., 
f(in), so erhalt man eine Determinante 


B= A f(t) f(t). - +f (tn), 


n 


deren Elemente 


b;; — IT (t; — Ty) 
eA 


homogene Formen (n— 1)-ten Grades in den GréBen 7,,..., t,t, +--+» ty sind; (die r;, t; werden 
also vorlaufig als Veranderliche angesehen). Die Determinante B ist dann in diesen GréBen eine 
Form vom Grade n(n— 1), die iibrigens verschwindet, wenn ein t,t, oder ein r,=T, gesetzt wird, 
weil dann in der Determinante zwei Zeilen bzw. zwei Spalten iibereinstimmen. Mit Hilfe eines 
bekannten Satzes der Algebra iiber die Darstellbarkeit eines Polynoms als Produkt seiner Wurzel- 
faktoren folgt vollends leicht, da8 

n(n—1) 


(-1) * 2(@—t)/) I @—r)) 
a tj i<j 
F(t) Ff (te) « - -f (tn) 


ist. Wegen (2) ist also sicher A von Null verschieden. 
Man kénnte nun ebenso die Unterdeterminanten A;; berechnen und die G; mit Hilfe der 
Cramerschen Regel aus (1) bestimmen. Da es aber wegen 4+0 genau eine Lésung von (1) 


ise (15) 


* Der Beweis 14Bt sich auch rein algebraisch unter Beniitzung von gewissen Reihen aus Binomial- 
koeffizienten fiihren. Den oben ausgefiihrten, wesentlich kiirzeren Gedankengang verdanke ich einem Hin- 
weis von Herrn Dr. K. Zoller, Stuttgart. 
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gibt, liefert der Vergleich von (1) mit (12) sofort 


uD 1 g(r) . 
G; => nn f(r) (i= the ae elieliels n) . (16) 


Nach (10) und (11) wird nun 


und damit 


Lex 
oe (Qres8\G.e= 1, 
Die G; (16) geniigen also den Bedingungen (3) und (4). 
b) Das Kreisbogenprofil. Hier sind die Teilzirkulationen Liésungen des Gleichungs- 
systems 


> Sere (pets) ee (17) 


t;—1; n2 
P= J i 


dessen Determinante wieder A (15) ist. Durch geeignete Zusammenfassung von (12) und (13) 
ergibt sich jetzt 


eat 2 SICA «Oyen 2r;—3 g(7;) — 
ee Gee Gee gw ee oS 8 
Die G; (18) lefern den exakten Wert des Gesamtauftriebs, wenn sie wieder der Gleichung (3) 
geniigen. Die Bedingungsgleichung fiir das Gesamtmoment lautet dagegen nunmehr1! 


eee 2 Y On — 3G =1. (4’) 


Cm ex 2n = 


— 


Aus (10) und (11) folgt mit (18) 


wahrend 


wird. Die Bedingung (3) ist also wieder erfiillt, die Bedingung (4’) aber streng nur in der Grenze 
fiir n—> oo, praktisch in ausreichender Naherung von n= 3 ab’. 


(Eingegangen am 2.Marz 1951.) 
Anschrift des Verfassers: Paul F. Byrd, Palo Alto Calif. (USA) 2655 Second Street. 


1 Man fiihre in Gleichung (3.2) bei N. Scholz an Stelle der I; die G; ein und dividiere durch den exakten 


Wert des Momentenbeiwerts Cm. 
2 Vel. Zahlentafel 1 bei N. Scholz. 
22% 
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Behandlung ebener Elastizitatsprobleme mit Hilfe 
hyperkomplexer Singularitaten. 


Von Kurt Schmidt. 


1, Einleitung. Die ebenen Probleme der Elastizitatstheorie werden meist* mit Hilfe der 
Airyschen Spannungsfunktion, des Inversionsverfahrens oder durch komplexe In- 
tegration gelést. Weniger bekannt diirfte sein, daB L.Sobrero® einen anderen interessanten 
Weg gefunden hat, die gesuchte Spannungsverteilung aus einer hyperkomplexen Funktion ab- 
zuleiten. Diese Methode steht zu der Fragestellung der ebenen Elastizitatslehre in einem ahn- 
lichen Zusammenhang wie die Algebra der komplexen Zahlen zur Potentialtheovie. Nach einer 
Darstellung dieses Verfahrens, das in der vorliegenden Arbeit® auf den ebenen Spannungs- 
zustand angewendet werden soll, wird die hyperkomplexe Spannungsfunktion fiir nachstehende 
Belastungsfalle entwickelt : 


a) Punktlast in der unendlich ausgedehnten Scheibe, 

b) Punktlast im Innern einer Halbscheibe, 

c) Punktlast in einem Parallelstreifen, 

d) zwei diametral entgegengesetzt wirkende gleich grofe Krafte im Innern einer Kreisscheibe. 
Soweit fiir diese Falle Liésungen bereits bekannt sind, werden sie mit den hier gewonnenen Ergeb- 
nissen verglichen. 


2. Die Gleichungen des ebenen Spannungszustandes, Wird eine Scheibe in der (x, y)-Ebene 
eines rechtwinkligen kartesischen Koordinatensystems nur durch Krafte beansprucht, die in der 
Ebene der Scheibe wirken, so gelten, sofern vom Einflu8 des Eigengewichts und damit der Massen- 
krafte abgesehen wird, bei Verwendung der iiblichen Bezeichnung der Spannungskomponenten 
mit 6,, oy und t die Gleichgewichtsbedingungen 


007%. Or Oo Crs 
a Oy =e an =0. (1) 


Mit den Verschiebungen &, 7 in der x-, y-Richtung lauten die Gleichungen fiir die Dehnungen 
Ex, €y und die Gleitung y 


at 1 
eae (ee Oy)’, | 
Oy eel 
oer (0, — VG), (2) 
on ide 1 
Ox ' Oy G 


Hierin sind FE der Elastizitatsmodul, G der Schubmodul und » die Querzahl, die durch die Glei- 
chung E—2(1-+-») G miteinander verkniipft sind. Die Elimination der Verschiebungen aus (2) 
liefert die fiir die Verzerrungskomponenten geltende Vertraglichkeitsbedingung 


Ot, | Oey O'y 
dy? | Ox? dAxdy~ 


Werden in dieser Gleichung die Verzerrungskomponenten durch die Spannungen nach den Glei- 


chungen (2) und die Schubspannung mit Hilfe (1) durch die Normalspannungen ersetzt, so folgt 
die Gleichung 


OG, 0, eo ao 
dat 1 Oe 1 Ge Taye 40. 1 9,) =0, 


* Vel. etwa L. Féppl, Drang und Zwang, Bd.3. Miinchen 1947. 

2 ZL. Sobrero, Memorie della Reale Accademia d’Italia 6 (1934/135), S.1; oder auch L. Sobrero, Theorie der 
ebenen Elastizitat. Hamburger Math. Einzelschriften, Heft 17. Leipzig 1934. 

® Gekiirzte Fassung einer von der Technischen Hochschule Hannover genehmigten Habilitationsschrift 
(D 89); Referent: Prof. Dr.-Ing. 0. Flachsbart; Korreferent: Prof. Dr.-Ing. habil. A. Pfliiger. 
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die die Vertraglichkeitsbedingung durch die Spannungen ausdriickt. Die Elimination von 7 aus 
(1) liefert 
Go, doy 


ot ye 
womit aus der vorhergehenden Gleichung fiir o, schlieBlich folgt 
of Ox o4 Ox of Ox 
Oxt | Ox? ay? | dyt a A Ao, = 0 a 


Mit den in gleicher Weise abgeleiteten Gleichungen fiir oy und t ergeben sich fiir die Spannungen 
die drei Bipotentialgleichungen 
AAG =O, AAgy— 0, AAr=0. (3) 
3. Das Sobrerosche Verfahren mit hyperkomplexen Funktionen, a) Allgemeine theo - 
retische Grundlagen. Die GréBe 
foatjotjret jd 
soll als hyperkomplexe Zahl definiert werden, deren Elemente a, b, c, d beliebige Funktionen der 
vier reellen Veranderlichen x, y, u, v sind. Die Konstante j und ihre Potenzen sollen die Ordnung 


der reellen Funktionen innerhalb der hyperkomplexen GréBe kennzeichnen; dann ist f eine Funk- 
tion der hyperkomplexen Verdnderlichen 


s=a+jytPutfo, 


d.h. es ist f=f(z). Fir den hier interessierenden ebenen Fall verschwinden u und v. Uber die 
Konstante j wird derart verfiigt, daB die Differentialgleichung 


OLN ity ero tnee Cw nee, 
axt | “ax dy? | Ay4 
mit z=x-+/y erfillt ist. Dies fiihrt zu dem Multiplikationsgesetz 
I 25 ae J 0) (4) 
das die Potenzen j*, j°,... und héhere durch die ersten vier Potenzen j°= 1, j, j? und j® auszu- 


driicken erlaubt. Von Interesse ist weiter die leicht nachzurechnende Relation, die die imaginare 
Einheit i mit den hyperkomplexen Termen verkniipft, naimlich 


at a ey). 


Beziiglich des Rechnens mit hyperkomplexen Zahlen und Funktionen sei auf die erlauternden 
Ausfithrungen des Anhanges unter Ziffer 6 verwiesen. Die Voraussetzung der Differentiierbar- 
keit der Funktion f fihrt zu 


af dfdz df af df az df. 


Qn dzdx dz? dy dz dy dzl° 
Die Unabhangigkeit des Differentialquotienten df/dz vom Wege ergibt 


oder 


A+ PtP +P a= (bib tPothay), 


wobei die Indizes die partiellen Ableitungen kennzeichnen. Diese hyperkomplexe Gleichung liefert 
vier reelle Gleichungen zwischen den ersten partiellen Ableitungen der Funktionen a, b, c, d. 
Werden die Elemente gleicher Ordnung mit Hilfe von (4) zusammengefaBt, so folgen die Cauchy- 
schen Bedingungen 

—ce,+d,=0, a,+d,=0, a,—b,=0,  b,—¢,—2¢d,—0. (5) 
Werden in diesem Gleichungssystem je drei der Funkionen a, b, c, d eliminiert, so ergeben sich 
nacheinander aus der jeweiligen vierten iibrigbleibenden Gleichung die Bipotentialgleichungen 


Ada=0, AAb=0, Ade=0, AAd=0. (6) 


Der Vergleich zwischen den beiden Gleichungen (1) und den ersten beiden Gleichungen (5) zeigt 
die Aquivalenz zwischen den Beanspruchungskomponenten des Spannungszustandes und dem 


326 Schmidt: Behandlung ebener Elastizitatsprobleme usw. Ingenieur-Archiv 


ersten, dritten und vierten Element einer Funktion der hyperkomplexen Veranderlichen z= x+-Jy, 
und zwar gilt 
O,=—C, Cy ==, aC 


Es darf somit gesetzt werden 
f(2)=oyF jO—-P ort. (7) 
Gleichzeitig entsprechen auch die Differentialgleichungen (3) und (6) einander. Dieser Zusammen- 
hang zwischen den Elementen einer hyperkomplexen Funktion und den Spannungskomponenten 
fiihrt die Lisung spezieller Probleme der ebenen Elastizitatstheorie auf die analytische Aufgabe 
guriick, eine Funktion f(z) mit z—=-x-+ ] y eines ebenen Systems zu 


~ Oar ermitteln, fiir die bestimmte Singularitaten und Randbedingungen 
; vorgeschrieben sind. 
= ae b) Speziellere Ausfiihrungen. Sind innerhalb einer Scheibe 
ij he fir ein Element dz—dx-+-jdy einer Kurve die den elastischen 
Ordy Spannungszustand entsprechenden Komponenten einer Kraft dP 
a £ und dQ parallel zur x- bzw. y-Achse, so verlangt das Gleichgewicht 


Abb. 1. Zur Gleichgewichtsbetrach-  (vg]. Abb. 1) 

tung innerhalb einer Scheibe. pains of dy, dQ Les oy ages dy ; 
Dabei sind die beiden Komponenten auf die Dicke der Scheibe bezogen, die gleich der Langen- 
einheit gesetzt wird. Andererseits stellt die hyperkomplexe Funktion f(z) die elastische Bean- 
spruchung des betreffenden Punktes dar; dies liefert mit Riicksicht auf (7) und Beachtung von (4) 


f (2) dz= (a, + jb —j°o, + Jt) (dx +-j dy) 
= (o,dx —rdy) +-j(0,dy + bdx) + j?(bdy —2tdy —o,dx) -j*(tdx —o,dy). (8) 


Das erste und das letzte Element dieser Gleichung entsprechen den Kraftkomponenten dQ und 
dP. Das Integral [f@ dz langs eines endlichen Bogens der Kurve entspricht dann den Kompo- 


nenten Q und P, die infolge des elastischen Spannungszustandes von der einen Kurvenseite auf | 


die andere iibertragen werden. Ist der Integrationsweg geschlossen und wirkt innerhalb des ein- 
geschlossenen Bereiches keine auBere Kraft, so hat das Integral den Wert Null. Wird dagegen in 
das umschlossene Gebiet eine aufere Kraft eingeleitet, die fiir die Funktion f(z) eine Singularitat 
bedeutet, so verschwindet das Integral ff (2) dz nicht mehr. Das erste und das letzte Element 
entsprechen dann den Komponenten Q und P der angreifenden Kraft. Der Cauchysche Satz und 
der Residuumbegriff behalten damit fiir die hyperkomplexe Funktion f(z) ihre Giltigkeit. Greift 
nun im Punkte z eine auBere Kraft mit dem Komponenten P und Q parallel zur x- bzw. y-Achse 
an, so 1aBt sich der durch diese Kraft hervorgerufene elastische Spannungszustand durch die 
Funktion 


darstellen, wobei k eine willkiirliche hyperkomplexe Konstante 


k= a+jp+Py +70 
ist; a, f, y, 6 sind dabei reelle GréBen. Der Wert I des Integrals langs eines geschlossenen In- 
tegrationsweges fiihrt fiir die gewahlte Funktion zu 


I= Of(z)dz= 2nik. 
Einzelheiten iiber die Auswertung dieses Integrals sind in Gleichung (23) des Anhangs unter Ziffer 6 


wiedergegeben. Wird die imaginare Einheit i= S (37 +-j*) gesetzt, so folgt nach einfacher Um- 
formung 


T= $ f(s) dz = x[— (B+ 6) +j (8a — y) +7? (6 — 36) +7 («+ y)). (9) 


Dieses Ergebnis liefert mit der integrierten Gleichung (8) fiir die Komponenten der Kraft 


P=nx(aty), Q=—x(6+8). (10) 
Die Funktion f(z) entsprach nach Gleichung (7) in ihrem ersten, dritten und vierten Element den 
Komponenten des elastischen Spannungszustandes. Diese Elemente sind daher physikalisch un- 
mittelbar anschaulich. Dies gilt jedoch nicht fiir das zweite Element, das durch die Funktion b 
ausgedriickt wurde; der Zusammenhang zwischen diesem Element und den Spannungskompo- 
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nenten geht aus Gleichung (5) hervor. Danach ist 


Clom E are 92 
b,, = 2 ee ee RS as sll L 
pasate dy Be 1— (sey +758) aes 
Be Oty ae eer O2E 0”) 
Cate es ral? oe Y Be iy) ; 


wobei die Spannungen durch die Verschiebungen ausgedriickt wurden. Aus den Ausdriicken fiir 
b,, und 6, 1a8t sich nun die Funktion b ermitteln. Aus (1) folgen zunachst unter Beachtung von (2) 
fiir die Verschiebungen die Differentialgleichungen 


Cn 1 are 1 ; On 
Ox? oy ) Oy? 2 ai) Ox dy’ 
OP apes 1 pire ozé 1 ony. 
oy? e) Us a 5a5, D (Cu Ox? ’ 
damit wird 
dg oe an 
— = | | 
= 21+ ae ) Ox oy 53| Pad,, 
es dB: aé an 
Ly 2) ee » ay? oF an ay 
Integriert folgt 
= PE reoti en. os an| | ” 1__E 0& Lay On 
en aie a, a ao Seer (on pq) 


Dabei sind die auftretenden willkiirlichen Funktionen, die hier unwesentlich sind, fortgelassen. 
Dieses Ergebnis liefert fiir eine Integration nach (8), wenn z,= §&,-+-j 7, ein Punkt mit verschwin- 
denden Verschiebungen sein soll und z= €-+-j7 ein beliebiger Punkt der Ebene ist, 


Z 
did 13a | eal 0& | OM t jeé —0n 
f{e\d= ee Fr | real Ox | 5y )a oy (ay | a 


zy zy 
E 1 0& 1 on e 0& 2+ v On ; 
J 1+ Ai Fear Da * ay) i) 2 dy 2 524 
a (11) 
pan dD; yn 2+ 48 1 0g on 
2 ¥ = See ie | 
J rept |(s Ox 2 33) oy el Ox | 5) a) 


wn) Wo; ; If GE, Of) i y/ees on 
3 wei f 
a ee ale reall: Hoy) a 
zy 

Aus dieser Gleichung lassen sich nun die Verschiebungen (€— &,) und (7—7,) ableiten. Da fiir die 
Verschiebung in x-Richtung nur das erste und dritte Element der rechten Seite von (11) Beitrage 
liefern, mu ein hyperkomplexer Multiplikator A+-j?B, wobei A und B reelle Zahlen sind, fiir 
die Gleichung (11) so bestimmt werden, daf fiir das erste Element das totale Differential 


S205 ag 
fe pp att 
entsteht. Diese Bedingung fiihrt zu der Gleichung 


icy i ea oy ae ee 
= (7544 B) 5a (34 B) dy, 


1—y» 2 Oy 


1 : 
die Ee Sel und —————— liefert. 
l—y»v 9 


Die Verschiebung in x-Richtung ergibt nunmehr das erste Element der Gleichung 
zz 
I , 
SS ae (12a) 
oat 


Wird in entsprechender Weise Gleichung (11) mit C-|-j? D multipliziert, um das totale Differential 
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jdy, fiir das jetzt nur das zweite und vierte Element von (11) Beitrage liefern, zu erhalten, so 
nehmen die reellen Konstanten C und D die Werte 

vot ee 

or 1+’ 2 I+» 


an. Die Verschiebung in y-Richtung folgt dann aus dem zweiten Element: 


SFG) ders tad emer iced (12b) 


Mit den beiden Gleichungen (12a) und (12b) lassen sich die Verschiebungskomponenten errech- 
nen. Weiter folgen damit fiir die reellen GrofSen a, B, y, 6 der hyperkomplexen Konstanten k die 
Beziehungen, die zwischen diesen Gréen untereinander bestehen. Werden namlich in (12a) 
und (12b) die Integrationen itber eine geschlossene Kurve erstreckt, so miissen die Verschiebungen 
£—& und 7—y7, verschwinden. Dies fiihrt mit Hilfe von (9) zu den Gleichungen 


Bl—»r)—é6(3+7)=0, a(3+7)—y(l1—»v)=0. 


Mit Gleichung (10) ergibt sich schlieBlich die Abhangigkeit der GréBen a, B, y, 6 von den Kom- 
ponenten P und Q einer auBeren Kraft aus 


a= -—— P, 


Damit ist die hyperkomplexe Funktion f(z)= k/(z— 2), die den elastischen Spannungszustand 
fiir eine im Punkt z angreifende Kraft mit den Komponenten P und Qin x- und y-Richtung dar- 
stellt, vollstaindig bestimmt. Mit den Ergebnissen fiir a, 6, y, 6 nimmt die hyperkomplexe Kon- 
stante k den Wert 


k= [1 — 9) + 723 L »)] ERIE HY) ard (13) 


an, wobei besonders hervorzuheben ist, daB die zugrunde gelegte Singularitat einen Pol ersten 
Grades besitzt. 

c) Das Prinzip der analytischen Spiegelung. Fir die Anwendung auf spezielle Pro- 
bleme ist das Verhalten der Funktion f(z) am Rande eines ebenen Systems von besonderem In- 
teresse. Wird z.B. ein freier Randteil betrachtet, der zur x-Achse parallel verlauft und an dem 
keine auBeren Krafte angreifen, so verschwinden in den Punkten dieses Randes die Spannungen 
oy und t, d.h. die Funktion f(z) darf hier nur aus dem zweiten und dritten Element von (7) be- 
stehen. Im Hinblick auf zu lésende Randwertaufgaben wird zum Aufbau einer Funktion f(z), 
die in der oberen Halbebene des (x, y)-Systems im Punkte z=] a, a sei eine gegebene reelle Zahl, 
einen Pol besitzen mége, die Singularitat vom Typ 1/(z— j a) vorgegeben; dabei ist zunichst der 
Zahler, der im allgemeinen eine hyperkomplexe Konstante sein wird, gleich der reellen Einheit | 
gesetzt. Um die Funktion f(z) zu bestimmen, liegt es nahe, sich der Methoden der Potential- 
theorie in analoger Weise zu bedienen. Befindet sich namlich eine Einheitsmasse im Punkte 
%9—=jades Feldes und wird weiter lings der reellen Achse eine Massenbelegung, die durch eine un- 
bekannte hyperkomplexe Funktion w(x) darstellbar sei, angenommen, so ist das ,,Potential* 
dieser Massen 


+ 00 


hen ng Ey whee E) 
fol) = 35 J | Ee dx. 
Fitr einen Punkt ¢ der reellen Achse nimmt mit z—> € die Funktion f,(z) den Wert 
ty fe [as 
fol pag | See ete 


an. Da die Funktion (x) auf der x-Achse nur zweite und dritte Elemente enthalten soll, wird 
auch der mittlere Term dieser Gleichung nur aus solchen bestehen. Fiir einen Ansatz 

[u(x) = jA(x) + J? B(x) 
geniigt es dann, die reellen Funktionen A(x) und B(x) derart zu bestimmen, daB das erste und 
letzte Glied der Gleichung die vorgeschriebene Randbedingung, keine ersten und vierten Ele- 
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mente zu enthalten, erfiillen. Aus 


_ 


a ieweuee, eae oe ia, 


x%—>— 0 


folgen damit fiir den Ansatz (¢)—=jA(¢)-+-j?B(¢) und der in seine Elemente zerlegten Funktion 
1/(¢—ja) unter Beachtung von 


xz (35+ 7) u(t) =% (—A—jB+ fat pB) 


die Gleichungen 
Ci eee | 2 a 
f2+ g2 | (24 a®)2 “9 A=0, (C2 + a®)? begs: 
Nach Eliminierung von A und B ergibt sich 


u()=j4@)+ PB) =—izlatat rep 


und weiter 


sss: en Newer yreesetee eer ee 


 2—ja 


Die Integration kann mittels Partialbruchzerlegung oder bequemer mit Hilfe des Residuensatzes 
ausgefiihrt werden. 

Wird namlich statt x voriibergehend die im allgemeinen aus vier Elementen bestehende hyper- 
komplexe Veranderliche w eingefiihrt und der Integrationsweg auf der reellen Achse durch einen 
Halbkreis derart geschlossen, daB die reelle Achse und ein vorgegebener singularer Punkt w= wy, 
eine im vierdimensionalen Raum liegende Ebene festlegen, so gibt die Integration auf dem in dieser 
Ebene liegenden Halbkreis keinen Beitrag, da fiir beide Terme des Integranden der Nenner von 
héherem Grad ist als der Zahler. Unter Beachtung der beiden singulaéren Punkte w)=z und 


o. = (37 + j°)a = ia liefert der Residuensatz fiir den ersten Term 


z La 1 z+ita 
—» = d tip ieee == = 

Res fr (W).w,—« 2+ a2 Vee Res fr (W)w,—i0 2 ia(z—ia) 2 22+ a? 

Der zweite Term des Integranden wird an denselben Punkten singular. Das Residuum bei w= =z 

nimmt den Wert an 


ja—z 


Res fr (W).9=s =a (2 +02)? ° 


Zur Berechnung des Residuums an der Stelle wy—=ia wird eine Potenzreihenentwicklung fiir die 
ein wenig umgeformte Funktion 
a® (j — i)— a? (w— ia) 
(w— ia)? [(w +i a)? (w— 2)] 
durchgefiihrt. Die Entwicklung der eckigen Klammer im Nenner fiihrt mit 
f(ia)=40?(z—ia) und  f’(ia)= 2(2 ia) (ta — 2) + (214)? 


zu der Darstellung 


(0 fia) (0 —2) lapse = fa) + EFM (ote) = 4a%(2—ia) 1 — FEAT (wo —ia)| 
und damit zu dem Residuum 
Res fir Wain = 5 eG sin Pa a (ores) aie 
WeGa Gees 9 ot SF 
~~ A(z— ia)? 4(z—ia) ° 


Zusammengefaft ergibt sich fiir das Integral 
mw a+ (2t-+ j)2 ai 


oy ot a? 


Zits Kes=—= 
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Somit wird 
1 A Ps 22 I 3a2z 8 


4 a 3a L 2 ee. a | 
fo(2) = Pra (2+?) ed 2a TF + a?)?} I leg | (Geaiman)e J Cems || 


ee 


Da die Funktion f(z) in der oberen Halbebene nur im Punkte z ja einen Pol haben kann, liegt 
die Vermutung nahe, daf sich das Ergebnis mit Hilfe von Funktionen, die in der unteren Halb- | 
ebene Pole besitzen, umformen aft. Ein Ansatz mit hyperkomplexen Konstanten fiihrt durch — 
Koeffizientenvergleich zu 


et 277—1 iE Be ee ack 
Jol eoaeier faa + 2 Cay 2a Gagan 


Dieses Resultat findet in der Potentialtheorie eine gewisse Analogie in dem ,,Schwarzschen 
Spiegelungsprinzip“, mit dessen Hilfe sich bekanntlich durch Hinzunahme von zum Rande 
symmetrisch aufgebauten Funktionen vorgegebene Randbedingungen erfiillen lassen. Aus der 
letzten Gleichung ist ersichtlich, da® dem Pol ersten Grades der Funktion 1/(z—j a) im sym- 
metrischen Punkte ein Pol vom dritten Grade zugeordnet werden mu, um die Randbedingungen 
auf der x-Achse zu befriedigen. Das in der komplexen Funktionentheorie benutzte Spiegelungs- 
prinzip erhalt somit hier eine erweiterte Fassung. Fiir die Anwendungen der Theorie sind weiter 
die Funktionen von Bedeutung, die im Punkte z ja einen Pol ersten Grades besitzen, wenn die 
Singularitaten innerhalb der Scheibe durch die Funktionen 

ea Hass te ee 

ofa * 3—— 4) (5 (Si) (0) 
vorgegeben sind und weiter verlangt wird, daB auf der reellen Achse das erste und das letzte Ele- 
ment verschwinden sollen. Die Durchfiihrung der Rechnung erfolgt in gleicher Weise wie die bei der 
Singularitat 1/(z—ja); es soll daher hier nur das Ergebnis der Rechnung zusammengestellt werden: 


xi YH 2j2—1 Tie es aid 
fo(2) ee 22 an eee (z+ jaye” 

waacher | ‘vl igen pales 
fi (2) ery z+ja 1 Oe gaan (2-=- ja)?” 14 
fi =— F420 2 pete a 
aja "a+ ja (@+ja) Sa 

ei 1 2 ITP 

aed Ii 2 le 0 

Js(2) ecanee aaa rin ce air sas 


Mit Hilfe dieser Gleichungen gelingt es, lineare Kombinationen mit reellen Koeffizienten dieser 
Funktionen anzugeben, die im Punkte x0, y—a der (x, y)-Ebene eine Singularitat ersten Gra- 
des besitzen und auf der x-Achse nur zweite und dritte Elemente enthalten. 


4, Lésung spezieller Spannungsprobleme, a) Punktlast in der unendlich ausgedehn- 
ten Scheibe. In der unendlich ausgedehnten Ebene einer Scheibe greife im Punkte x=0, y=0 
des Koordinatensystems eine Kraft P an, deren Wirkungslinie mit 
der negativen Halbachse von x zusammenfallt (vgl. Abb.2). Die 
Kraft P ist wieder auf die Dicke der Scheibe bezogen, die gleich 
der Langeneinheit gesetzt sei. Die dieser Kraft das Gleichgewicht 
haltenden auBeren Krafte sollen im Unendlichen liegen und sind 
daher fiir den Spannungszustand im Endlichen bedeutungslos. Nach 

AVN oO actin dee -Vollschaibe. den Ausfiihrungen des vorhergehenden Abschnittes kann der durch 

P hervorgerufene Spannungszustand durch die Funktion f(z) = k/z 

mit der hyperkomplexen Veranderlichen z—.x-|j y beschrieben werden. Nach Gleichung (13) 
ist die hyperkomplexe Konstante fiir die hier vorliegende Belastung 


k= —F (a—» + PO +). 

Damit wird nach Gleichung(7) die Spannungsfunktion 

Pid= vy) 77 (359) 
A mt Cela: 
deren Aufspaltung mit Hilfe (22) die bekannten Werte der Spannungen!? liefert. 


f(@) = oy + jo—fPo, + Ppr= 


2) 


* Vel. etwa K. Girkmann, Ing.-Arch, 11 (1940), S.418 oder auch K. Girkmann, Flachentragwerke, S.108, 
2. Aufl. Wien 1948. 
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Die Komponenten & und 7 der Verschiebung eines Punktes der (x, y)-Ebene berechnen sich 
nach den Gleichungen (12a) und (12b). Danach wird 


z 


EE ps4 ( 


eer 1 7 
VB) f $G)de a (— 9+ Pyke 
Nach Kinsetzung des oben angegebenen Wertes fiir k und der Zerlegung des Logarithmus in seine 
Elemente nach Gleichung (23) ditrfen nur die reellen Terme herangezogen werden. Die Zu- 
sammenfassung liefert 


e Tee lty [ y? an G3 
5 i= ie c | 0) Fale | 5 In(x?+- y?)|. 


In entsprechender Weise ergibt sich fiir die Komponente der Verschiebung in y-Richtung 


“ PUES?) a ae Soe 
ie Ween ar was 5 xy’ 


dabei sind &, und 7, Integrationskonstanten, die eine Verschiebung der Scheibe als Ganzes be- 
riicksichtigen. Die beiden erhaltenen Ausdriicke stimmen mit den Angaben von A. und L. Féppl! 
iiberein. 

b) Punktlast im Innern einer Halbscheibe. Der Rand der Scheibe falle mit der 
x-Achse zusammen, so daf die obere Halbebene die Halbscheibe darstellt. Im Abstand a vom Ko- 


ordinatenursprung (vel. Abb.3) liege auf der y-Achse der Angriffs- ; 
prung (vg ) lieg y-Ac 8 IN NN 
CS \ ~ 


punkt der auBeren Last, die in zwei Komponenten parallel zu den \ 
SEQ 


Koordinatenachsen zerlegt werde. Die Komponente P sei parallel 
zum Rand gerichtet und weise wieder in die Richtung der nega- 
tiven x-Achse; die zum Rand senkrecht verlaufende Komponente 
Q falle in die Richtung der negativen y-Achse. Mit Hilfe der 
Ergebnisse unter Ziffer 3 14Bt sich die hyperkomplexe Funktion 
f (2), aus der die durch die Komponenten P und Qim Punkte z=ja 
hervorgerufene elastische Beanspruchung abgeleitet wird, in einfacher Weise angeben. Unter 
Beriicksichtigung der fiir die Komponenten gewahlten Richtungen wird mit (13) und der durch 
(14) eingefiihrten Bezeichnungen die Spannungsfunktion 


fQ@=—Z|d—)fAO4+64+40) 4+ 216+)£04+ 0-946). 


Abb. 3. Last in der Halbscheibe. 


Die Aufspaltung dieser Gleichung in ihre vier Elemente liefert die Werte der Spannungen, die im 
Einklang mit den Ergebnissen von E. Melan® und K. Girkmann? stehen. 

c) Punktlast in einem Parallelstreifen. In einem unendlich langen Parallelstreifen 
von der Breite a greife parallel zu den Randern eine Kraft Q an. Der Ursprung des Koordinaten- 
systems falle mit dem Angriffspunkt der Kraft zusammen; Q wirke in Richtung 
der negativen y-Achse (vgl. Abb.4). Die Abstande der y-Achse vom rechten und 
linken Rand seien durch die Langen a, und a, gekennzeichnet. Die Aufgabe be- 
steht nun darin, die Spannungsfunktion anhand der allgemeinen Ausfiihrungen 
unter Ziffer 3 aufzubauen. Danach sind im Koordinatenanfangspunkt Pole ersten 
Grades anzunehmen, die durch die Funktionen 1/z, j/z, j?/z und j3/z dargestellt 
werden. Fiir den rechten Rand xa, wird weiter verlangt, dais die Spannungen 
o, und t verschwinden. Zur Erfillung dieser Randbedingungen ist auf der 


Abb. 4, Last in 
x-Achse im Punkte x — 2a, ein Pol dritten Grades einzufiihren; die Koeffizienten _ einer streifen- 


: : 3 : : férmigen Scheibe. 
der dadurch bedingten Funktionen sind derart zu bestimmen, daf die dritten und Be ok 


vierten Elemente fiir den Rand x =a, Null werden. Das Ergebnis der Rechnung ist in dem 
Gleichungssystem (15) zusammengestellt. 

Wird nun auch fiir den linken Rand x = —a, verlangt, da die dritten und vierten Elemente 
verschwinden, so lassen sich aus (15) vier analoge Gleichungen entwickeln, indem fiir a, der Wert 
—a, eingefiigt wird. Ein Vergleich zwischen (15) und (14) zeigt den einfacheren Aufbau der Glei- 


chungen (15); die durch Spiegelung erhaltcnen Terme in (15) weisen namlich in den Zahlern nur 


1 4, und L. Féppl, Drang und Zwang, Bd.1, $.273. Miinchen und Berlin 1941. 
2 E. Melan, Z. angew. Math. Mech. 12 (1932), 5.345 und 20 (1940), S. 368. 
3 K. Girkmann, Ing.-Arch. 11 (1940), 5.424 und 5.422. 
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jeweils zwei verschiedene Elemente der’ Koeffizienten auf, wahrend in (14) drei verschiedene Ele- | 
mente auftreten. Fiir die weitere Darstellung bedeutet dies eine wesentliche Vereinfachung. 


Teed eres Mee ee nen 
Chae piceresrremn ryry a IN Oe Le 
j j ; bigidil: 
LO Sar psangaae ait (s—2a,)? i (15) 
ip Ue 2a ae 
Up ere aan ns 1(@—2a,)? ” 
P27 j ed 
f(@) = Mae (2—2a,)? ! i 


Zz %2—2a, 


(E20, )e © 


Im folgenden soll die Spannungsfunktion f(z), die die von der Kraft Q bewirkte Spannungs- 
verteilung angibt, aufgebaut werden. Nach Gleichung (13) ist die hyperkomplexe Konstante fiir 
diese Belastung 


k=p[iBty+F0—9)]. 
Im Koordinatenanfangspunkt ist somit von den Singularitaten j/z und j?/z auszugehen. Werden 
nun mit F(z) bzw. F’,(z) die hyperkomplexen Funktionen bezeichnet, die im Punkte x =0, y =0 
wie j/z bzw. j?/z einen Pol ersten Grades haben und deren dritte und vierte Elemente an den beiden 
Randern x —a,, x = — a, verschwinden, so gilt 


R@aite, 


wobei die Punkte die noch zu entwickelnden Funktionen zur Erfiillung der Randbedingungen an- 
deuten sollen; vorausgesetzt ist, daB die Art der Lastaufnahme im Unendlichen nach der schon 
von B. de Saint- Vénant erkannten Erfahrungstatsache die Spannungsverteilung in der Umgebung 
der Lastangriffsstelle nicht stért. Die durch die Kraft Qin dem Parallelstreifen auftretende Span- 
nungsverteilung la8t sich dann durch die hyperkomplexe Spannungsfunktion 


f= [6+)F (2) + (1— ») F(z) | (16) 


Ae ide 


wiedergegeben unter Kinfithrung der dimensionslosen Veranderlichen 2 = xj y=z/a. Die ein- 
gefiihrten Indizes kennzeichnen die Abhangigkeit von den Verhaltnissen A,= a,/a bzw. A,= a,/a. 

Zum Aufbau der Funktion F, (z) wird zunachst die Funktion f,(z) von (15) herangezogen. Die 
Spiegelung der Funktion j/z am rechten und linken Rande liefert fiir die erste Naherung der ge- 
suchten Funktion 

= J Joos ips ween ef 
F(z) == | - -} 2a: = 2hi = | 
z 


Vy MAE eoiy hh 2 (z+ 2A,)? a(z ode 
9 2 9 2 


In dieser Gleichung sind der Kiirze halber die neu aufgenommenen Funktionen zusammengefaBt, 
und zwar kennzeichnen die oberen (unteren) Vorzeichen und Indizes die durch Spiegelung am 
rechten (linken) Rande erhaltenen Funktionen. Die im Punkte x =2 4, bzw. —2 A, auf der x-Achse 
singular werdenden Funktionen rufen am linken bzw. rechten Rand eine Stérung der Rand- 
bedingungen hervor, die eine erneute Spiegelung erfordern. Diese fiihrt zur Bildung der zweiten 
Naherung fiir F(z), die durch Aufnahme von weiteren Funktionen, die auf der x-Achse in den 
beiden Punkten x = -+ 2(A,-+-A,) = -+ 2 singular werden, entsteht. Dabei sind die vier neu hinzu- 
gekommenen Funktionen der ersten Naherung einzeln zu spiegeln. So liefert beispielsweise der 
Term j/(?—24,) mit Hilfe der zweiten Gleichung von (15) unter Beachtung des Ortes des vor- 
gegebenen Poles den Ausdruck 

; a Snes 
sty t20.+ 4) ete 20 +4 2 
Sollen Funktionen héherer Ordnung in gleicher Weise gespiegelt werden, so sind — L. Sobrero} 
wies bereits darauf hin — die Gleichungen (15) zu differentiieren, woraus sich dann die gesuchten 
Beziehungen ergeben. Die beispielsweise fiir die vorgegebene Funktion J?/(2— 21,)2 anzuschreiben- 
den Terme folgen aus der dritten Gleichung von (15) mit 

= orcs 

aay 4(1 i) eras 6(1 A,)? oe 


* a,a,O. in der deutschen Ubersetzung S.28. 
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Sind alle Funktionen, die in den Punkten x —-+- 2A: auf der x-Achse singular werden, in dieser 
Weise gespiegelt, so ergibt dies zusammengefaBt 


ip 


yA b+242) +P + 643) __ 


+ sg F 2, 


ery 


(2>- 2)* 


6 j (2 Az +242) +38 (2 Az +2 23) 


(z+ 2)4 


Dieser Ausdruck zusammen mit der oben genannten ersten Naherung stellt fiir F(z) die zweite 


Naherung dar. 


Die neu hinzugekommenen Funktionen, die ihrerseits wieder die vorgegebenen 


Randbedingungen stéren, werden erneut gespiegelt; dadurch werden Polynome aufgebaut, die 
auf der x-Achse in den Punkten x= -| 2 (1 +4n) singular werden. Auf eine Angabe dieser Funktionen 


darf verzichtet werden; sie bilden mit de 


n bereits angeschriebenen die dritte Naherung fiir F, (2). 


Hervorgehoben sei nur noch, dafs die zweite Naherung einen Pol vierten Grades, die dritte Nihe- 
rung einen Pol fiinften Grades enthalt usf. Die Fortfiihrung des Verfahrens liefert fiir die Span- 
nungsfunktion F(z) eine Reihenentwicklung, deren Besonderheit darin besteht, daB jede weitere 


Naherung einen Term mehr besitzt als die 


vorhergehende. Ein Vergleich der einzelnen Naherungen 


untereinander zeigt einen gesetzmaBigen Zusammenhang zwischen den Naherungen ungerader 
Nummer und denen gerader Nummer, und zwar bilden die Zahler der Funktionen gleicher Po- 
tenz in den aufeinanderfolgenden Naherungen jeweils eine arithmetische Reihe. Diese Tatsache 


vereinfacht die Berechnung der Naherun 
dem die rant Darstellung d 


ps Ly = Wat [on F% | 2 


is Be 


gn—2 


gen hoherer Ordnung wesentlich; sie erméglicht aufer- 
er Funktion 

j 
zm 


— 


ea eae apa hae Steet) 
[4 m?—(n—2)2-4 7 (n—1)(2n—5) 
j(m—2)} — > (2-1) (m—1)(1—24,) 
hetero oD [12a 
aoe Z s : TF 
z+ [(m— Ay|\" (m—n +2 ie ote 2 | Sa! 
ee = rk we or + m?—(n—2)?4 3 (n—2) (2n—5) 
+ P(n— 1)| — (0 — 2) (m—1) (1-24) 
Epes ome ee 
ir (ol 2) ae 0) (2) 
j (n— 2)| —>(n—1)(m—1) (1—24.) 
1 ‘ 
ie + (nay I= (2498 
>: Ss gn—3 e 2 ): ; 4 | : 2) | - = 
eel re al iia eg aie) 
+ j®(n— 1)| — > (n— 2) (m—1) (1—2.4,) 
(m2) [1 — (245) 


(17) 
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dabei sind wieder obere bzw. untere Vorzeichen und Indizes einander zuzuordnen. Auf eine noch 
straffere Zusammenfassung des Ergebnisses wurde im Interesse der Deutlichkeit verzichtet. Der 
Gedanke liegt nahe, durch Umordnen der einzelnen Glieder der Naherungen untereinander zu 
einer geschlosseneren Darstellung oder gar zu Funktionen zu gelangen, die denen der komplexen | 
Theorie analog sind. In Anbetracht der technischen Fragestellung der vorliegenden Arbeit soll | 
hierauf jedoch nicht weiter eingegangen werden. Beziiglich der Konvergenzfrage der erhaltenen | 
Reihe geniigt der Hinweis, das die Reihe auf Grund des Prinzipes der analytischen Fortsetzung — 
konvergiert. Die Anzahl der erforderlichen Naherungen fiir eine vorgegebene Genauigkeit bei der 
numerischen Anwendung der Reihe 1aBt sich aus den Beitragen ersehen, die die héheren Nahe- 


rungen liefern. 


In gleicher Weise wird die Funktion F(z) fiir (16) entwickelt. Die Spiegelung der Funktion 
A,,4 


j?/z an den beiden Randern liefert fiir die erste Naherung 


A slg tdas Jee J | gu pb thedl <i 
Es g | BFA, oA GEIAY | (BF 2A) 
Die Spiegelung der in den Punkten x = + 2A, singular werdenden Funktionen fithrt zu 
ofa all (eens i 2 “3 
SIF EP EINECS 2) ae i ee a ee 
hele (G2) (2 ¥ 2)3 (FF 2)8 ; 


dieser Ausdruck liefert mit der ersten Naherung zusammen fiir die gesuchte Funktion die zweite | 
Naherung. Die Fortsetzung des Verfahrens fiithrt schlieBlich zu der folgenden Reihenentwicklung 
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Damit ist (16) vollstandig bestimmt, so daB die Spannungsverteilung in einem Parallelstreifen fiir 
vorgegebene A, , 4,-Werte numerisch ermittelt werden kann. 

Von besonderem Interesse sind die Sonderfalle 4,=/,—0,5 und j,=1,4,—0, die sich leicht 
aus den vorliegenden Ergebnissen ableiten lassen. Auf eine Wiedergabe der Spannungsfunktionen 
fiir diese beiden Belastungsfalle soll hier der Kiirze halber verzichtet werden. Der erste Fall, der 
die Beanspruchung eines Parallelstreifens durch eine in seiner Mitte angreifende Kraft Q betrifft, 
wurde bereits auf andere Weise von E.Melan! und K. Girkmann? behandelt. Die Angaben dieser 
beiden Arbeiten stehen in guter Ubereinstimmung mit den numerischen Ergebnissen, die anhand 
der hyperkomplexen Spannungsfunktion erzielt wurden. Dabei konnte die Reihenentwicklung 
fiir Punkte, die auf einer Parallelen zur x-Achse einen Abstand a/4 (bzw. a) haben, nach der 6. 
(bzw. 12.) Naherung abgebrochen werden. 

Die bisherige Untersuchung behandelte den Spannungszustand in einem Parallelstreifen, 
dessen Belastung durch eine Kraft parallel zu den Randern erfolgte. Der Fall, da8 ein unendlich 
langer Streifen senkrecht zu seinen Randern beansprucht wird, ist praktisch nicht realisierbar 
infolge der unendlich gro® werdenden Biegespannungen o,. Erst der Ubergang zu einem Parallel- 
streifen von endlicher Linge fiihrt zu einer neuen bedeutsamen Fragestellung, namlich zur Be- 
handlung von Scheiben rechteckiger Umrandung, die durch Punktkrafte beansprucht werden. Es 
liegt nahe, auch hierfiir eine Liésung mit Hilfe hyperkomplexer Singularitaten aufzubauen. Auf- 
gaben dieser Art sind bislang mit Hilfe Fourierscher Reihen und der Differenzenrechnung be- 
handelt worden. Daneben ist die Arbeit von J. Fadle? zu nennen, der mit Hilfe komplexer Eigen- 
wertfunktionen den Spannungszustand einer quadratischen Scheibe bestimmte. 

d) Zwei diametral entgegengesetzt wirkende gleich groBe Krafte im Innern 
einer Kreisscheibe. Im folgenden soll der Spannungszustand einer Kreisscheibe vom Radiusa 
bestimmt werden, die in ihrem Innern durch zwei gleiche, entgegengesetzt gerichtete Krafte P 
auf Zug belastet wird. Die beiden Krafte greifen nach Abb.5 auf einem Durchmesser an; ihre 
Angriffspunkte haben vom Kreismittelpunkt 0 die 
Abstande b’ und b’’. Weiter werden die beiden in der 
Abbildung eingezeichneten rechtwinkligen Koordi- 
natensysteme x’, y’ und x’’, y’’ eingefiihrt. Der Ko- 
ordinatenanfangspunkt 0’ (bzw. 0’’) liegt auf der 
reellen Achse und ist vom Kreismittelpunkt b’-+- c’ 
(baw. b’’--c’’) Langeneinheiten entfernt. Die Werte 
fiir c’ und c’’ ergeben sich aus Gleichung (24) und (25). Abb. 5. Belastung im Inner einer Kreisscheibe. 

Fir die zu ermittelnde Spannungsfunktion besteht 

nach (13) die hyperkomplexe Konstante aus ersten und dritten Elementen. Wird zunachst das 
(x’, y’)-System betrachtet, so sind im Punkte x’=c’, y’—0 die beiden Singularitaten 1/(z’— ec’) und 
J?|(# —c’) anzunehmen. Statt der rechtwinkligen Koordination x’, y’ werden nun die Bipolarkoordi- 
naten J’, eingefihrt. Fiir Punkte der Scheibe ist A’ >0. Die beiden Singularitaten miissen an der 
gegebenen Kreiskontur A’) = const. so gespiegelt werden, da die auf das neue Koordinatennetz be- 
zogenen Spannungen Ope, und Ty’ pang verschwinden. Dies erfordert nach dem Spiegelungsver- 


fahren die Annahme eines Poles dritten Grades im Punkte x’-=—c’, y’=0. Die Funktionen, aus 
denen sich dieser Pol zusammensetzt, enthalten wieder Koeffizienten, die nur aus ersten und 
dritten Elementen aufgebaut sind. Zur Bestimmung dieser Koeffizienten werden die Beitrage 
der einzelnen Singularitaten fiir die nach den Transformationsformeln (26) errechneten Span- 
nungen Oped, und Tu! pan nach Potenzen von cosj’ geordnet. Die Koeffizienten sind dann so 
festzulegen, da alle Terme gleicher Potenzen von cos’, die sich aus den beiden Randspannungen 
ergeben, verschwinden. Diese Bedingung ist mit den angegebenen Singularitaten allein nicht zu 
erfiillen. — Ahnliche Verhaltnisse kommen in der komplexen Potentialtheorie vor. Beispielsweise 
ist zur Bestimmung der Stroémungsfunktion fiir das Feld einer Quelle, das durch eine Kreiskontur 
gestért wird, auBer der zu der Quelle spiegelbildlich gelegenen Singularitat im Kreisinnern noch 
eine Singularitat im Mittelpunkt des Kreises anzunehmen, die der im Unendlichen gelegenen 
Senke entspricht ; die Bedingung, daf der Kreis zur Stromlinie wird, 1aBt sich dann erfiillen. Wird 
derselbe Gedanke auf das hier vorliegende elastische Problem angewendet, so ist zu den bereits 
vorhandenen Singularitaten noch im Punkte 0 der Abb.5 ein Pol anzunehmen. Weiter kann eine 


1 E. Melan. Z. angew. Math. Mech. 5 (1925), S. 314. 
2 K. Girkmann, Flachentragwerke, S. 109, 2. Aufl. Wien 1948 oder auch Ing.-Arch. 13 (1942/43), 5.273. 


8 J. Fadle, Ing.-Arch. 11 (1940), S. 125. 
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hyperkomplexe Konstante von der Art const - (1—j?) aufgenommen werden, die den von einem | 
hydrostatischen Druck herrithrenden Spannungszustand darstellt. Mit einem solchen Ansatz ge- 
lingt es, die an der Kreiskontur vorgeschriebenen Randbedingungen zu erfiillen und damit die — 
noch unbekannten Koeffizienten der einzelnen Funktionen festzulegen. Die beiden gegebenen 
Singularitaten 1/(z’—c’) und j?/(z/—c’) fiihren dann zu den Kombinationen 


1+)? ye aya eaey 


1 2f—1 pe, Lj 20! CofA, a 
(z’ + ¢’)3 ) z’ — (b’ +c’) ] 


ere, | 5 ae c (2’ +-¢’)? 
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if P ee eee Lea pee ee Joe 
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re = e” Sind, (1—j?) . 


ee 
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Hierin sind besonders die im Punkte 0 singular werdenden Funktionen zu beachten, deren Vor- 
handensein sich dadurch anschaulich erklart, da®B die Scheibe infolge einer im Punkte x’=c’, 
y’=0 angreifenden Kraft P sich nicht im Gleichgewicht befindet und deshalb einen weiteren 
singularen Punkt im Innern der Kreisumrandung erforderlich macht, an dem eine gleich groBe 
aber entgegengesetzt gerichtete Kraft Panzunehmenist. Fir den Sonderfall A’= 0 folgen aus (19) 
unter Beachtung der verschieden gewahlten Koordinatensysteme die erste und die dritte Zeile 
der Gleichung (15); die im Punkte 0 vorhandenen Singularitaten von (19) riicken fiir diesen Fall 
ins Unendliche, so daB sie fortgelassen werden kénnen. In analoger Weise kann nun von den 
auf das (x’’, y’’)-System bezogenen Funktionen 1/(z’’+-c’’) und j?/(z’’-+ c’’) ausgegangen werden, 
um die zugehérigen Singularitaten zusammenzustellen. Eine Rechnung ist dazu entbehrlich; 
anhand (19) 1a8t sich unmittelbar das Ergebnis angeben. Unter Einfiihrung der Bipolarkoordi- 
naten A’’, ’’ werden in (19) die Verdnderliche z’ durch z’’, die Werte A’ durch —d”, b’ durch — b’’ 
und ¢’ durch —c” ersetzt. In dem hier interessierenden Bereich der Halbebene x’’< 0 ist dabei 
A’<0. Aus diesen Betrachtungen folgt unter Beachtung von (13) fiir das in Abb.5 veranschau- 
lichte Problem die Spannungsfunktion, die aus den auf die beiden Koordinatensysteme bezogenen 
Anteilen zusammengesetzt ist, 
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In (20) wurden beide Koordinatensysteme beibehalten, da so die Ubersicht am ehesten ge- 
wahrt bleiben diirfte. Die nach (19) im Punkte 0 singular werdenden Funktionen heben sich 
gegenseitig auf; sie treten daher in (20) nicht mehr auf. Fiir Punkte des horizontalen Durch- 
messers der Scheibe (Abb.5) lassen sich die Spannungen o, und o, leicht angeben, da der Nenner 
der Spannungsfunktion (20) reell ist; fiir Punkte des vertikalen Durchmessers folgen die Span- 
nungen aus der Zerlegung der einzelnen Funktionen in ihre Elemente. Die beiden Normalspan- 
nungen werden auf beiden Durchmessern zugleich die Hauptspannungen. Die Auswertung des 
Integrals 
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4, Element 
kontrolliert die Gleichgewichtsbedingung. 


Fir den Sonderfall c’—=c” +0 greifen die beiden Krafte am Kreisrand an. Mit A’=—A”=0 
und Beachtung von (24) folgt aus (20) 
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Dies ist die Spannungsfunktion fiir die Spannungsverteilung einer Kreisscheibe vom Radius a, 
die durch zwei gleiche einander gegeniiberliegende Krafte P auf Zug belastet wird. Die auf Druck 
beanspruchte Kreisscheibe wurde bekanntlich erstmalig von H. Hertz! untersucht. Auch fiir den 
speziellen Fall, da eine Kraft im Innern der Scheibe, die andere auf dem gleichen Durchmesser 
aber am Kreisrand angreift, laBt sich die Spannungsverteilung aus (20) ohne weiteres anschreiben. 

Der in Abb.5 dargestellten Belastung einer Kreisscheibe ist der Fall der Beanspruchung einer 
unendlich ausgedehnten Scheibe verwandt, die durch ein Zugpaar gespannt wird, wenn zwischen 
den Angriffspunkten der beiden Zugkrafte eine Bohrung vorhanden ist. Dieser Fall wurde von 
R. Sonntag? unter der Voraussetzung behandelt, da® der Durchmesser der Bohbrung klein ist 
gegeniiber dem Abstand der Kraftangriffspunkte. Die Anwendung der hyperkomplexen Theorie 
auf diese Aufgabe gestattet nach den bisherigen Uberlegungen die Spannungsverteilung in ge- 
schlossener Form anzugeben. Greifen naimlich die beiden Krafte in den Punkten z—=—c’ und 
z—=- e’’ an, so ergibt sich fiir diesen Belastungsfall mit den alten Bezeichnungen der Abb.5 die 
hyperkomplexe Spannungsfunktion 
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Die dabei aufgenommenen Funktionen, die fiir den Mittelpunkt der Bohrung einen Pol zweiten 
Grades liefern, driicken hier den Einflu8 eines hydrostatischen Druckes aus. Fiir den Grenzfall, 
daB die beiden Krafte unmittelbar am Rande der Bohrung auf einem Durchmesser angreifen, 
folgt aus (21) die dieser Belastung entsprechende Spannungsverteilung, die bereits von R. Sonn- 
tag® mitgeteilt wurde. 

' Bei den Untersuchungen dieses Abschnittes stand die Ermittlung der Spannungsverteilung 
im Vordergrunde. Um auch die Verschiebungen zu bestimmen, sind die Gleichungen (12a) und 
(12b) heranzuziehen. Anhand (20) und (21) laBt sich weiter die Spannungsverteilung fiir ein 
kreisringformiges Kettenglied entwickeln, dessen Randkreise konzentrisch oder exzentrisch an- 
geordnet sein kénnen und das am Innen- oder Au®enkreis durch zwei gleich groBe und entgegen- 
gesetzt gerichtete Krafte auf Zug oder Druck belastet wird. Die Lésung dieser Aufgabe fiihrt 
auf eine Reihenentwicklung, die in ahnlicher Weise wie das Ergebnis fiir die Spannungsverteilung 
einer Punktlast in einem Parallelstreifen aufgebaut ist. 


5. Zusammenfassung und Schlufbetrachtung. Ausgehend von einer von L. Sobrero gegebenen 
Darstellung des Zusammenhanges einer hyperkomplexen Algebra mit der Aufgabe der ebenen 
Elastizitatstheorie wurden fiir denebenen Spannungszustand einer Scheibe aus einer hyper- 
komplexen Funktion die allgemeinen Gleichungen fiir die Spannungen und Verschiebungen ab- 
geleitet. Die weitgehende Analogie zur komplexen Potentialtheorie in Verbindung mit einem er- 
weiterten Spiegelungsverfahren erleichtert erheblich das Verstindnis der Zusammenhinge im 
Hyperkomplexen. Die durchgefiithrten Anwendungsbeispiele machen die sich damit bietenden 
Vorteile gegeniiber den in der Einleitung genannten Methoden deutlich, wie ein Vergleich mit 
den aus dem Schrifttum bekannten Lésungen zeigt. Bei den untersuchten Belastungsfallen han- 
delt es sich vorwiegend um Randwertprobleme erster Art. Diese sind dadurch gekennzeichnet, 
daB der Scheibenrand spannungsfrei ist (sofern nicht an Randpunkten Krafte angreifen). Zur 
Ermittlung des Spannungszustandes werden bei der hyperkomplexen Theorie der EinfluB der an 
den Scheiben angreifenden Krafte wie auch die Bedingung der vorgeschriebenen Spannungs- 
freiheit des Randes durch Singularitaéten dargestellt. Mit der Auffindung der hyperkomplexen 
Spannungsfunktion, die eine besonders geschlossene Lisungsform liefert, ist die Aufgabe geldst ; 


1 H. Hertz, Z. Math. Phys. 28 (1883), S.125. ; 
2 R. Sonntag, Mitt. Mech.-Techn. Labor. der TH Miinchen; 3.Folge, Heft 34, 5.37. Miinchen 1930. 
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die Spannungen und Verschiebungen selbst errechnen sich dann auf ganz elementare Weise. 
Infolge des iibersichtlichen Aufbaus der gewonnenen Spannungsfunktion werden etwaige Fehler 
bei der numerischen Auswertung leicht erkannt. 

Die vorliegende Untersuchung 1aBt es wiinschenswert erscheinen, die hyperkomplexe Theorie 
auf eine ganze Reihe weiterer konkreter Einzelprobleme anzusetzen, da dieses Verfahren auch fiir 
schwierigere Aufgaben hei ertraglichem Arbeitsaufwand numerische Ergebnisse von ausreichender 
Genauigkeit verspricht. Vorzugsweise diirften hier neben den Randwertproblemen erster Art 
Randwertprobleme zweiter Art, bei denen die Werte der Verschiebungen am Rande vorgeschrie- 
ben sind, wie auch gemischte Randwertprobleme in Betracht kommen. 

Bei AbschluB der vorliegenden Arbeit machte mich Herr Prof. A. E. Green, University of 
Durham (England), auf verschiedene Arbeiten! aufmerksam, die sich auch mit der Ermittlung 
der Spannungsverteilung in elastischen Scheiben befassen und die nicht unerwahnt bleiben sollen. 
In diesen Aufsaitzen werden die Spannungs- und Verschiebungskomponenten aus zwei komplexen 
Potentialfunktionen abgeleitet. Auch dieses Verfahren fihrt zu eleganten Liésungen und es wird 
interessant sein, die weitere Entwicklung zu verfolgen. 


6. Anhang. a) Uber das Rechnen mit hyperkomplexen Zahlen und Funktionen. 
Um eine Einfiihrung in die hyperkomplexe Algebra zu geben, sollen noch einige Hinweise folgen. 
Dabei ist nicht an die Darstellung einer systematischen Theorie der hyperkomplexen Zahlen und 
Funktionen gedacht. Die folgenden Ausfithrungen sind deshalb in ihrem Umfang nur darauf be- 
schrankt, verschiedene in der Arbeit benutzte Satze sowie einige Besonderheiten zu erlautern. 
Fiir ein weitergehendes Studium muB auf die Originalarbeit von L. Sobrero und das dort zitierte 
Schrifttum verwiesen werden. 

Die Addition bzw. Subtraktion von hyperkomplexen Zahlen erfolgt durch die Addition 
bzw. Subtraktion entsprechender Elemente. Das Produkt zweier hyperkomplexer Zahlen liefert 
im allgemeinen ein Polynom bis zur sechsten Potenz von j; die Potenzen j*, j> und j® lassen sich 
mit Hilfe von (4) und weiter, wenn diese Gleichung mit j und j? multipliziert wird, auf die ersten 
vier Potenzen von j zuriickfiihren, wie es unter Ziffer 3 angedeutet wurde. Die Aufspaltung eines 
Quotienten in seine vier Elemente ist etwas umstandlicher; fiir die praktische Rechnung ist 
die Schreibweise von T. Levi-Civita besonders geeignet und tibersichtlich. Danach 1a8t sich eine 
hyperkomplexe Zahl 


z= a+ jb+ jet jd 


— a, b, c und d sind reelle Zahlen — in die Form umschreiben 
z= (0’ +-ia’’) + (B’ + 7B”) (1 + j?). 


Die GréBe i ist die imaginare Kinheit und, wie schon angefihrt, durch i— ye j+J*) mit den 


hyperkomplexen Einheiten verkniipft. Ferner bestehen die Zusammenhange 
O=a—c; o'=b—d; pe; B’ => (3d—2). 
Die Zahlen «= a’ -+-ia’’ und 6 = f’-+ if” sind dann komplex; die hyperkomplexe Zahl z1aBt 


sich somit in der Form 
s=atpti 
darstellen. Nach Multiplikation mit «—{(1-+ j?) wird der Quotient 
il 1 : 
s=5 804+. 


z a 


Die Erweiterung mit dem konjugiert komplexen Wert von « bzw. o? fithrt nach einigen Umrech- 


nungen zu der in seinen vier Elementen dargestellten Funktion 1/z. Nicht méglich ist die Division, 
wenn 


oF? —= w/2 + //2 — (a — c)® + (bh —d)? 


verschwindet. 


1 N. Muschelisvili, Z. angew. Math. Mech. 13 (1933), S.264; vgl. weiter die hier an b Li 

> 4. ow. 5 Dns 5 teratur, 
K. J. Schulz, Over den spanningstoestand in doorboorde platen. aie! Delft 1941. yh 
A, C. Stevenson, Philos. Mag. (7) 34 (1943), S.766. 
A. C. Stevenson, Proc. Roy. Soc., A, 184 (1945), $.129 u. 218. 
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Fir den Sonderfall a= x, b—y und c—d=0 wird z.B. 


: = : “ | ? amy | yy | 15 & YY" 3 2 
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mit r°— x*+ y*. Steht statt zder Wert z—2z, mit z= x%)+j V9, so ist x durch x— x, und y durch 
Y— Yo ZU ersetzen. 

Fiir die durchgefithrten Anwendungen ist die Behandlung der Funktion 1/z" mit reellem n 
bedeutsam. Dafiir wird mit r? = x?-+ 2 
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Die Anwendung des binomischen Satzes fiihrt zu der Reihenentwicklung 
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Mit der Beziehung i (1-+-j?)—j +-j® und des bekannten Wertes fiir i liefert diese Gleichung fiir 
ungerades n nach Aufspaltung in ihre vier Elemente die Darstellung 
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und entsprechend fiir gerades n 
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Dabei sind die Potenzen von x durch solche von y ausgedriickt worden, um das Rechnen fir 
eine numerische Auswertung der Formeln mit Hilfe der Rechenmaschine weitgehendst zu ver- 
einfachen. 

Fir spezielle Anwendungen kann es niitzlich sein, einen entsprechenden Ausdruck fir 1/z” in 
Polarkoordinaten darzustellen. Mit der Einfiihrung von ?=arctg y/x unter Beibehaltung der 
alten Bedeutung fiir r folgt 
seats LD ging 4 nt (14 2) e- t+ gin 
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und zerlegt 


ones 
== cosn§ + + sind sin(n | no els sin P cos(n-+1)) — J sin n9| + 


gn 


; ; Tare 
j= = sind sin(n 1) a- oe sin) cos(n + 1) — = sin nv 


Von den hyperkomplexen Funktionen ist weiter der Logarithmus einer hyperkomplexen | 
Variablen bedeutsam. Mit der Levi-Civitaschen Schreibweise und unter Anwendung der Taylor- | 


schen Formel wird 


nz =In|(w $y) ~ 5 +H) 9] =In(@ +39) — 5 TF) ep 


x+iy’ 


Nach Aufspaltung der komplexen Ausdriicke folgt schlieBlich die Gleichung 


co 2 
lng = sin (x? | y2) iaretg 2 : (1-97) ee 5 (i577) - (23) | 


x2 y? xe ty?’ 
aus der sich unmittelbar die vier Elemente ergeben. 

b) Zur Geometrie der Bipolarkoordinaten. Der Zusammenhang zwischen recht- 
winkligen Koordinaten und Bipolarkoordinaten la8t sich unmittelbar aus einer hyperkomplexen 
Transformation in Anlehnung an die in der komplexen Theorie geltenden Beziehungen’ ableiten. 
Mit der Apolloniusschen Kreisschar A=const. und der dazu orthogonalen Kreisschar « = Const. 

folgt 
a. on oo! Sind : sin | 
oat — © €oj At cos wu | J° GofAFcos ph” 


Dabei bedeutet der Wert c einen Parameter, der das vorliegende Bipolarkoordinatennetz be- 
stimmt. Mit den Bezeichnungen der Abb.6 ist weiter fiir alle Punkte eines Kreises der Schar 
A = const. das Verhdltnis 
0./0, = e’ und fiir Punkte 
der Schar wu = Const. der 
Winkel wu = x — 0, + O,. 
Die Trennung der Elemente 
in der Gleichung fiir z 
liefert 
ons 
Sind sing” 

Weiter folgt nach einigen 
--£ | Zwischenrechnungen 


Se 
Ss 


See | yl 
: aa Raa 
Abb. 6. Bipolare Koordinaten. NG oy 
tg = 


Foy 
Diese beiden Gleichungen und die Transformationsgleichung fiir z gestatten, fiir jeden Punkt der 
Ebene die Koordinaten des einen Systems durch die des anderen auszudriicken und umgekehrt. 
Fiir die Kreise A—const. und 4s = Const. folgen die Gleichungen 
BEA olin ee oe i 

(x c &tg A) miro ~ SinrA ’ ONY eG) a 
Die Mittelpunkte der Kreise A= const. liegen auf der x-Achs i i i 

. - e und sind durch die Ab 
x—=cetgA bestimmt. Die Radien dieser Kreise betragen ua aed me 


es 
"— Cina (9) 

Der Abstand } ermittelt sich aus b+c—ctgd zu 
Diem (25) 


1 Vel. etwa W. Fliigge, Strenge Berechnune i i i 
Aig Hes mas i fore anes Ce ee unter Einzellasten, $.3, Berlin 1928 
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Die Gleichung fiir «= Const. liefert differentiiert 


dy ‘pee Cex ees Gin A sin uw 
dx SY ytectgu = 1+ CofAcosy 
und weiter 
pane 1+ Go} Acos Sin A sin 1 


GojA-+ cos wu ” CS A Soret cose 


wenn mit y der durch die Tangente an die -Kurven mit der x-Achse gebildete Neigungswinkel 
im betrachteten Punkt der Ebene bezeichnet wird. Mit diesen Ausdriicken lassen sich dann die 


bekannten Transformationsgleichungen 


0, =0,cos*y + o,sin*y + 7sin2y | 


O, == 6,sin’y + o,cos*y — rsin2y 


Tay —= (y— Gx) siny cosy + tceos2y 


(26) 


benutzen, die die in rechtwinkligen Koordinaten ausgedriickten Spannungen o,, oy und 7 in die 


auf Bipolarkoordinaten bezogenen Spannungen umzurechnen gestatten. 
Umrechnungen sei noch auf den Zusammenhang 


(CojA + cos)? = (1 + CofA cosm)? + Sin?A sin? w 


hingewiesen, der mit der Gleichung sin?y + cos? y = 1 identisch ist. 


[Aus dem Festigkeitslaboratorium der Technischen Hochschule Hannover.] 


(Eingegangen am 15.Marz 1951.) 


Anschrift des Verfassers: Dr.-Ing. Kurt Schmidt, (20a) Hannover, Technische Hochschule. 
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Beitrige zur Theorie der achsensymmetrisch belasteten dicken 
Kreisplatte insbesondere bei elastischer Lagerung. 


Von I, Szabo. 


1. Einleitung und Problemstellung. In einer fritheren Mitteilung! habe ich das Problem der 
elastisch gebetteten dicken Kreisplatte unter Vernachlassigung der zwischen Platte und Unter- 
lage auftretenden Reibungskrafte behandelt. Der Zweck der folgenden Ausfiihrungen ist es 
einerseits, dasselbe Problem bei Beriicksichtigung der letztgenannten Schubkrafte zu lésen, 
andererseits auf weitere Anwendungsméglichkeiten der hier angewandten Methode hinzuweisen. 
Eine spater folgende Mitteilung wird die praktische Auswertung der Ergebnisse bringen. 

Es sei co=f(r) die achsensymmetrische Last einer Platte von dem Halbmesser a, der Hohe h, 
dem Schubmodul G und der Querkontraktionszahl y (Abb.1). Die Platte ruhe auf dem Halb- 
raum als Unterlage und fiir das elastische Verhalten beider seien die Hookeschen Gesetze mab- 
gebend. Dann bestehen zwischen den Verschiebungen ¢ und g (Abb.1) bzw. der aus ihnen ab- 
leitbaren Gesamtdehnung ¢ und den Spannungen o,, 0,, 0,, T (Abb.2) folgende Beziehungen: 


ee as 

Saal rg AOE (1) 
= OG VE a de | VE 

5. 26(5: i) a 26(a, | tale 


(2) 


= 0 VE ey ess 
Me =26(- T=) t= 6+ 55) | 


und dieselben Gleichungen fiir die Unterlage, wobei fiir diese alle GréBen mit dem Index | ver- 


sehen sind (Abb. 1). 


Abb. 1. Rotationssymmetrisch belastete Kreisplatte auf Abb. 2. Volumenelement der Platte bzw. Unterlage. 
elastischem Halbraum. 


Unter der Annahme, da die zwischen Platte und Halbraum auftretenden Reibungskrafte 
als Schubspannungen aufgenommen werden, lauten die Randbedingungen, wie folgt: 


12, 6) —— 0, (3) 

0,(0, r)= —fi(r),— (4) 

T(O;-r) ==0> (5) 
jo.(h,r) fir rsa, 

1.00.7) = Ig oe (6) 

C1 (0, r) = C(h, r) eae (7) 


1 J, Szabé, Ing.-Arch, 19 (1951) S. 128, Druckfehlerberichtigung hierzu S. 354. 
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|r (h, r) funeers << ae 


7, (0, r md eee, (8) 
01 (0, r) = (h, r) 3 r <a > (9) 
On (z,. 4) == 0. (10) 


Die Randbedingung (8) bedeutet, daB zwischen Platte und Unterlage beliebig groBe Reibungs- 
krafte auftreten kénnen (also ,,unendlich rauhe Berithrungsflaichen“’), wahrend der ersten Mit- 
teilung vollkommene Glattheit zugrunde gelegt wurde; zwischen diesen beiden Grenzfallen 
sind offenbar alle méglichen Faille eingeschlossen. 

Die Randbedingung (9) kénnte man auch durch die physikalisch gut vertretbare Forderung 


x (htt) = fu. (hr) (9a) 
(fo= Haftreibungszahl) ersetzen; die Erledigung dieses Falles bereitet keine neuen Schwierig- 
keiten und er wird hier in kurzen Andeutungen mitbehandelt. 
Die Bedingung (10) bezieht sich auf die freigelagerte Platte; fiir die eingespannte wire sie 
durch 
0 (2,4) = 0 (10a) 


zu ersetzen. 


2. Die Verschiebungen in der Platte und Realisierung der ersten drei Randbedingungen. Aus 
(1) und (2) und den Gleichgewichtsbedingungen 

Oor , Or—or , OT 00: 

or r Oz 


Se ee 0 


ergeben sich! unter Vernachlassigung des Eigengewichtes (dessen Beriicksichtigung naherungs- 
weise in der Weise geschehen kann, dafi man das Eigengewicht als gleichmaBige Last zu f(r) 
schlagt) folgende Differentialgleichungen 


(1—2»)4¢+22—0, (1—29)(40—4) +3 =0, AAC =0, 


wobei A= = ae + den Laplaceschen Operator bedeutet. Der Randbedingung (3) ge- 
niigende Lésungen sind nach Szabd?[Gleichungen (32) und (33)! 
pee Pe ess ream Bata we**)| Jy (ter) +-Az+B, (11) 
Eile eaters: Be (3 —49) a ra 
2 | Pe Oar le Qs sq sma, |° 


(12) 
—qe*)| J, (der). 
Hierbei sind P,, Q;, Px. q, 4 und B noch frei verfiigbare Konstanten, Jy(A,r) und J; (A,r) die 
Besselschen Funktionen erster Art und schlieBlich sind die A, vermége der Beziehung J, (A,a) =9 
durch die unendlich vielen positiven Nullstellen t, von J, (t) festgelegt: 2, = t,/a. 

Die Randbedingung (4) liefert mit (2), (11) und (12) 


26) > ear x AP, An Qe) Jo (te) | t=, | Sear 


k=1 


z 


oe: —Ay, 2 
2(1— 2») (Pre 


Diese Gleichung mit rdr multipliziert und zwischen 0 und a integriert, ergibt wegen 


a 


AEG rar 0 


den Wert von A: 


f(r)rdr. (13) 


1 A, Nadai, Elastische Platien, $.310. 
2 Siehe Funote auf S. 342. 
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Die Multiplikation mit J (A,r) rdr und Integration liefert wegen der Orthogonalitat der 
Besselschen Funktionen 


5 Se (r) Jo (A,r) r dr 
k, 0 a 
SS Pt a? G J? (A, a) Ch (14) 
Die Forderung (5) ergibt 
Ax( Pr Qe) 4 = q) = 9. (15) 
Durch (14) und (15) kénnen P, und Q, auf p, und q, zuriickgefithrt werden: 
es (3— 47) p,, . Cy, 
R= 4 (=o, eae, ae 
a , (16) 
0,=—. Pr Le ct (Ga i) Wh ae ae. | 


4(1— 29) A, 4 (1— 2) A, 2A, 


3. Die Verschiebungen im Halbraum und die Erfiillung der Ubergangsbedingungen zwischen 
Platte und Unterlage. Mit den noch frei wahlbaren Funktionen ®(A) und Y(A) sind die Verschie- 
bungen des Halbraumes C, und 9, nach Szab51 durch die uneigentlichen Integrale 


=i upd D(A) Jp(Ar) e241 da + [ WA) Jo(Ar) e241 dA, (17) 


0 


eorieed | i ie 3—4,)D(A)| _ 
a= saan) | CMA Are* da + [ HA) ee 5 J,(aran | 08) 
0 0 


gegeben. Aus der Randbedingung (6) flieBt nach Elimination von P, und Q), vermége (16) 


ce ie $y BA Arya = 
0 
=o ae | “ > pots py (Ae Wem Sin A, h) — 20 eae (Ay hee” — (19) 


Pireaians Jo (tur) fie re a 


0 ia iP SPUho 
Vermége der Hankelschen Umkehrformel fiir « =0, 1, 2, . 


/ F(A) J,(Ar)AdA = g(r) F(a) = [ p(t) J, (Ar)rdr (20) 
und der Integrale 
i Jo(ar)rdr= * J,(aa), fr Jo har) Jo (Ax) dr = *-Fo@ee) Ji(ha) 
0 eek 
erhalten wir aus (19) 
1 
P (a) + 3 O(a) = — 1 ady( ( ) 4 ee Es arteries é (21) 


wobei zur Abkiirzung 


Me fini : q 
C= 2(1—2») (A, he “*" — Gin A, h) — eee het = Gin A, h) + ¢,Coj Ah (22) 
und 
raeG. 
be (23) 


eingefiihrt wurden. 


1 Siehe FuBnote auf S, 342. 
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Die Randbedingung (7) fordert das Bestehen der Gleichung 


[¥@sanda— Siar [(8 —4y + 2A,h)e %*" 4 e%#™) 4 


4, ae ; 
| Teor [(3 —4— 2a, h)e** 4+ 6%") 4 26, Gin Ay | Ser ney eg: 


woraus sich nach Multiplikation mit rdr bzw. J(A,r)rdr und Integration von r—0 bis ra die 
Beziehungen 


P(A) J,(Aa) = 


} 7 = (Ah + B) > (24) 


bzw. 


(ABS, add — eJy(h0) 
i #B—fe CA Dn (25) 


ergeben, wobei 


CS Ss [(3—49 + sicuaneticd ee | 
a (26) 
to zy hae 4y —2),h)e™" + e”*] + 26, Gin dah | 
bedeutet. Aus (8) ergibt sich mit der Abkiirzung 
Ip — 
Mreieas ehh De jel | 
(27) 


+5 ay [2 Ah + Ie —e*"] — 26, Sindh | 
die Beziehung 
Ap Jay fii ee 
—2{ [va ee ne r)dA = “> k wlAur) ir or a 
: ° ra. 


Wenden wir auf diese Beziehung die Hankelsche Umkehrformel (20) — fiir x1 — an, so er- 
halten wir 


d=7)90) 2. wN r med oe oe 
Oe (1—2»,) i oe all r Jy (Ax 7) Jy (Ar) nara FD Ah . (29) 
Setzen wir den hieraus entnommenen Wert fiir ®(A) in (21) ein, so bekommen wir 
pee dh Ay, @ 
P(A) = + a(l— 2) A ) Ak ee? as | 


f (30) 
> (A, a —v) AJ, (Aa 
2u(1 —%,)a adh (Aa) eee c 7 i ) 5 | 


Diese Gleichung mit J, (Aa) ad multipliziert und zwischen 0 und ©0 integriert, ergibt mit Riick- 


A 

sicht auf (24) 

( di __ a Le [a TEMS 

| P@AG0 = = (Ah+ B) =F a(1—2n) dAdo (A; a) i(e—1) 

; : (31) 

Ie 7 , (—»)A [Ji (a) )da | 
S27 a Ds Gi. Jy (ya) | ae Aa ae | ay | 
as: eA Jala dir =p a 

Die Gleichung (30) mit — multipliziert und zwischen A=0 bis /= 00 integriert, liefert 


2} 
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unter Heranziehung von (25) 


To(ana Mu . PF  AsR(Aa) da 
ee) D,=+¢ 2) >) Avda Jo(toa) | ay a8 


0 


(32) || 


oS st RDA aie (1—»)A “Ja ADE: 
—2n—a)| > Ci, Jo (Ax @) | (22— Ap) (2 ; —72) | 129" aA 
k=1 4 “0 


Denkt man sich hier D,, A, und C, (gemaB (26), (27) und (22)) eingesetzt, so ergibt sich fiir 
n—1,2,3,... ein unendliches, lineares Gleichungssystem fiir die noch unbekannten Koeffizienten | 
Pi: Po» »+- und q,, q, +--+, zu dessen Auflésung offenbar noch ein lineares System benétigt wird 
und das muB aus (9) folgen. Zunichst liefert diese Forderung 


co 


-—_—_—_———-— —— 


| (8 —49,) P(A) ane 
[|v | etna 
sid) Saleen (33) 
=S'lfp+! aD sie 39 Pile NE 
(3 —4?) h e 
OQ, 2 = eke 2(1—2) (p,e pe —9,¢*)| Jy (Avr) 9 (r <a). 
Setzt man hier nach (16) P, und Q, ein, so ergibt sich 
| ay ! ar J, (Ar) a= > B,J, (ur), (r<4), (34) 
wobei zur Abkiirzung 
B= rae a 2 ae etl | 


(35) 


Terre L(3 
gesetzt wurde. Die Gleichung (34) mit J,(A,r)rdr multipliziert und tiber den Giltigkeitsbereich 
(r=0... a) integriert, liefert 


ie —40,) B(A)| Ana J, (Ana Aa) di a 
[ |” ) Seer nee ea = Be Tina) - ae) 


dy 2 Ah) lh 0%") Gof ah | 
k 


Setzen wir hier den aus (29) entnommenen Wert fiir (A) ein, so bekommen wir 


Ji(aa)dA a Bn Jg (Ana) (1—»,) (3 —4,) ap rs SR (Aa) dA 
[vats — jz (rare Sa > Ache (As @) Tap aay a) 


0 
Nun wird die Gleichung (30) mit 
man mit (37) 


Jo (An @)Bn [  AJ3(Aa) da 
(Ee 2n| S Ci Jy (Ae 4) | (A? — 23) (22— 22) ! 
Meo erneny yo lon, ie en 
. ae Paes Dy Adee Aa i, aay | 
0 


das zweite lineare Gleichungssystem fiir die Px baw. q gewinnt. Benutzen wir anstatt (9) die 
Randbedingung (9a), so erhalten wir 


Dia Jat) 2 tae See 


we ry (A, he — Gin Ah) + e oj oe a 


Und diese Gleichung mit rJ, (A,r) dr bzw. r dr multipliziert und zwischen r—0 und r— aintegriert, 


dh aa 


rT pa cae von A= 0 bis A = © integriert, woraus 


(38) 


eee (Ape ‘e® Sind, h) — 


| (33a) 
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liefert das in diesem Falle an Stelle von (38) tretende Gleichungssystem 


Ay, frJy (A,r) JolAnt) oh == a’ fy sates (An hem es Sind, h) 


k=1 hie) (38a) 
Gn 
son) (A, hen — Git Ay h) + ¢, Go) A,h| J2 (A, a) . 
esti) 
> fa USAC Torey tl 
Nach Auflésung von (32) und (38) — woriiber noch zu sprechen sein wird — and auch die 


anderen Koeffizienten bestimmbar: P, und Q, aus (16), Baus (24), nachdem man Y(A) aus (30) 
berechnet hatte. Zum SchluB kann @(A) aus (29) ermittelt werden. Damit ist aber der Verfor- 
mungs- und Spannungszustand in der Platte und Unterlage festgelegt, ohne daB man die Rand- 
bedingung (10) erfiillt hatte; also derselbe Umstand wie in Szabé!. (Fiir die fest eingespannte Platte 
mit der an Stelle von (10) tretenden Randbedingung (10a) ware dagegen das Problem schon gelést, 
da (12) der Forderung (10a) geniigt.) Die Erfiillung der lokalen Spannungsfreiheit des Platten- 
mantels fiir die freigelagerte Platte ist auch hier mit den verwendeten Funktionen nur im Sinne 
des de Saint-Venantschen Prinzips méglich und dieser Frage wollen wir uns jetzt zuwenden. 

4, Radialspannung und Radialmoment am Plattenmantel. Da die Randbedingung (10) 
— o,(a, 2)=0 — nicht erfiillt ist, herrscht am Plattenmantel eine mittlere Radialspannung, 
die sich aus (2), (11), (12), (19) und unter Beriicksichtigung von (27) zu 


guards 1 fg, (aa o| > Ay Zo as os (39) 
errechnen la®t, wahrend san das Radialmoment 
m= fate “edz = Cn | Daeg Os + eg SO deh — spe Sin ah) 
—— bes qep OO ht a Sin A n) iz (40) 
Ee (a Sin Es (1 — Gof Az m)) J, (ce) | 
erhalt. Am Plattenmantel haben wir also das aquivalente Atel 
sd Pia a de 
ig? ees ey (A, he 1 Sin Ah) — Gin At i a = To , ms 
fig = me— 3, = ow > Aes (ti? o** 4.24, Col Ah — Ach Gin Ah — 
—2 Gin ark) + re: vy (Ai 2 Ag h Co Ah — Ane Gi hy (42) 


A 

4 26indA,h) — are a A,h + 2(1 — Coj A.h)] Bike 

5. Die Erfiillung der Randbedingung am Plattenmantel im Sinne des de Saint-Venantschen- 
Prinzips. Das durch (41) und (42) gegebene System der Resultierenden der .,Restspannungen™ 

in radialer Richtung kann durch den gleichmafigen Zug (oder Druck) — o, und durch die reine 


Biegung — img zum Verschwinden gebracht werden, wodurch der Randbedingung (10) ,,im 
Gesamten‘* — wenn auch nicht lokal — geniigt wird. Die zu den Lasten gehérigen Defor- 
mationen — bzw. Verschiebungen — sind nach Szab6' 
,_ tt fh ner, 
PS a ele OY el GR | 
F : 2 / h A 
3m |20 (2— s) +(—1)r| -- 6mp(1—»)|2— s|r (43) 
hae (1-10) Gh8 oh a. (1+) Gh | 


1 Siehe FuBnote auf S. 342. 
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und somit ergeben sich nach (2) die Spannungen zu 


Ce ’ =0, O.=6,=— G6, | 
ZA N WY wv 12m h 44, 
op — We if == 0. oO, = 0; => Bale r)-| ( ) 


Durch diesen zusitzlichen Spannungs- und Deformationszustand werden die aus den Rand- 
bedingungen (7) und (9) flieBenden Formeln gedndert. An Stelle der aus (7) folgenden Glei- 
chung (24) tritt jetzt 


/ (2) J, (aa) = Ah | B- 


wahrend die zur Berechnung von B dienende Beziehung (31) durch 


Gr hy 3mR[vh? + a?(1—»)] | 
2G(1 +7) F 2G(1+)h  [’” ce, 


a | vorh 3mr[vh?+a(1—»)]) __ 
a ek nal x60) 26 oR | 
= ie = " Ji(aa)da 
= (1 20) DS) AedeJo(dea) i (2 — ay i) (46) 
2 ha 1p AL adi (eda 
—2y0t—n)] Ci, Jo (Ax @) f See i | vies | 


zu ersetzen ist. Als letzte der aus (7) flieBenden Gleichungen ist (32) in 


Je Oia) Dem One (I= 9)) 2)  Aye(Aa) dd 
4A, G(1 +r) h322 ! > 2, Ahi So(Ana) | (B— ky (22 Tay 
j (47) 


>, CuJo(Ax a) | OL oa Hoe 
k=1 : 


YO ih Wee 
ey =—jgy(P—a) | 1—20 Ae 


0 


umzuandern. An Stelle der sich aus (9) ergebenden Gleichung (34) tritt 


Way Sa ie (1—3, 3p (1—»)r 
{| geet este 2, Besa (ter) SC(le) i) Clo 


womit man genau so wie mit (34) verfahrt und schlieBlich 2u der an Stelle von (38) tretenden 
Gleichung 


Jo (A,,2)B,, (Gh ome laa [ F(a) di 


(1—2)/, | G(x) (1—2%) 22 h? 1—2y | A(a?—22) | 
Neel out is oe fx “_ Ji(do)da te? 
— k; Fall i @ | a Eee [=o = Ax hy, Jo( (Aya “| (22— Az) al 


kommt. 


6. Das System der unendlich vielen linearen Gleichungen fiir die Koeffizienten p, und q,. 
Dieses System ist durch die Gleichungen (47) und (49) gegeben und wir wenden uns der Frage 


sas ene zu. Die in den Gleichungen (47) und (49) vorkommenden uneigentlichen 
ntegrale . 


fede dh P Jz (1) dt 
EET terme Saat | ere 


le 92.32 (A.a) dd #2J2(t) dt 
(2? 9 = Glee 


— Hi) (PR) | (Pt) (@ 8) 


0 
JiGad, se el), [oe aN Le ( 1J3(t) dt a 
t eS 
( 


ae) ?—t2) (? —t?) = Ones , (50) 


Me) Fe 


Ji We 
[exe =p) 


a? == a 3 1?) 


SSS 
2 


— Ai) (B= 
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erweisen sich — wie man durch Heranziehung der Potenz- bzw. asymptotischen Reihen leicht 
zeigen kann — als existent und I, bzw. I,,,, sind in Szabé1 schon berechnet worden. Auch fiir dje 
neu hinzugekommenen 7, if) und 1 lassen sich ahnliche Formeln wie fiir I, und I;,, angeben, 
sie erweisen sich aber bei der numerischen Berechnung der Integrale gegeniiber der Simpson- 
schen Methode so unterlegen, da8 auf die Mitteilung dieser Formeln hier verzichtet werden kann. 


Die numerischen Auswertungen ergaben folgende Resultate?: 


n i 2 3 4 

I”) | —0,03405 | —0,01018 | —o0,00478 | —0,00282 
a 0,034.02 0,00001 | —0,00002 | —0,00003 
i 0.00001 0,01015 | —0,00007 | —0,00001 
1 | —o,00002 | —0,00007 0,00482 | —0,00001 
1) | —0,00003 | —0,00001 | —0,00001 0,00282 
7 0.01125 0,00083 0.00041 0,00024. 
ie 0.00083 0,00172 0,00014 0,00008 
ih 0.00041 0.00014 0,00054 0.00005 
rs 0.00024 0,00008 0.00005 0,00024 


Setzen wir in (47) und (49) fiir 4,, B,, C,, D,,¢,und mg die Gleichungen (27), (35), (22), (26), 
(41) und (42) ein, so erhalten wir nach Einfihrung der Abkiirzungen 


N™ (,) thes 


4 An(1— 21) 


Jo(An@) re 9 v 


NO (An) = 5 72(1—2») 


F® (4,) = 


(1+ v)h8 72 22 


or + el or 2) (A,he *#" Es Sind, h) A, aI”) - wl —,) (A, he 7#* Es Sindh) aly, ; 


2(1—yv) 
(1+ v) h?22 AP 


FO) (A,) = 


+1 (L— 2m) (Abe #* — Sin Ach) a2 12 — 2 (1—m) (eho *™-+ Gin Ach) Aa TY 


das Gleichungssystem 


(1—2v—A/,h)e’n" + CofA,h 


2(1—2)Gind,h 


Jon ® [1 —2y-+-A, hye ™*+ Col Aah], 


Anh) e+ Coj Anh], 


3(1—») (i he emt ah [oe *** + Cf ah] — 2Gind,h) + 


[2 Azh2e **" + A,h (eR 2 Co} A, h) — 3 Sind, h] + 


Pn 


Jo (Ak a) 


4 (1—2+A,h)e~7n*4 Cof A,h 


ae Pk NO (An)(1— 2 


4-2 (1 — Cof A, h)] — w (1 — 2%) Ana? Ibn Sin Ah — 2 (1 — 4) aT, Coj Anh 


2 u(1—r,) (1—») 


v) 


0, 


(1—2r) NM (A,) 


1 Siehe FuSnote auf S.342. 


nm 


Fy” (Ax) + > Ck 
k= 1 


Jone | 


6(1—») 


(28 pAh) ene Col ALk 


N®) (A,) | (1 +») h? 22 Ap 


[Ax h Gin At h + 


(51) 


(52) 


2 Die Mittel, zur Durchfiihrung dieser und noch weiteren Aufgaben, wortiber in einer — Georee oe 
wahnten — Mitteilung berichtet werden wird, sind mir von der Notgemeinschaft und den Freunden der 
TU Berlin-Charlottenburg dankenswerterweise zur Verfiigung gestellt worden. 


] 
| 
| 
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(129 + Anh) e’n* + Cof Anh Y, 2(1—27)GofA,h 
% (1—20—A,h)e*n* 4+ CofA,h |" (1—2v—A,h)e—*9* 4 CofA,h 


Joh) pO ay 4 Jo49) _ pl _ yy 
De Pk NG, \(1—2) F, (Ax) Dy "Uk (—2») NO(A,) n ( t) i 


53) 
SoG 8) | 2 emai el nr ann 
+>) 4 * nN@(A,) | Se vy he AD AR a Sy eae Ades 


4 2 (1— 2) a2 TL? Cof Ah + 4 (1 — ) 2 A, I Sin Ach | ae 


Pn 


2 u(1—2,)(1—») Wiad 2(1—»)vA —_ 
(1—2») NP) (A,) "(+2942 NOG,) 
Durch die Transformation 
Pr=2eq,(1—2r)+6,, =e (54) 
ergibt sich aus (52) und (53) 


(1—2v—/,h)e’m*®+4+CofA,h Ds | ; GinA,h ith. 
Bos (1—2y+A,h,)e 79" 4+ CofA, h pe se (1—2y+ A, h)e7n* +4 CofA,h 
e Jon) Fo Jo (Ay a) (1) | 
0 ry jee rh) 
re ‘aan + Dra Sanita? © (55) 
oa 
+> Log (2 Fs + ae yaaaga Aah Sin Aah-+ 2 (1 — Co} A) — 


(1) |__ 2u—)(1—») 1a 
—p (1 — 2) aA. Tin SinA,h—2u (1— Mi) Genel Aer, aA NOCG Aal==0; 


A, he 
ere (1—2v+A,h)e jacked h Barat 2») eos ee Tee 
(1—2v—A,h)e—*»*#4 oj Ah M27) hem eo ne 4 
oa 3 (A, 4) Jo (i, 2) (2) 
re A ey ARO) (a2 )y2 
ie ser Te a) + Daz SNOT es oe 


Jo(4x@) { 2(1— ) 


0%) Jo pl) = (56) 
we NMG OP eats Cee rues es -Cof Ae h)) 


+2 (12) af) Coj Agh-+ 4 (1 — »,) AeaP IL) Sin Ay h| rE 


L 2 uw (1—2»,)(1—») a AT) . 2(1— vr) vA 
(1— 20) N(A,) "  (+9) = 29)An NO An) 
Beziiglich der Auflésbarkeit der Gleichungen (52) und (53) 1a8t sich aus dem durch die Trans- 
formation (54) entstandenen System (55) und (56) folgendes aussagen: Sowohl in (55) wie in (56) 
gehen die Konstanten (d. h. mit 6, und ¢, nicht behafteten) Glieder von einem gewissen Zeilen- 
index N an gegen Null. Von diesem Zeilenindex an gehen fiir beide Systeme die Koeffizienten 
der 0, und ¢, bis zu einem gewissen Spaltenindex K < N gegen Null; folglich erhalt man von 


der N-ten Zeile ab ein homogenes Gleichungssystem fiir 6x,dx41,... und éx,€x41,... welches 
durch 


On = Oxy = = ex EK+1 Sasa 0 


gelést werden kann, d.h. fiir groBes k kann 

Pi=2 (1 — 29) 4, “Gg. 0 (57) 
gesetzt werden, wahrend die ersten Koeffizienten p, und q, so gewonnen werden, daB man aus 
(52) und (53) in der itblichen Weise ein endliches, lineares Gleichungssystem herleitet. Im An- 
schlu8 an das vorher Gesagte bedeutet das, dafs wir aus den nach Abspaltung der homogenen 
Gleichungen iibrig bleibenden 2N inhomogenen Gleichungen mit 2 K(N > K) Unbekannten 
die ersten 2 K Gleichungen mit den ersten 2 K Unbekannten herausnehmen und auflisen. 
Dieses Vorgehen ist, was die ¢ anbelangt, um so mehr erlaubt als doch jede Belastung f(r) 
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durch eine endliche Besselreihe (im wesentlichen von der Form Sey Jo(dut) beliebig genau approxi- 
k=1 


miert werden kann, was das Verschwinden der c, groBer Nummer nach sich zieht. Da weiter 
DN? (,) und 1/N(A,) fiir groBe n sehr stark abnehmen, ist auch ihr Ersetzen durch Null er- 
laubt. Es bleiben dann endlich viele Gleichungen fiir endlich viele Unbekannte 6,, &, ibrig. 
Das Nichtverschwinden ihrer Determinante, das ihre Lisbarkeit sichert, erscheint dadurch ge- 
wahrleistet, daB die zugehérigen homogenen Gleichungen (A—0) nach (13) der unbelasteten 
Platte (d.h. verschwindenden Verschiebungen) entsprechen. Daher sind dann p,—q,—0 die 
einzigen Lésungen der homogenen Gleichungen. Diese Uberlegungen entheben uns der Mihe, 
der mathematisch schwierigen Frage nach den Liésungsbedingungen fiir das System der unendlich 
vielen linearen Gleichungen (55) und (56) mit unendlich vielen Unbekannten nachzugehen. 
Mit dem Kroneckerschen Symbol 

[OM vktee ah, 

_ 58 
Onk \1 es k ==, ( ) 


1aBt sich das Gleichungssystem (52), (53) auf die Normalform 

co fo.) \ 

>, Onk Pk + Bs papte = Vx 

ie se Nn seias) (59) 

> Gn k Pk + ») bat, = dn | 
k=1 k=1 


bringen, wobei on%, Prk» Yn> %,,> Ong und d, folgende Bedeutung haben: 
Pree Jo(An a) I —A,h h 
On = a —3») 7, 2v+A,h) e + CofA, A] Onn + 


Jo%,a) { 3(1—») ie ne Sauk eS > 
} ne eee [ante **® + a,h(e**L Coj Aph) —2 Gin Ah] 


—£(1— 2m) (Ake + Sin Anh) Aya? Ike + (1—») (A, he *" — Sin Ah) Teal 


Prk = ee [((l — 2 »—A,h) en® + Gof Anh] Ong + 


Jo(Ax 2) { 3(1—») 
1—2» \(l+»)An Ay 


—£ (L— 2m) (A he*" + Gin Ah) A, a 1 — np (1—») (A, he — Sin Ah) ain! 


va [Azhee*™ — A, h (eo + Coj Ah) + 2 Sin Agh] — 


eee 2 ee) et 


| 
| 
| 
| 
t= — lt Fie Jol, 4) en + | 
| 
| 
| 
| 
| 


$2 >? Solda) ae ge apps * Ua Gin Aah} 2 (1 — Gof deh) — 
k=1 ie 


—£ (1-2) Gin Ah- Aa? tin — (1 — %) Cofaah- ali Ck» 


Pot = d= [(1—2 »—Ayh) 0 """ + Cof Anh] Ona + 


pF) (oars ARM et Auk (M+ 2.6of Ask) — 8 Gin Aah] — 
—ay nk 


SSA re e** _ Sin Ah) a? I) + ps (1—%) (4, he" + Gin Ay h) Aa? In ? 
2 


Son) fie Srl ehlO.. 6 
bak = aq 9) 73 LUE 2y-+A,h) en4 + Coj Anh] Onn + 


well (A a) 1—» 242 Ah h Ay, h 2p j, h 3 Sind h a 
+ Sate Ceres [2 Azhe Ayh (e**" +. 2€of A,h) + | 


Aye . (2)| 
+4 (1-2) (ap he** — Sin Aph) a? n+ je (L—m) (eh e*” + Sitt Aeh) Ana? Tne 
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2(1—r) vA 
~ (1+) (1—29) 22 


= = a : =a) A a2 1 (1) vee Hae a) Cn ¢, C0} A,h 


+ > hima)! re ee 5 (2 deh Sin Agh + 3 (1 — Coj dehy] + (64) 
k=1 
424 (1—2) a HY - CofAgh + 4 (1 — Aca IP Gin Aph | ae 


Damit kann das gestellte Problem als gelést angesehen werden. 


7. Hinweise auf andere Lagerungen von achsensymmetrisch belasteten dicken Kreisplatten, 


Die in dieser Mitteilung angewandte Methode 14Bt sich auch fiir fest gestiitzte oder fest ein- | 
gespannte dicke Kreisplatten unter rotationssymmetrischer Last verwenden, wenn man die | 


Randbedingungen entsprechend andert. Die Resultate werden gegeniiber der elastischen Lage- 
rung einfacher zu gewinnen sein. 


a) Die frei gestiitzte dicke Kreisplatte. Das Problem ist fiir spezielle Belastungs- 


falle von Timpe! und Ohlig? behandelt worden. Die von diesen Verfassern formulierten Rand- | 


bedingungen, insbesondere fiir den Plattenmantel, fiihren zwangslaufig zu dem Resultat, daB 
die Auflagekraft durch die am Plattenumfang auftretenden Schubspannungen aufgenommen wird. 


Da eine solche freie Lagerung physikalisch schwer realisierbar ist, soll hier eine andere — durch | 
das de Saint-Venantsche Prinzip gerechtfertigte — Formulierung der Stittzung erfolgen. Die | 


gesamte achsensymmetrisch verteilte Last 


a 


P=2n{ f(r)rdr = (@&— a?) 7 Po, (65) 


soll in der Plattenebene z— h auf eine schmale, 


ringférmige ,,Stiitzflache“*, die in r= a starr ist 
gleichmaBig iibertragen werden (Abb. 3). Die 
Randbedingungen lauten 
t (25-0) == 0, (66) 
o.(0, 1) = — f(r). (67) 
TOs ria Oe : (68) 
(l) site 2 < @ ; 
BO! 1 69 
(h, r) (2 pa ee (69) 
Abb. 3. Die freigelagerte achsensymmetrisch belastete dicke t(h, r) =n) ° (70) 
Kreisplatte. 
G(2. 0) 0 (71) 


Mit Riicksicht auf die Forderung ¢(h, a) = 0 und die Randbedingung (66) erhalten wir fiir die 
Verschiebungen aus (11) und (12) 


Sy Ay Ay. % Eine : 
c= Di Pee + Oe — ae (ee M+ ey] Jy(hee) + Ae 
Ss = Ajit Ay h h Sos ah . (72) 
~ > [Pee + Q.e —— 213) (Pee k + 90%)) Jy(dua) — Ab 
= ST IP; isis Ay, (3 —4») q, ie 
: Pa = 20a Qe 4 xI—ayh| ° 


z —Ayz Ay | We) 
SESH) (Pre — Ke ) Jy (Anr) « 


Die Randbedingungen (67) und (68) fithren hier genau so wie bei der elastisch gelagerten Platte 
A. Timpe, Z. angew. Math. Mech. 4 (1924), S.361. 
2 R. Ohlig, Ing.-Arch. 13 (1942), S.155. 


1 
2 


a I 
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zu den Gleichungen (13) bis (15), bzw. (16), mit denen die Verschiebungen in 


Sd renee, [(3 Ay 2 A, 2) e tk + ok | as 


q fe aoc z : 
4A, (a5) le ef (3 — 4 — 2A, 2) e* ae 7. Sin dae Jo(Anr) 4- Az— 


—>| ans [((3— 4+ 2A,h) eh 1 eh » (72a) 
ry 4440. —2) zh)e +e aie 
q ae Fu 
au “Lege *Y + (8 — 40 — 2A,h) ol) 4 7- Sin dhl Jy(&ua) — Ah, 
=> Thana | Av — 2 A,2) yt ca ofk*| aie | 
a tia) 
| I [(3 4 24,2) elk? — Ape c;, l | 
4A, (1— 2) PS pepe” +e |— os Coj Mi J, (Ax) 
itbergehen, wahrend aus (69) zunichst [s. Gleichung (19)] 
1— 
26 1 Ae > lots 5 (A he “**_ Gin A, h) — Saw (ehe'* — Sinan 
0 e . (74) 
+ c, Co} Ak h we uri —— | eagh Whe IF § <i a, 
| l= Oia 7(a®—a3) 29 ay —< Tr <= a 


folgt. Multiplizieren wir diese Gleichung mit J,(/,7r)r dr und integrieren von r= 0 bis r= a, 
so bekommen wir wegen der Senne. der Eigenfunktionen 


Bee eet Arh) Fay) Mu he’*” — Sin Ah) + 0, Co Agh = | 


ian) | 
ae Sins ae yg Pay Jy (Ay 4) ay | 
A,PGCI3(A,a) a? 1, (a? —ai)CI3(A,a) 


(75) 


Der Ubergang zur schneidefirmigen Lagerung, d.h. a,—> a, lift sich mit Hilfe der L’ Hospital- 
schen Regel fiir die Koeffizienten vollziehen und man erhalt fiir diesen Fall 
AQ) P 
Gira 2a*2GJo(A;, 4) 


(76) 


Daf der Grenziibergang a,— a fiir z<h mit dem der Partialsummen vertauschbar ist, folgt aus 
der gleichmafigen Konvergenz der auftretenden Reihen beziiglich a,. Dieses bedeutet, daB die 
fiir a, --a gefundenen Lésungen fiir a, a wenigstens im Inneren der Platte (z<h) gleichmaBig 
gegen Grenzwerte streben. Man ist versucht, diese als Lésung des Problems fiir schneideférmige 
Lagerung anzusprechen. Freilich treten fiir z=h divergente Reihen auf, so daB der AnschluB 
an diesen Plattenrand (zh) so hergestellt wird, daB man die durch (74) gegebene untere Be- 
lastung vor dem Grenziibergang fiir sich betrachtet und in diesem Bestandteil a,>a 
vollzieht. 


Aus der Randbedingung (70) ergibt sich 


pr(Aphe *" + Sindh) + gy (Aphe“** + Sind, h) = 2(1— 2) Sindh. (77) 
Aus (75) und (77) lassen sich p, und q, berechnen: 
1—2y I sie : (1) Ay h ; | 
= | Aph + + Gin 2A,h Gin* fal) — Aphe" + Sind,h)) , 
Pk Gin? A, h— Ah? oe ( xh + q CM aA Sr Cle AL cy (Aphe® == a) | a 


— a ; 1 ‘ ; 
es oat ae |° ek (1) (A, he dich Sind, h) — cy, (ish + oi Sin 2h — Sint? Ayh)), | 
qe 
womit das Problem bis auf die Erfiillung der Forderung (71) gelést ist. Dieser noch nicht er- 
fiillten Bedingung kann man auch hier — genau so wie bei der elastischen Lagerung — im Sinne 
des Saint-Venantschen Prinzips gerecht werden, worauf jedoch an dieser Stelle nicht naher 
eingegangen werden soll. 
24 
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b) Die am Rande eingespannte dicke Kreisplatte. Auch in diesem Falle kénnen 
die Verschiebungen (72) und (73) herangezogen werden und neben: (66) erscheinen die Rand- 
bedingungen (67) bis (70) — letztere insbesondere fiir a, a — sinnvoll, wahrend (71) hinfallig 
wird. Da nach (73) die Radialverschiebung @ fiir r= a verschwindet, erleiden die Punkte des 
Plattenmantels, mit Ausnahme von z= h, nur vertikale Verschiebungen und dies entspricht ge- 
rade dem Wesen der Einspannung. [Eine andere Bestatigung fiir diesen Sachverhalt folgt aus 
(72): Es ist (Fr) =0]. Damit ist auch dieses Problem grundsatzlich erledigt, da nach Auf- 
lésung von (75) und (77) nach p, und g, vermége (13) und (16) alle Koeffizientenin den Ver- 
schiebungsfunktionen bekannt sind. 


8. Zusammenfassung und eine SchluBbemerkung. In einer vorangehenden Mitteilung* wurde 
das Problem der achsensymmet isch belasteten dicken Kreisplatte auf elastischer Unterlage 
untersucht und unter Vernachlassigung der zwischen Platte und Unterlage auftretenden 
Reibungskrafte die Lésung angegeben. Die vorliegende Mitteilung greift noch einmal dasselbe 
Problem auf unter Beriicksichtigung der erwahnten Krafte. Auch in diesem Falle fihrt das 
Randwertproblem auf ein — praziser gesagt zwei — Systeme von linearen Gleichungen mit 
unendlich vielen Unbekannten. Das unendliche Gleichungssystem ist das Charakterstische fiir 
den Spannungs- und Deformationszustand einer auf dem Halbraum aufliegenden dicken Kreis- 
platte, wahrend z. B. bei fester Einspannung die Koeffizienten in den Verschiebungsfunktionen 
sich direkt berechnen lassen, worauf hier abschlieBend hingewiesen wurde. 

Zum AbschluB noch eine Bemerkung: Den Randbedingungen (7), (8) und (9) liegen die Vor- 
aussetzungen zugrunde, daB zwischen Platte und Unterlage eine stetige Beriihrung besteht, 
was offenbar nur solange zutrifft, bis zwischen ihnen Druckkrafte auftreten; dies wird bei den 
iiblichen Belastungsarten meistens der Fall sein. Bei gewissen — in der Praxis allerdings kaum 
auftretenden — Belastungen kann es vorkommen, da in einem bestimmten, meist schmalen, 
am Rande gelegenen oder achsennahen Bereich neben den Druckkraften auch Zugspannungen 
auftreten, und das bedeutet in Wirklichkeit, daB die Platte sich von der Unterlage abhebt. Fiir 
die Praxis heift ein solches Resultat, die Belastung so zu andern, daB kein Zugbereich auftritt. 
Aus diesem Grunde hat es wenig praktisches Interesse, die einem solchen Fall entsprechenden 
wirklichen Spannungs- und Deformationsgleichungen zu ermitteln, was mit Hilfe einer ,,Ite- 
ration” prinzipiell méglich ware. 


(Eingegangen am 27.Marz 1951.) 
Anschrift des Verfassers: Prof. Dr.-Ing. Istvén Szabé, Berlin-Charlottenburg 2, Technische Universitat. 


Druckfehlerberichtigung zur ersten Mitteilung, Ingenieur-Archiyv 19 (1951) S. 128: 
Seite 129, Formel (11): = statt = 


Seite 130, Formel (21): fiige rechts hinzu: +A. 
Seite 131, Formel (31): in der zweiten Zeile e** statt e”*. 


fo) 
Seite 132, Zeile 6 von unten: >) statt ie 
k=1 k=v 


Seite 134, Zeile 1 von oben: o;(h, r) statt o(h,r). 

Seite 134, Zeile 11 von oben, zweite Formel rechte Seite: a@ statt @. 
Seite 135, Formel (54): 3—2y» statt 3—4y». 

Seite 136, Formel (57): mp statt mp. 


Formel (59), (60) und (61): statt 


o,h D Ma 
2G(1+ 7) 2G(1—y) ° 
Seite 137, Formel (63): (l—v) Jo(Ana). statt (1—y Jy) (Ana). 
Seite 138, Zeile 8 von obzn: > statt .») . 
a, B=1 a, B=1 
Seite 139, Zeile 3 von unten: fiige den Faktor open ersten Glied im Zahler hinzu. 
Seite 142, Ziffer 12, zweite Zeile: Hookeschen statt Hockeschen. 


1 Siehe FuBnote auf S. 342. 
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Neue Konstruktion der Wirkungslinie des resultierenden 
Massenwiderstandes eines eben bewegten Getriebegliedes. 


Von O. Tolle. 


Ausgehend von der Méglichkeit, bei Ermittlung einer bestimmten Wirkungslinie die gesamte 
Massenwirkung eines eben bewegten Getriebegliedes durch eine Einzelkraft P——mbs zu er- 
setzen, wobei m die Masse und bs der Beschleunigungsvektor des Schwerpunktes S ist, wurden 
von H. Alt’, K. Federhofer® und dem Verfasser® zeichnerische Verfahren angegeben, mit denen 
diese Wirkungslinie der resultierenden Massenkraft P ermittelt werden kann. Die nachstehende 
Konstruktion ist eine Anwendung des Beschleunigungsplanes von M.Tolle* unter gleich- 
zeitiger Benutzung des Tragheitsmittelpunktes® (Schwingungsmittelpunkt). Die resul- 
tierende Massenkraft P—=—  mbs (Abb. 1), Ersatz fiir — ¢ Js und — mbs im Schwerpunkt S angrei- 
fend, ist parallel bs, entgegengesetzt gerichtet 
und liegt von S aus im Abstand 


2 
—eJs = ts bisk tk bisx 
—mbs Sk bs —_ bs : 


e= 


Hierin bedeutet: b,sx die relative Tangential- 


beschleunigung S gegen K, ¢ die Winkelbeschleu- 
nigung, Js das Massentraégheitsmoment bezogen 


auf Schwerachse 1 KS, is = | Js/m den Triag- 
heitshalbmesser und t,= ig/sx den Abstand des 

auf K bezogenen Tragheitsmittelpunktes Tx SBBS A: puclrteeee Deetisiateng, der UY ehasgsiins der 
von S aus. 

Die Wirkungslinie der resultierenden Massenkraft P kann andererseits auch durch 
einen auf der Verbindungslinie KS liegenden Punkt A festgelegt werden, dessen Entfernung von 
S aus gemessen gemaBs Abb. 1 betragt 

SA 3 e tk bisk 
~~ > sin & ae bg sind) 


Wird fiir b,sx/sinx = x gesetzt, so ist a=txx/bs. Diese Entfernung a kann nach Abb. 2 in einfacher 
Weise mit Hilfe des Beschleunigungsplanes nach M. Tolle* folgendermafen gefunden werden. 


Gegeben sind die Be schleunigungen bx eines Punktes K und bysK (7) =ty 
die Schwerpunkthbeschleunigung bs, damit auch die Kompo- 


nenten der relativen Beschleunigung bsx von S gegen K, 


(k) 


und zwar die relative Normalbeschleunigung b,sx von S 
nach K gerichtet und die relative Tangentialbeschleunigung 
bisx L KS. Die Anfangspunkte (S) und (K) der Beschleuni- 
gungen bs und bx liegen im Abstand 6,sx parallel KS, die 
Endpunkte[S] und [K] auf einer Senkrechten zu KS; Strecke 
[K][S] =bsx. Man trigt im Beschleunigungsplan von (S) 
aus in Richtung (K) die Strecke t,x=i§/sx auf, Endpunkt 
(Tx), und verbindet (Tx) mit [S],zieht durch [K] die Parallele 
zu (K)(S), dann liegt auf dieser die gesuchte Strecke a 
zwischen bs und der Verbindungslinie (Tx)[S]. Aus Abb.2 
folgt namlich »= b,sx/sina und aus der Ahnlichkeit der Drei- 15] 
ecke a/x = tx/bs d.h. a= txx/bs. Ohne Aufzeichnen des Be- 


Abb. 2. Beschleunigungsplan nach M. Tolle mit 


schleunigungsplanes nach Abb. 2 ergibt sich hieraus in Anwendung zur graphischen Lagenbestimmung 


der resultierenden Massenkraft. 


1 H. Alt, Z. angew. Math. Mech. 6 (1926), 5.58. 
2 K. Federhofer, ZVDI 74 (1930), 5.234. 
3 OQ. Tolle, Ing.-Arch. 1 (1930), $.377, ZVDI 74 (1930), $.1453, AWF-Getriebeblatt 637/638. 


4 M. Tolle, ZVDI76 (1932), S.799. 
24* 
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Abb.3 folgende Konstruktion fiir die Lage der resultierenden Massenkraft P im 
Lageplan selbst. 

Sind wiederum bx und bs gegeben und mit ihren Anfangspunkten K und S eingetragen, so. 
wird von S aus die Strecke tx i3/sx = S(Tx) nach K hin aufgetragen (und nicht in Richtung |, 
zum Trigheitsmittelpunkt Tx). Der Endpunkt (Tx) wird mit Endpunkt [S] von bs verbunden. | 
Die Parallele zu KS durch den Endpunkt [K] von bx schneidet zwischen der Verbindungslinie | 
(Tx) S und ‘bs die gesuchte Strecke a ab. In Abb.3 ist dann die so gefundene Strecke a—=SA 


von Saus nach aufen aufzutragen. 


| 1s / 
: / Abb. 4. Konstrukt’ 
= wie in Abb. 3 fir ; 
[S] bk Hi dere Werte der 
Abb. 3. Zeichnerische Ermittlung der Wirkungslinie 4 4 Me Sea eo i schleunigungen 
des resultierenden Massenwiderstandes im Lageplan a bo und to 
(ohne Aufzeichnen des Beschleunigungsplanes). a K Ss 


Dieselbe Konstruktion in Abb.4 durchgefiihrt, in welcher andere Werte fiir bx und bs an- | 
genommen wurden, ergibt einen Wert a, der von S aus nach K hin aufzutragen ist. Es gilt 
folgende Vorzeichenregel fiir a: 

1. Liegt die gefundene Strecke a von [S] aus auf derselben Seite wie die Strecke S(7Tx), dann 
ist a= SA von S aus auf der Verlangerung von KS aufzutragen (Abb. 3). 

2. Liegen die gefundene Strecke a und Strecke S(7'x) auf verschiedenen Seiten von [S], se 
ist a= SA von S aus nach K hin aufzutragen (Abb. 4). 


Abb. 5 und Abb. 6 
zeigen eine weitere Ver- 
einfachung der Kon- 
struktion: Die Verbin- 
dungslinie des auf K be- 
zogenen Tragheitsmit- 
telpunktes Tx mit dem 
Endpunkt/S]derSchwer- 


punktsbeschleunigung bs | 
(bzw. deren Verlange- 
f rung) schneidet sich mit | 
Abb. 5. yy ‘ era Abb. 6. der Parallelen zu KS, | 
Ky "A'S durch den Endpunkt[K] | 
Abb. 5 und Abb. 6. Ejnfache Konstruktion der Wirkungslinie des resultierenden Massenwider- : aR 
standes, wobei sowohl die Ubertragung der Strecke a als auch die Vorzeichenregel entbehriich ist. der Beschleunigung bx 
gezogen, im Punkt 4’, | 
durch welchen die resultierende Massenkraft P— —mbs hindurchgeht (Strecke S’A’ = 
a= SA). Liegen die Verhaltnisse wie in Abb. 4, so zeigt Abb.6, da®B auch hier die letzt- 
genannte Konstruktion durchfihrbar ist. Es braucht hierbei also weder die Strecke a iber- 
tragen zu werden, noch ist eine besondere Vorzeichenregel zu beachten. 
Die vorstehend angegebene Konstruktion ist nicht nur von allen bisher bekannten die ein- 
fachste, sondern sie ist dann, wenn sonst schleifende Schnitte vorkommen, die genauere. 


Poke 


(Eingegangen am 27, Marz 1951.) 
Baurat Dipl.-Ing. 0. Tolle, Konstanz (Bodensee), Bahnhofplatz 10. 
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Massenreduktion und Massendruckdrehkraft fiir die Schubstange 
eines Kurbelgetriebes. 


Von O. Tolle. 


1, Einleitung. Das zur Schwungradberechnung allgemein benutzte Massenwuchtdia- 
gramm?ist eine Zusammenstellung der Kurve der reduzierten Masse und der Arbeits- oder 
Wuchtkurve. Neben dieser Anwendung kann aus der reduzierten Masse auch die Massendruck- 
drehkraftkurve entwickelt werden. Soweit es sich um den Einflu® der Schubstange auf 
Massenreduktion und Drehkraft handelt, werden nachfolgend Gleichungen aufgestellt, welche 
den bisher tiblichen Ersatz der Stangenmasse durch je eine im Kreuzkopfbolzen und eine im 
Kurbelzapfen angreifend zu denkende Masse noch verbessert, und aus denen dic von der Stange 
herrithrende Massendruckdrehkraft fiir jede Kurbelstellung mit gréBerer Genauigkeit als bisher 
bestimmt werden kann. Es wird ferner eine einfache Konstruktion angegeben, mit welcher 
sich die Drehkraft bestimmen l48t, wenn nach einem der bekannten zeichnerischen Verfahren2 
GréBe und Wirkungslinie der resultierenden Massenkraft der Schubstange gefunden worden ist. 
AbschlieBend wird noch eine neue graphische Lésung angefiihrt, wie mit Hilfe des Massenwucht- 
diagramms die Kurbelzapfengeschwindigkeit fiir die verschiedenen Getriebestellungen gefunden 
werden kann. 


2, Reduzierte Masse der Schubstange. Unter der reduzierten Masse M, eines Getriebegliedes 
bezogen auf einen beliebigen Punkt versteht man die Masse gleicher Wucht. Als Reduktions- 
as wahlt man im Kurbelgetriebe fast immer den Kurbelzapfen K. Im Folgenden bedeuten 
s. Abb. 1): 

r den Kurbelhalbmesser, / die Schubstangen- 
lange, A= r/l das Stangenverhaltnis, x den 
Kurbelwinkel, 6 den Erhebungswinkel der 
Stange, w= dax/dt= a’ die Winkelgeschwindig- 
keit der Kurbel, m,—df/dt—/p’ die augen- 
blickliche Winkelgeschwindigkeit der Stange 
um den Kreuzkopfbolzen B, v= rw die 
Geschwindigkeit des Kurbelzapfens K, 


aa) [sin o+ 54 sin 2 a) die Kreuzkopfge- 


a 
schwindigkeit, s, den Abstand des Schwer- baa 


Ki 
punktes S,der Stange vom Kreuzkopfbolzen B, ie ae: 
a 1— So den Abstand des Schwerpunktes S5 Abb. 1. Geschwindigkeitsverhiltnisse im Kurbelgetriebe. 
der® Stange vom Kurbelzapfen K, M, die 
Schubstangenmasse, J, das Tragheitsmoment der Stange bezogen auf Kreuzkopfbolzen- 
mittellinie, M,, die auf Kurbelzapfen K reduzierte Masse der Schubstange. 


Fir die Schubstange eines Kurbelgetriebes gilt® 


WY Mim (M- 1" sin? (1+24 cosa) +72 (1) 
oder 
Nee Ape econ? hie wewivenee an 
a =| »— FIL —cos2 4+ A (cosa— cos ONT ae 


Hiernach setzt sich die gesamte auf den Kurbelzapfen reduzierte Masse M,, der Schubstange 
aus zwei Teilen zusammen, 


1 F, Wittenbauer, Z. VDI (1905), S.471, sowie Graphische Dynamik, 5.759. Berlin 1923; M. Tolle, 
Regelung der Kraftmaschinen S.104, 3. Aufl. Berlin 1921, C. B. Biezeno und R, Grammel, Technische Dy- 
namik §.935. Berlin 1939. ; 

2 H, Alt, Z. angew. Math. Mech. 6 (1926), S.58; K. Federhofer, Z. VDI 74 (1930), 5.234; O. Tolle, Ing.- 
Arch. 1 (1930), $.377; Ebenda XIX (1951), 8.355; Z.VDI, 74 (1930), $.1453, AWF-Getriebeblatt 637/638. 

3 M. Tolle, Regelung der Kraftmaschinen 5.107/119 Gl. (34), Gl. (84a) u. Gl. (41). 
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1. aus dem konstanten Wert M;,—= = im Kurbelzapfen K und 


2. aus dem veradnderlichen Wert 


i (m— 2) sin?« (1-+2/ cosa) = ; (Ms 22) [1—cos2 aA (cosa— cos3 a) | 5 


im Kreuzko pfbolzen B zu denken. Vorstehende Gleichungen sind unter Fortlassung der Glieder 
mit J? und mit Hilfe von Ersatzmassen aufgestellt worden. Es folgt nun eine Ableitung fiir die 
reduzierte Masse einer Schubstange ohne Massenersatz. Nach Abb.1 ist die Wucht E, 
der Schubstange! 
VEN CAME) Os : 
Eo 5 z 2 _ Myc Wes, sin B 
Mae? | J2(6')* Pele 
=< 5 ano — M,c B’s, sin B . 
Mit Kurbelzapfenmittelpunkt K als Reduktionspunkt, dessen Geschwindigkeit v ist, wird 
Mr. v 
c= = 
und somit 
2E, 
NM Mes a . 

: A fe! v Hud, : 

Es‘ist o/ = ee, Cob fi== = =,, ferner sinB=A/sina, somit 
dsinB _—, dsina ___, dsina da dsinB dB ; 
ee Us pe re ap de os On as ae cos B 
und damit 
,_.4008% = cosa sv COS 
dia leosB lL cosp 
a) Erste Annaherung: Fir kleine Winkel 6 kann cosf =1 gesetzt werden; dann wird 


i + cosa 


und 


2 all 
Ma=M,(<) Jolt cosa | = —2M,¢ (+ cosa) ie sind 
/ 


ONE bh 
=m, (=) + 53 cosa — 2M, — =P Asina Cosa. 


b) Zweite Annaherung J?=0. Hier ist c=v(sina ++ sin 2a) =vsina (1+ Acosa), 


Cc 


— =sina(l-+Acosa) und oles sin? « (1+2A cosa); 


v v 


folglich wird 


M,.= My sin?« + 22M, sin? « cosa + as cos*a — 2A M, = sin? o cosa. 
Mit 


1 l 


kommt 


M,. = Mysin?« + 22M, + sin? « cosa +- = cos? a. 


1 Dn 
Nachweis: Es erfolgt Drehung der Schubstange um Momentanpol $$ mit , = = = f’, folglich ist die 


Jz BYP 
W a ee ; 
ibe Es Ge op RS J% = Jsz + M, e? = J, + M,(e?—s3). 
ach dem Cosinussatz ist e?—s? = 0? 29s,sinf und somit Jy = J, -+ M,(0?—2¢s, sin), und mit der 
Kreuzkopfgeschwindigkeit ¢ = ef’ =o, wird E,= Jn} a Macs 


9 "ou —M,co,s,sinf. 
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Erganzt man die rechte Seite durch 
H bre ade ; 
Ta sin? « — Ja sin? . 
so wird 


M,,= (M.— #2) sin? a -+ 2A 4 M, sin? « cosa + a 


oder 


M,2= LY pee sin?a /1-+ 24“) - Ms cosa\ + 22, a 
P I ip P 


Hier ist a) = M,, der im Kurbelzapfen zu denkende Teil der reduzierten Schubstangenmasse, 


wahrend der erste Summand des Ausdruckes fiir M,, der auf den Kreuzkopf entfallende Anteil 
M2 ist. Dieser kann nach (2) rechnerisch oder nach der in M. Tolle ,,Regelung der Kraft- 
maschinen“ in Abb.93 Seite 108 fiir hin- und hergehende Teile angegebenen Konstruktion auch 
fiir die reduzierte Masse der Schubstange, wie folgt, zeichnerisch ermittelt werden (Abb. 2). 


Man zeichne mit p= 2/ ; M, als Halb- 


messer einen Kreis, von dessen Mittel- 
punkt M die Kurbelstellungen ausgehen 


und errichte auf 0— 0 im Abstand M, — ah 


12 
von M aus die Senkrechte. Abstand 
AC= Z, von M aus auf der Kurbelstellung « 
aufgetragen, liefert die Strecke DE=Z’, 
welche von M aus aufgetragen mit Ab- ) 
stand FG = Z” die gesuchte reduzierte 
Masse M/, ergibt, soweit diese im Kreuz- 
kopfzapfen zu denken ist. Denn nach 
Abb. 2 ist 


db 
M,- se 


Abb, 2. Zeichnerische Ermittlung des auf den Kreuzkopf 
entfallenden Anteils der reduzierten Schubstangenmasse. 


Z =(Ma— 7) +225} Macosa., 


Li — Jin 
2 == 7’ sing — Z sin* 6 = (M.— = sin? a 2A ; a ae = Me. 
if pee ees 
cage 
Fiir die prismatische Stange mit ; = 2 und 2 = M, gilt nach (2) 
fir AA.) = M, sin? a(1 + = Aeoss' —- a ; 
fir J—0: Mee M, sin? a + "2 ; 
Mit einem Stangenverhaltnisse 1 = a werden die fiir die graphische Lisung erforderlichen 
Werte 
; 2 ie Sk 
a) bei prismatischer Stange M, — 2 = Vise en gs M,=0,2 M,, 
b) bei einer Schubstange mit =} = 0,35 und = 0,60 M, 
dagegen 
,—2=04M, und 9=20,35M,=0,14M,. 


Abb.3 zeigt die Kurven fiir die reduzierte Schubstangenmasse M,, und zwar Kurve a fiir eine 


; J. : 
prismatische Stange, Kurve b fiir eine Schubstange mit 3 = 0,35 und M,— 7 — 0,4M,. Der im 


Ingenieur-Archiv 
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Kurbelzapfen zu denkende Teil der reduzierten Masse ist im ersten Fall i — 3 M,, im “wei- 
ten Fall M7, = 0,6 M,. 


M; 


C2 


0 4 x 3a 2 


Abb. 3. Kurven der reduzierten Massen fiir 2 = 1/5. 
1. prismatische Stange Kurve a; 2. Schubstange mit J,/l? = 0,6 M, und s,/l = 0,35 Kurve b. 


Es 148t sich nun nachweisen, da die durch (2) ausgedriickte Kurve der reduzierten Masse 
fiir die Schubstange sich aus drei Sinuskurven mit den Frequenzen w, 2m und 3m zusammen- 
setzt, wie dies auch bei der Kurve der reduzierten Masse nach (la) der Fallist. Setzt man 


1 ay nore afi 
sin? « = (1 — cos 2) und sin* a cosa = -> sin 2a sina = (cosa — cos3a) 


so ergibt sich aus (2 


BE 


2 


Nimmt man den letzten Ausdruck der Gleichung (3) in die Klammer hinein, so kann fiir die 
reduzierte Masse der Schubstange auch geschrieben werden 


M,.= : (Ms dls ) (Ms B) cos2a+A = M, (cosa — cos 30)| 5 (4) 


aay (a ) (Ms fh) cos 20+ At Mz (cosa — cos 3)|-+ 2 : (3) 


2 P 
3. Massendruckdrehkraft T, yon der Schubstangenmasse herriihrend, Die Wucht der Schub- 


stange ist 
| L ge i 


E,= + M,,° => M,, (ro)? 


am 
AX und nach dem Wuchtsatz gilt (Abb. 4) 


7, M 
dE dy 
Sey 


- “ 
Abb. 4. Massendruckdrehkraft dt i dt 
infolge Schubstangenmasse. 7 0 

. folglich ist 


(es 1 dE, __ 1 dE,d« _ 1 dE, Sine dM,, _ro* dM,, 


rao dt ro da dt r da ree dx 2 dx 
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Differentiert man M,, in (3) nach a, so erhalt man 


fa eae 2 (M. — ) sin 24 —A + M, (sina — 3 sin3q) 


und damit 
: (6) 


Diese Gleichung fiir T, enthalt nur reine Sinusglieder mit den Frequenzen w, 2@ und 3@. Man 
bestimmt die Massendruckdrehkraftkurve daher am einfachsten graphisch wie bei der Zusammen- 
setzung harmonischer Schwingungen mit verschiedenen Frequenzen. Abb.5 zeigt fiir eine pris- 


rw? ), 
= Te (ma #) sin 20% — + “4 M, (sina — 3 Sin 34) 


A 
0 x a 3a 2a 


Abb. 5. Massendruckdrehkraftkurven fiir 1. prismatische Stange Kurve a; 
2. Schubstange mit J,/l? = 0,6 M, und s,/l = 0,35 Kurve b. 


: : 1 1 eas 
matische Stange aie My, und p= 7 beid= = die Massendruckdrehkraftkurve a, 


ie 3 l 
die durch algebraische Addition der Kurven c, d und e gefunden wurde. Fiir die prismatische 
Stange gilt nach (6) 


9 : 
T= "T |g Me sin 2 « 3 Ma(sina—3 sin 3a) 
rw” ; : 9 
== M, 0 Et y 
- (sin2a 0 sINA-T 4H sin3.) 


In Abb.5 ist sin2 «durch Kurve c, = sing durch Kurve d und a sin3a durch Kurve e dargestellt. 


Die Kurve 6 ist die Massendruckdrehkraftkurve fiir eine Schubstange mit a = 0,6 M, und 


7 == (Ousa hel A == = und in gleicher Weise gefunden wie fiir die prismatische Stange aus 


2 
T, = M, (0,6 sin2.x— 0,0525 sina +-0,1675 sin 3a). 


Ist (Abb.6) die resultierende Massenkraft P—-— M,b, der Schubstange (wo b, die Schwer- 
punktsbeschleunigung der Schubstange) nach GréSe und Wirkungslinie nach einem der zeichne- 
rischen Verfahren (FuBnote 2 auf S$. 357) ermittelt worden, so kann ihr Einfluf auf die Dreh- 
kraft T, in einfacher Weise wie folgt bestimmt werden. Man zerlegt Pin zwei Komponenten, 
von denen die eine P’ im Kreuzkopfzapfen senkrecht zur Gleitbahn angreift, mithin keine 
Arbeit leistet, also auch keine Drehkraft liefert, wahrend die andere P’ im Kurbelzapfen an- 
greift (also durch den Schnittpunkt 0 von P und P” hindurchgeht). Am Kurbelzapfen sind 
dann im Gleichgewicht: P’, T, und die in die Richtung der Kurbel fallende Kraft R (s. Krafte- 
plan Abb. 6a). 
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Ist, wie in Abb.7, der Schnittpunkt der Normalen zur Gleitbahn mit der Wirkungslinie von P | 
nicht zuganglich, so findet man die Richtung von P’ nach dem Strahlensatz, indem man mit K 
als der einen Ecke eines beliebigen Dreiecks beginnend, die anderen Ecken auf die Wirkungs- 
linien von P und P” legt. Ein hierzu ahnliches Dreieck (durch Parallelenziehen) ergibt einen 
weiteren Punkt der Wirkungslinie P’ und damit auch deren GréBe nach Krafteplan Abb. 7a. 
Fiir den Kurbelzapfen findet man dann wie oben aus P’ die Drehkraft T, und R. 


4, Vergleich der neuaufgestellten _ 

Gleichungen fiir die reduzierte Masse 

und fiir die Massendruckdrehkraft | 

pu einer Schubstange mit den bisher | 
Abb6a %ebriuchlichen Formeln. 


a) Reduzierte Masse. Beim 
Vergleich des in (3) gefundenen Wer- 
tes fiir M,, mit dem nach (la) findet 
man, daB in beiden Gleichungen die 
konstanten Glieder und diejenigen 
mit cos 2a vdllig iibereinstimmen. 
Der iibliche Ersatz der Stangen- 


masse durch = direkt an der Kurbel 


aid ae und M,— & am Kreuzkopf ist also 
es vollstandig gerechtfertigt, soweit die 
gréBeren Glieder in Betracht kom- 
R men. Bei den Gliedern mit 2 ware 
y an Stelle von M,— S der Gleichung 
, 
wee s a 
Abb. 6 und Abb. 7. Zeichnerische Abb.6a und Abb. 7a, Die zu den (1a) richtiger M, 7 zusetzen. Priift 
i e r r, i er Abb. und Abb. 7 zugehéri is os 7 
gre tie et ees eerie far Kiaftecke, So Man den Unterschied z.B. fiir eine 
Stange mit at = 0,6 M, und ; 0,35 
: I 
bei A = —, so ware an Stelle von 
Ley V ip Solel 2 
ERM ENT SD yi 
SO TAURI ; M,— 0,6M,= 0,4 M, 
4G My) S935 ay genauer 0,35 M, zu setzen; dieser Unter- 
“Qo =r) los schied von 0,05 M, ist noch mit 4 zu mul- 
: tiplizieren, so da® ein Fehler von 
\. on Ah Ss ih 
VS Vi ; i (Mm — = M, ‘) = = 0,05 M,=0,01 My, 
\ i und gegen 0,4M, der Glieder ohne 4 ein 
b / solcher von 0,01/0,40 = 2,5°% besteht. Bei 
\ co prismatischer Stange mit af = = M, und 
0 x a 3a en 1 
2 2 Sie : —_ 1 Al fes 2, 
Abb. 8. Prozentuale Abweichungen der auf den Kreuzkopf ent- we pea tas A ma Si g cane UE ra me 2 3 me 
fallenden Anteile der reduzierten Schubstangenmasse nach (1) und Sy 1 : Lipa 
(2) in Abhingigkeit yom Kurbelwinkel. und 7 C M,; => 9 M, ist der Fehler = “Tx M, 
= 39 Mz und damit gegen 3 M, der Glieder ohne A dieser 36/5 = 5%, d.h. doppelt so groB 


wie oben. 


_Bestimmt man aus den Gleichungen (1) und (2) den auf den Kreuzkopfzapfen entfallenden 
Teil M,» in beiden Fallen, so ergeben sich fiir die Schubstange je nach der Kurbelstellung ver- 
schiedene Werte. So ist z.B. fiir eine prismatische Stange mit 

JP tt Sy 1 


1 0 
1B =3 ™M,, [= Sr? Ne und « = 30° 
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mach (1) M,,= 0,2243M,, nach (2) Mr3=0,2099M,. Das Verhiiltnis 0,2243 M,/0,2099 M, 
= 1,0695 besagt, daB im ersten Fall Mf, um 6,95% gréBer ist. 


Abb.8 zeigt die prozentualen Abweichungen der Werte M/5 in den verschiedenen Kurbel- 


stellungen fiir eine prismatische Stange und eine Schubstange mit ; = Op ud J3 0,6 M, 


: i 
bei A = =: In den Totlagen sind die Abweichungen am griBten, fiir x = > und “== =  liefern 
(1) und (2) dieselben Werte fiir M7}. 


Fir unendliche Stangenlinge (A—0) stimmt der gebrauchliche Ersatz, namlich bei pris- 


matischer Stange 3 M2. bei beliebiger Schubstangenform 2 im Kreuzkopfzapfen in 


allen Kurbelstellungen mit den Werten der neuen Gleichungen genau iiberein. 


b) Massendruckdrehkraft. Nach Gleichung (5) gilt allgemein 


ro? dM, 
et veer 


Differentiert man M,, in Gleichung (la) nach «, so wird 


dM,, 1 
de om 2 


2 
(— 7) [2sin2«— A (sina — 3sin3«a)] 


und somit 


Tw? 


i ¢ (Me i) [2sin2a—A (sina —3sin3a)}. 


Hieraus ergibt sich z. B. fiir die prismatische Stange (M— T= z M) 


T,="@" = M, [2sin20—A (sina —3sin3a)] 


rw? , : ; 
are 3[8sin2a—4A (sina —3sin3a)}. 


Der entsprechende Wert fiir T, nach Gleichung (6) lautet 


5 : 1 : ; 
T= + M,sin2a— > 5 M, (sina — 3 sin 3a) 


rw" 2 f : 
= oe M,[8sin2 « — 3/A(sina— 3sin3)]; 
der Unterschied liegt also auch hier nur bei dem Glied mit A. 


5. Konstruktion der Kurbelzapfengeschwindigkeit v mit Hilfe des Massenwuchtdiagramms. 
Eine einfache graphische Lisung ist bereits in ,,Regelung der Kraftmaschinen“ von M. Tolle, 
3. Auflage Seite 121 u. f. in den Abb. 101 bis 103 angegeben. Wie dort ausfiihrlich beschrieben, 
gilt fir das Massenwuchtdiagramm, in welchem die auf den Kurbelzapfen reduzierten Gewichte 
G,—= M,¢ als Abszissen, die zugehérige Wucht E— $ M, v2 als Ordinaten aufgetragen sind, 
folgendes: Zieht man in Abb. 9 von dem Koordinatenanfangspunkt O aus Strahlen nach der 
Massenwuchtkurve, so gilt fiir den Winkel — eines solchen Strahles 


E v r? w? ah r 7 am, 
beG te 2 27 Le 


d.h., die augenblickliche Kurbelzapfengeschwindigkeit v (bzw. m und n) ist dem Wert |/tg p 
proportional. Fiir ein rechtwinkliges Achsenkreuz (x, y) mit x= v und y=tgy= v*/2g = 0/28 
und Koordinatenanfangspunkt 0’ (Abb.9a) ergibt sich eine Parahel, deren Ordinaten die Ge- 
schwindigkeitshéhen v?/2 g und deren Abszissen die Kurbelzapfengeschwindigkeiten v sind. 
Bei Beriicksichtigung der Zeichenmafistabe fiir Wucht E mit 1 em A ¢emkg und reduziertes 
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Gewicht G, lem “4 y kg wird 


pee te pee ares 
v= 2s By | tgp baw. oe 2g —itgy—y, 


und mit v= x entsteht wiederum eine Parabel (Abb. 9). 

Sollen fiir die verschiedenen Getriebestellungen die Geschwindigkeiten v (w oder n) auf- 
gesucht werden. so zeichnet man diese Parabel am einfachsten und genauesten nach der be- 
kannten Hiillkonstruktion, indem man zunichst fiir einen beliebigen Winkel y die Kurbelzapfen- 
geschwindigkeit v = Y(e/y) 2 g- jtey 
ausrechnet. Legt man eine beliebige 
£ Ordinate y=O’A und tragt in 
einem beliebigen Mafstab 1 cm 
‘ §6m/s von A aus die errechnete 
Kurbelzapfengeschwindigkeit v als 
Abszisse AB auf, dann sind die fiir 
die Hillkonstruktion der Parabel 
erforderlichen Stiicke bekannt: der 
Scheitel O’ und die Scheiteltangente 
x, ferner fiir die gewahlte Ordinate 
y = tg y = (y/e) (v?/2 g) = OVA, die 
zugehorige Abszisse v= AB, d.h. 
QO Orin TAS ein zweiter Parabelpunkt B. Der 

a Schnittpunkt T der Parabeltangente 
Abb. 9 und 9a. Konstruktion der Kurbelzapfengeschwindigkeit mit Hilfe in B mit der Scheiteltangente x liegt 


des Massenwuchtdiagramms. Ast 
in 1/2 AB von O’ aus. 


Fiir jeden weiteren Punkt der Massenwuchtkurve, z. B. P,, findet man die zugehérige Ge- 
schwindigkeit v,, indem man durch P, den Strahl von O aus mit der Ordinate yin A, zum Schnitt 


bringt. A, B, ist dann die gesuchte Geschwindigkeit v,. In Abb. 9 sind auch die GréBen von Unax 
und Uni, eingezeichnet, die bekannterweise aus den Fahrstrahlen gefunden werden, welche oben 
und unten tangential an die Massenwuchtkurve gelegt sind. Werden an Stelle der Geschwindig- 
keiten v als Abszissen in beliebigem Mafstab die Winkelgeschwindigkeiten @, bzw. minut- 
lichen Drehzahlen n aufgetragen, so kénnen gleichfalls mit der Parabel fiir jede Getriebsstellung 
die zugehérigen Werte w, bzw. n, abgelesen werden. 

Die vorstehend angegebene Konstruktion kann auch fiir die um gekehrte Aufgabe be- 
nutzt werden, d. h. fiir gegebene Werte v (@ oder n) die zugehérigen Fahrstrahlwinkel ~ zu be- 
stimmen (Schwungradberechnung). 


Abb. 9 


(Eingegangen am 27. Marz 1951.) 
Anschrift des Verfassers: O. Tolle, Konstanz (Bodensee), Bahnhofplatz 10. 
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Die Berechnung der Profilform bei vorgegebener Geschwindigkeitsverteilung. 
Von E, Truckenbrodt, 


1, Einleitung. Das Problem, fiir eine vorgegebene Profilform die Geschwindigkeitsverteilung 
an der Profiloberflache zu berechnen, ist durch die von F’. Riegels! weiterentwickelte Singulari- 
tatenmethode auf sehr einfache Weise gelést. Nach diesem Verfahren laBt sich die Geschwindig- 
keitsverteilung dadurch berechnen, da8 man an bestimmten Stellen der Profilsehne die Profil- 
ordinaten mit fiir diese Stellen charakteristischen Konstanten multipliziert und dann iiber diese 
Produkte summiert. Auch die aus der Geschwindigkeitsverteilung sich ergebende Gesamtkraft 


und das Gesamtmoment lassen sich nach einer erganzenden Arbeit von E. Truckenbrodt? durch 
einfache Summenformeln bestimmen. 


Haufig interessiert aber auch das umgekehrte Problem, daB die Geschwindigkeitsverteilung 
vorgegeben ist, und die dazu gehérende Profilform gesucht wird. Diese Frage ist von Bedeutung 
beim Entwurf von Profilen, fiir die besondere Forderungen beziiglich Laminarhaltung, Verhal- 
tens bei hohen Mach-Zahlen oder Kavitationssicherheit gestellt sind. 

Die Behandlung dieses Problems ist nicht neu. Sie wurde erstmalig von A. Betz? auf der 
Grundlage der konformen Abbildung in Angriff genommen und spater von W. Mangler* weiter- 
entwickelt. Beiden Verfahren haften aber gewisse Schwierigkeiten an, die F’. Riegels®® veran- 
laBten, die Aufgabe mit Hilfe der bei der Berechnung der Geschwindigkeitsverteilung so erfolgreich 
angewandten Singularitatenmethode zu lésen. Nach diesem Verfahren ist es méglich, das Pro- 
blem allgemeiner als mit der konformen Abbildung zu behandeln und insbesondere die Schwierig- 
keit der Erfiillung der SchlieBungsbedingung fiir die Profilkontur durch einen Fourier-Reihen- 


ansatz fiir die Profilkontur zu beheben. 


Mittels des Quell-Senken-Verfahrens wird von Riegels der Fall des symmetrischen Profils ohne 
Anstellwinkel gelést. Dabei wird die Aufgabe zuriickgefiihrt auf die schrittweise Berechnung einer 
nichtlinearen Integralgleichung fiir die Profilneigung dy/dp, wobei y die Profilordinate und 
x= tl cos@ die Profilabszisse (Sehne) bedeuten. Die einzelnen Koeffizienten der Fourierschen 
Reihe fiir y erhalt man durch eine Analyse aus der vorgegebenen Geschwindigkeitsverteilung. 
Da fiir die Fourier-Analyse ein Analysator zur Verfiigung stand, ist Riegels bei der Fourier-Reihe 
stehen geblieben. Es liegt natiirlich nahe, ahnlich wie es Riegels schon bei der Berechnung der 
Geschwindigkeitsverteilung aus einer vorgegebenen Profilkontur getan hat, die Bestimmung der 
einzelnen Fourier-Koeffizienten durch Anwendung der Gaufschen Quadratur zu umgehen und ein 
Summenverfahren zu entwickeln, bei dem nur iiber die Produkte des vorgegebenen Geschwindig- 
keitsbetrages an gewissen festen Stellen der Sehne mit festen ein fiir allemal berechneten Kon- 
stanten zu summieren ist. Hieriitber wird im Abschnitt 2 berichtet. 

Auch der Fall des Profils mit Wélbung und Anstellwinkel wurde von Riegels behandelt. Ab- 
gesehen davon, das Riegels in seiner neueren Arbeit’ zur Berechnung der Geschwindigkeits- 
verteilung die Erfassung des Anstellwinkeleinflusses anders vornimmt als in seiner alteren Ar- 
beit® zur Berechnung der Profilkontur, 1aBt sich besonders dieser Abschnitt seiner Arbeit noch 
erheblich erweitern, was ja wohl auch schon von Riegels beabsichtigt war, wie aus einer Bemer- 
kung auf S.7 seiner alteren Arbeit hervorgeht. Auch bei der Behandlung des unsymmetrischen 
Profils bleibt Riegels bei der Fourier-Reihe stehen. Es ist daher Gegenstand der Abschnitte 3 
und 4 auBer der oben erwahnten Erweiterung und der Erfassung des neuen Anstellwinkelein- 
flusses auch die Ableitung des entsprechenden Summenverfahrens. 

Bei der Ableitung unserer Formeln wollen wir uns nicht auf die alteren Arbeiten von 
F. Riegels beziehen, sondern sie aus der neueren Arbeit herleiten. Dieses geschieht einmal, weil 
die alteren Arbeiten schwer zuganglich sind, und zum anderen mit Riicksicht auf die neue Er- 


1 F., Riegels, Ing.-Arch. 16 (1948), $.373 (I. Mitteilung), 17 (1949), 5.94 (II. Mitteilung) und 18 (1950), 
S.329 (Berichtigung). 

2 KE. Truckenbrodt, Ing.-Arch. 18 (1950), 5.324. 

3 4. Betz, Luftfahrtforschung 11 (1934), S.1. 

4 W. Mangler, Jahrbuch 1938 der deutschen Luftfahrtforschung, S.1I 46. 

5 F. Riegels, Z. angew. Math. Mech. 24 (1944), 8.273. 

6 F,. Riegels, Untersuchungen und Mitteilungen der deutschen Luftfahrtforschung Nr.3019 (1943). 
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fassung des Anstellwinkeleinflusses. Wir werden also stets auf die Formeln der I. bzw. II. Mit- 
teilung der neueren Arbeit zuriickgreifen. 


2. Symmetrisches Profil ohne Anstellwinkel, Ist die Geschwindigkeitsverteilung auf Ober- 
und Unterseite des Profils an gleichen Stellen der Profilsehne die gleiche, so entspricht dieses 
einem symmetrischen Profil! (Profiltropfen) y(), dessen Sehne in Strémungsrichtung verlauft, 


d.h. das nicht angestellt ist. 
Formelmafig kénnen wir fiir diese Geschwindigkeitsverteilung schreiben 


W=V-+w,. (1) 
Hierin bedeutet V die Anstrémgeschwindigkeit und w, die durch die Verdrangung des Profils 
entstehende Zusatzgeschwindigkeit. An der Vorderkante des Profils bildet sich ein Staupunkt, 
d.h. es muB dort sein W=0. 

Die Berechnung von w, geschieht mittels einer auf der Profilsehne angenommenen Quell- 
Senkenverteilung. Aus der Bedingurg, daB an der Profiloberflache keine Normalgeschwindig- 
keiten herrschen diirfen, findet man nach F. Riegels, Gleichung (9) der I. Mitteilung, die folgende 
Integralgleichung fiir die Profiltangenten dy/dx der Profiloberseite : 


12 
LGU =: ae : [a@e) | (2) 


dy\2 qt x—% 
ya ae (2) —I2 
Fihren wir den Zusammenhang 
ee! 
x= on cosp (3) 
ein, dann wollen wir den gesuchten Profiltropfen durch die Reihe 
l 5 
ee (4) 


beschreiben. Dieser Ansatz bietet die Gewahr dafir, da® sich das Profil an der Vorderkante 
(p=) und an der Hinterkante (p—0) schlieBt. 

Die Lésung der Gleichung(2) kénnen wir unmittelbar von F. Riegels, Gleichung (12) der 
I. Mitteilung, itbernehmen: 


fsing = 3) vb, sinyp. (5) 
I 
Die neue Abkiirzung f wollen wir als Geschwindigkeitsfunktion bezeichnen; sie berechnet sich zu 
W 
pee yal, (6) 


wobei 
x dy\? 
x= 1+ (2) (7) 
einen Faktor darstellt, der nur an den Stellen groBer Profilkriimmung (Vorderkante) von Be- 
deutung wird. 


Ist die linke Seite der Gleichung (5) bekannt, so lassen sich die Koeffizienten » b, durch eine 
Fourier-Analyse bestimmen. Sie liefert 


yb, = - [ £@ sing’ sinyq’ dy’. (8) 
0 


Differentiieren wir Gleichung (4) nach @ und setzen den Ausdruck der Gleichung (8) ein, dann 
finden wir, wenn wir den Zusammenhang der Gleichung (3) beachten, nach einiger Zwischen- 
rechnung als Lésung der Gleichung (2) 


dy Wat ee 1—4(F) 
dy ae wily | (dy dex’ 
Ix (x) == a i ee (3) Va (x’) — 1 tz ; — (9) 


Der Einfachheit halber wollen wir zuweilen schreiben yH=y. 
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Bei der Lésung dieser Gleichung wollen wir so vorgehen, da wir zundchst den Wurzelausdruck 
in der Klammer gleich eins setzen und entsprechend der vorgegebenen Geschwindigkeitsvertei- 
lung W/V die nullte Naherung fiir dy/dx berechnen. Diese Werte setzen wir dann in die Wurzel 
ein und berechnen die nachste Naherung. Die Iterationen setzen wir solange fort, bis keine Unter- 
schiede mehr zwischen den einzelnen Schritten auftreten. Wie sich bei der Durchrechnung von 
Beispielen gezeigt hat, konvergiert dieses Verfahren verhiltnismaBig schnell. 

Nehmen wir fiir die Beschreibung des Profiltropfens nach Gleichung (4) eine endliche Reihe 
v<N an, dann 148t sich fiir die Gleichung (9) mittels der Gaufschen Quadratur die folgende 


Summenformel angeben: 
—7 1 dy(*) mn SS 
ine ih dp Cae = >? Grn fa (10) 
m=1 
Hierin bedeutet N die Anzahl der iiber den Bereich 0< y< = verteilten Punkte m, deren Lage 
sich entsprechend Gleichung (3) zu 


ered! mm Em z 1 COU 
Xm = gy C08 Sr baw. je ee 2! -++ cos Ww) (11) 


berechnet. Es gilt dabei m= 0 fiir die Hinterkante und m= N entsprechend fiir die Vorderkante. 
Der Ausdruck f,, = f (¢ = ist die an den Stellen m herrschende Geschwindigkeitsfunktion. 


Sie berechnet sich fiir den p-ten Naherungsschritt zu 


f2) = x8) fO4+ x@) —1, (12) 
wobei 
W 
fO= (F), 1 (13) 
und 


=| 4) as) 


am 


ist. a,, und a,,, sind konstante Zahlen. Die Formeln zu ihrer Berechnung sind im Abschnitt 6 
angegeben. Sie sind fiir N12 berechnet und in Tabelle1 mitgeteilt. 


abellemlewar ea... Nee 
m 1 2 3 4 5 6 | 8 9 10 11 
Em/! 0,98297 0,93302 0,85356 0,75000 0,62941 0,50000 0,37059 0,25000 0,14645 0,06699 0,01704 
am 0,12941 0,25000 0,35356 0,43302 0,48296 0,50000 0,48296 0,43302 0,35356 0,25000 0,12941 
n=0] 0,16383, 0 0,14226 0 0,10490 0 0,06177| 0 0,02441 0 0,00284 
1 0 0,20854 0 0,13414 0 0,08627 0 0,04264 0 0,01137 0 
2 —0,05588 0 0,26219 0 0,12805 0 0,06912 0 0,02649 0 0,00305 
3 0 —0,13109 0 0,30178 0 0,11785 0 0,05178 0 0,01324 0 
4 —0,01198 0 —0,20119 0 0,32237 0 0,10246 0 0,03452 0 0,00381 
F 0 —0,03431 0 —0,25915 0 0,32197 0 0,08237 0 0,01852 0 
ann } 6 —0,00578 0 —0,05893 | 0 —0,30040 0 0,30040 0 0,05893 0 0,00578 
| ( 0 —0,01852 0 —0,08237 0 —0,32197 0 0,25915 0 0,03431 0 
| 8 —0,00381 0 —0,03452 0 —0,10246 0 —0,32237 0 0.20119 0 0,01198 
9 0 —0,01324 0 —0,05178 0 —0,11785 0 —0,30178 0 0,13109 0 
10 —0,00305 0 —0,02649 0 —0,06912 0 —0,12805 0 —0;26219 0 0,05588 
li 0 —0,01137 0 —0,04264 0 —0,08627 0 —0,13414 0 —0,20854 0 
12 | —0,00284 0 —0,02441 0 —0,06177 0 —0,10490 0 —0,14226 0 —0,16383 


Es ist bemerkenswert, daB die Werte a, verschwinden, wenn |m— n| geradzahlig (0,2,4...). 
Diese Tatsache ist besonders wertvoll fiir die Durchfiihrung der bereits oben erwahnten Iterations- 


methode. 
Um die nullte Naherung zu berechnen, setzen wir x9) =1, was fiir f) den Wert der Glei- 


chung (13) ergibt, der unmittelbar aus der vorgegebenen Geschwindigkeitsverteilung entnommen 
werden kann. Diese Werte setzen wir in Gleichung(10) ein und bestimmen fiir die Punkte 
n=1,3... die Profilneigungen y/(. 

Mittels der Gleichungen (12), (13) und (14) errechnen wir dann fir dieselben Punkte 
m=1, 3... die Funktionswerte f, die dann, in Gleichung (10) eingesetzt, die erste Naherung 


fiir die Profilneigungen y‘() an den Stellen n =2, 4... ergeben. 
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Die zweite und folgenden Naherungen gehen dann entsprechend vor sich. Hat man gute 
Ubereinstimmung zwischen den letzten Naherungen erzielt, so empfiehlt es sich, noch die Werte 
an der Vorderkante yy und an der Hinterkante yo zu berechnen. Die erforderlichen Quadratur- 
konstanten sind ebenfalls in Tabelle 1 mitgeteilt. Kennt man jetzt den Verlauf y’(g), so 
ergibt eine Integration iiber p die gesuchten Profilordinaten y, die man lediglich an den n zu- 


geordneten Stellen , aufzutragen hat. Mit Riicksicht auf die SchlieBungsbedingung mu sein 


Die Durchfithrung des eben beschriebenen Verfahrens wird im Abschnitt 5 an Hand eines Bei- 
spiels gezeigt. | 

Der Vollstandigkeit halber sei noch eine Beziehung zur unmittelbaren Berechnung der Pro- 
filordinate abgeleitet. Durch Einsetzen der Gleichung (8) in die Gleichung (4) erhalten wir nach 
einiger Zwischenrechnung 


(t) il Pie Uaes sey ; 
we aa | P(e) sing! nee de (15) 


0 
Dieses Ergebnis erhilt man natiirlich auch, wenn man Gleichung (9) iiber x integriert. Zur nu- 
merischen Auswertung dieser Formel dient die Summenformel 


N—1 
ye eee eee 
= (9 =F Sas ee (16) 


Die Formel fiir 6,,,, ist in Abschnitt 6 und die Zahlenwerte fiir N= 12 sind in Tabelle2 mitgeteilt. 


Mabelle2anb777, INT 


m ul 2 3 4 5 6 a 8 9 10 11 


0,98297 | 0,93302 0,85356 0.75000 0,62941 0,50000 0.37059 0,25000 0.14645 | 0,06699 0,01704 
0 0 0 | 0 


n=0 0 0 0 0 0 0 0 
1] 0,02616 | 0,02372 | 0,01892 | 0,01699 | 0,01382 | 0,01111 | 0,00811 | 0,00552 | 0,00320 | 0,00148 | 0,00037 
2] 0,01228 | 0,06392 | 0.04742 | 0,03556 | 0,02968 | 0,02270 | 0,01689 | 0,01120 | 0,00660 | 0,00299 | 0,00076 
3 | 0,00692 | 0,03353 | 0.10052 | 0,06770 | 0,04777 | 0.03711 | 0,02631 | 0,01770 | 0,01016 | 0,00466 | 0,00117 
4| 0,00508 | 0,02053 | 0,05528 | 0,13038 | 0,08167 | 0,05399 | 0,03869 | 0,02483 | 0,01445 | 0,00647 | 0,00165 
5} 0,00370 | 0,01537 | 0,03497 | 0,07323 | 0.14979 | 0.08750 | 0.05354 | 0,03469 | 0.01926 | 0,00874 | 0.00217 
ee 6 | 0,00288 | 0,01135 | 0.02624 | 0,04675 | 0,08452 | 0,15652 | 0.08452 | 0,04675 | 0.02624 | 0.01135 | 0.00288 
] 7] 0.00217 | 0.00874 | 0,01926 | 0,03469 | 0,05354 | 0,08750 | 0,14979 | 0,07323 | 0,03497 | 0,01537 | 0.00370 
| 8] 9.00165 | 0.00647 | 0,01445 | 0,02483 | 0,03869 | 0.05399 | 0,08167 | 0,13038 | 0,05528 | 0,02053 | 0.00508 
9] 0,00117 | 0,00466 | 0,01016 | 0,01770 | 0,02631 | 0,03711 | 0,04777 | 0.06770 | 0.10052 | 0.03353 | 0.00692 
10 | 0,00076 | 0,00299 | 0.00660 | 0.01120 | 0,01689 | 0.02270 | 0.02968 | 0.03556 | 0.04742 | 0,06392 | 0.01228 
11] 0,00037 | 0,00148 | 0,00320 | 0,00552 | 0.00811 | 0.01111 | 0.01382 | 0,01699 | 0.01892 | 0,02372 | 0.02616 
2] 0 0 0 0 0 0 0 0 0 Eto “0 


3. Skelettprofil mit Anstellwinkel, Ist die Geschwindigkeitsverteilung auf Ober- und Unter- 
seite an gleichen Stellen der Profilsehne voneinander verschieden und folgt der Beziehung 


W=Viw»,, (17) 


wobei das obere Vorzeichen fiir die Oberseite und das untere entsprechend fiir die Unterseite gilt, 
so entspricht dieses einem Skelettprofil, d.h. einem Profil verschwindender Dicke! y(), das unter 
einem bestimmten Anstellwinkel « angeblasen wird. Damit an der Hinterkante glatter AbfluB 
ne d.h. die Kutta-Joukowskische AbfluBbedingung erfillt ist, mu® an der Hinterkante 
w,= 0 sein. 

Die Berechnung von w, geschieht mittels einer auf der Sehne angenommenen Wirbelverteilung. 
Aus der Bedingung, da®B an der Skelettlinie keine Normalgeschwindigkeiten herrschen diirfen, 


findet man fiir die Skelett-Tangenten dy/dx nach F. Riegels, Gleichungen (1), (17) und (21) der 
II. Mitteilung, 


y2 

dy ies 1 OD ip Ge. 

da (*) — & = [ROS (18) 
= 12 


1 Wir wollen zuweilen schreiben y(s) = y. 
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Es bedeuten hierin y die Ordinaten der Skelettlinie und « den Anstellwinkel, beide von der 
Sehne aus gemessen. 


Zur Vereinfachung der weiteren Rechnungen wollen wir voriibergehend von dem profilfesten 
Koordinatensystem auf das strémungsfeste System iibergehen. Dann gilt 


yt=y—ar+K,, (19) 
_wobei K, eine willkiirliche Konstante bedeutet. 


Analog zu Gleichung(4) machen wir fiir die gesuchte Skelettlinie den trigonometrischen 
Reihenansatz 


l 
DES Se RRA IE (20) 


Wir kénnen die Lésung der Gleichung (18) unmittelbar von F. Riegels iibernehmen, 


Gleichung (37) der IJ. Mitteilung: 


gsing = 3) va; (cosyg —1). (21) 
1 
Hierin bezeichnet 
WwW 
e=F=+(p-}) (22) 


die vorgegebene Geschwindigkeitsverteilung. 
Durch eine Fourier-Analyse kénnen wir die Koeffizienten ya} aus Gleichung (21) wie folgt be- 
rechnen: 


4 
* 


ri = = | g(y’) sing’ cosy y’ dg’. (23) 
6 
Diesen Ausdruck setzen wir in Gleichung (20) ein und erhalten nach einiger Zwischenrechnung 
4 1 i 4 : 7 t. 4 , 
7-= — | g(¢’) sing’ In [(1 — cos(p + ¢’)) (1 — cos(p — 9’) ] dp’ + Ky, (24) 
f 
wobei K, wieder eine willkiirliche Konstante bedeutet. 
Den Anstellwinkel, das ist der Winkel zwischen der Sehne (Verbindungsgerade zwischen 
Vorder- und Hinterkante) und der Anstrémrichtung, finden wir leicht zu 
x ve i Y" 
al) = © (a) pO) == | sto) sing’ In tg dg’. (25) 
f 


Ist uns jetzt y* bekannt, so kénnen wir wieder auf das profilfeste Koordinatensystem durch 
die Beziehung 


cosp 


y= y* — 5 (y*(x) + ¥*(0)) + SE (yx) — 9*(0)) 0) 


iibergehen. 

Nehmen wir fiir die Beschreibung der Skelettlinie, Gleichung (20), eine endliche Reihe y< N an, 
dann kénnen wir fiir Gleichung (24) mittels der Gaufschen Quadratur die folgende Summen- 
formel herleiten: 


N—1 
ee) xn &m + aNnlim(gsin¢) . (27) 
Seema 2 


Hierin gilt g,, = g(~ == Tv . Das zweite Glied tritt auf, wenn die Geschwindigkeitsverteilung an 


der Vorderkante gegen Unendlich geht. d%,, sind konstante Zahlen. Die Formeln zu ihrer Be- 


rechnung sind in Abschnitt 6 angegeben. 
Unter Beachtung der Gleichung (25) finden wir dann sofort fiir den Anstellwinkel die Summen- 


formel 
N—1 


a(s) — dm 8m —+ dnlim(g sing) . (28) 
grr 


m=1 
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Den Ubergang vom strémungsfesten System auf das profilfeste System nehmen wir nach 


Gleichung (26) vor und finden 
N—1 


x) = se (v tet T) £3 » dgtigete dan ea (gsing) . (29) 


Die Bedeutung der Konstanten d,, und d,,, ist in Abschnitt 6 angegeben. Sie sind fiir N= 12 be- 
rechnet und in Tabelle4 mitgeteilt. 
Da gm = + (4 — 1) unmittelbar aus der vorgegebenen Geschwindigkeitsverteilung ent- 


nommen wird, geht die Berechnung des Anstellwinkels sowie der Skelettordinaten in sehr ein- 
facher Weise vor sich. Die gefundenen Werte y®) sind dann an den n zugeordneten Stellen £, 


aufzutragen. Die Durchrechnung eines Beispiels wird im Abschnitt5 gezeigt. 
Der Vollstandigkeit halber sei noch die Summenformel zur Berechnung der Skelettneigung 


angegeben. 


N—1 
1 dy* mun ; : 
~* es —_ 
y= on = ae > Cmn &m — €Nn lim (g sing). (30) 
2 Ea grr 
Die Formeln fiir c,,, sind in Abschnitt 6, die Zahlenwerte fiir N= 12 in Tabelle3 mitgeteilt. 

Mabelleis meni N=12 

mm | 1 2 3 4 5 6 7 8 9 10 11 12 

E Alt 0,98297| 0,93302| 0,85356| 0,75000| 0,62941| 0,50000| 0,37059| 0,25000) 0,14645| 0,06699| 0,01704| 0 

n=0 0 0 0 0 0 0 0 0 0 0 0 0 
| 1 0 0,10795| 0 0,04009| 0 0,02233| 0 0,01274; 0 0,00589, 0 0,00549 

2 —0,10795| 0 0,18540| 0 0,06628| 0 0,03578| | 0 0,01873| | 0 0,00589| 0 
3 0 |-0,18540] 0 0,24640| 0 0,08333| 0 0,04227| 0 0,01873| 0 0,01726 

4 ~0,04009/ 0 |-0,24640) 0 0,28904| 0 0,09187| 0 0,04228/ 0 0,01274| 0 
5 0 |-0,06628| 0 |—0,28904] 0 0,31100; 0 0,09187| 0 0,03578| 0 0,03197 

Ee 6 —0,02233/ 0 |—0,08333) 0 |-0,31100] 0 0,31100; 0 0,08333| 0 0,02233 | 0 
eee 0 |-0,03578| 0 |—0,09187/ 0 |-0,31100/ 0 0,28904| 0 0,06628| 0 0,054.30 

8 ~0,01274| 0 |-0,04228) 0 |-0,09187; 0 |-0,28904| 0 0,24640; 0 0,04009| 0 
9 0 |-0,01873| 0  |~0,04227| 0 |-0,08333/ 0 |-0,24640| 0 0,18540. 0 0,10059 

10 -0,00589/ 0 |-0,01873; 0 |-0,03578| 0 |-0,06628; 0 |—0,18540] 0 0,10795 | 0 
11 0 |-0,00589/ 0 |—0,01274) 0 |-0,02233| 0 |—0,04009} 0 |—0,10795| 0 0,31649 

12 0 0 0 0 0 0 0 0 0 0 0 0 
Tabelle 4. dy .dmn N= 12 

1 

m | 1 2 3 4 5 6 7 8 9 10 11 12 

eit 0,98297| 0,93302| 0,85356| 0,75000| 0,62941| 0,50000| 0,37059| 0,25000| 0,14645| 0,06699| 0,01704| 0 
ae —0,04530 |—0,05399 |-0,05248 |—0,03932 |~0,02147| 0 0,02147| 0,03932| 0,05248| 0,05399] 0,04530| 0,15652 

n=0 0 0 0 0 0 0 0 0 0 0 0 0 
1 | 0,01007| 0,00978] 0,00300| 0,00427| 0,00148| 0,00204) 0,00029| 0,00050 |—0,00054 |—0,00038 |—0,00066 |—0,00216 
2 |-0,00925| 0,04079| 0,02319| 0,01342) 0,00965| 0,00590| 0,00314| 0,00072 —0,00109 |—0,00207 /—0,00227 |—0,00905 
3 |—0,01716|—0,00696| 0,07116| 0,04083 | 0,02199| 0,01528] 0,00695| 0,00248 |—0,00241 |—0,00418 | —0,00498 |0,01947 
4 |-0,01831|—0,00703| 0,02494| 0,10089 | 0,05325| 0,02813| 0,01531| 0,00466 |—0,00284 |—0,00703 |—0,00815 |0,03316 
5 ]-0,01763|—0,01022] 0,00550| 0,04412| 0,12007| 0,05977] 0,02659]| 0.01065 |—0,00304 |—0,00925 |—0,01178 |—0,04817 
any 6 |-0,01494|~0,01111| 0 0,01925 | 0,05577| 0,12778] 0,05577| 0,01925| 0  |—0,01111|—0,01494|—0,06389 
7 |—0,01178 |—0,00925 |—0,00304] 0,01065 | 0,02659] 0,05977| 0,12007| 0,04412| 0,00550 |—0,01022 !—0,01763 |—0,07757 
8 |—0,00815 |—0,00703 |—0,00284| 0,00466 | 0,01531| 0,02813| 0,05325| 0,10089| 0,02494 |—0,00703 |—0,01831 |—0,08872 
9 |—0,00498 |—0,00418 |—0,00241| 0,00248 | 0,00695) 0,01528] 0,02199| 0,04083| 0,07116| 0,00696!—0,01716 |—0,09303 
10 |—0,00227 |-0,00207|—0,00109| 0,00072 | 0,00314| 0,00590] 0,00965| 0,01342| 0,02319| 0,04079|—0,00925 |—0,09061 
11 {—0,00066 |—0,00038 |—0,00054 0,00050 | 0,00029| 0,00204| 0,00148] 0,00427| 0,00300| 0,00978| 0,01007|—0,06584 

12 0 0 0 0 0 0 0 0 0 0 0 0 


4. Unsymmetrisches Profil mit Anstellwinkel. Wir wollen in diesem Abschnitt den allgemei- 
nen Fall des unsymmetrischen Profils mit Anstellwinkel behandeln. 

Die Geschwindigkeitsverteilung ist an gleichen Stellen der Profilsehne auf Ober- und Unter- 
seite verschieden. Wir wollen dafir schreiben 


W=V+w,+u, +aw,. (31) 
Uber die Bedeutung der Zusatzgeschwindigkeiten w, und w, wurde bereits in den Abschnitten 2 
und 3 berichtet. Die Geschwindigkeitskomponente w, entsteht dadurch, da® der endlich dicke 
Profilkérper unter dem Anstellwinkel a=1 angeblasen wird. Das obere Vorzeichen gilt wieder 
fiir die Profiloberseite, das untere entsprechend fiir die Unterseite. Es sei bemerkt, daB vom 


physikalischen Standpunkt aus die Geschwindigkeitsverteilung nur endlich groBe Werte haben 
kann. 
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Ist uns jetzt die Geschwindigkeitsverteilung auf der Oberseite mit W, und auf der Unterseite 
mit W,, gegeben, so berechnen sich die einzelnen Geschwindigkeitskomponenten wie folgt: 


V+w, => (W,+W,) und wy + 01105 = > (W, — W,). (32) 


Um die Geschwindigkeit, die an der Oberflache herrscht, auf die Geschwindigkeit auf der 
Sehne umzurechnen, ist nach Riegels die Oberflichengeschwindigkeit mit dem Faktor 


oul 


zu multiplizieren. Bei einem unsymmetrischen Profilkérper sind die Profiltangenten dy/dx an 
gleichen Stellen der Sehne auf Ober- und Unterseite verschieden. Da der Faktor x nur von 
Bedeutung ist an den Stellen groBer Profilkriimmung (Vorderkante), die im wesentlichen von der 
Tropfenform abhangt, wollen wir den SkelettlinieneinfluB vernachlissigen und machen die 
folgende Vereinfachung: 


“= eee (33) 


Hierin bedeuten y( die Ordinaten des Profiltropfens. 


Auf Grund der gemachten Vereinfachung kénnen wir die Tropfenform sofort berechnen, wenn 
wir in Abschnitt 2 fiir die Geschwindigkeitsverteilung 


bzw. fiir Gleichung (6) 
f= oo ee ih (34) 


setzen. Das dort beschriebene Iterationsverfahren, Gleichung(10), liefert die Profiltangenten 
yt’ in Abhangigkeit vom Winkel g, was durch Integration iiber g die gesuchten Ordinaten des 
Profiltropfens y() ergibt. 

Setzen wir die Komponente w,, die ja nur von der Profiltropfenform abhangt, als bekannt 
voraus, so kénnen wir die Skelettform und den Anstellwinkel nach Abschnitt 3 berechnen, wenn 
wir fiir die dort benétigte Geschwindigkeitsverteilung 


w= 1 + (2) (g (Me — M) — a4) 


bzw. fiir Gleichung (22) 


Ag (35) 


setzen. Es gilt 
us 
Ve 


Ag=— ox (36) 


Die Rechnung nach Abschnitt 3 fiihren wir durch, indem wir zunadchst Ag —0 setzen, d.h. mit 


. : Wo— Wu J dy(*) 
5 7 bzw. Be ee Vy (a) PEE (x) |. 


7m 


(35a) 


Den Faktor x entnehmen wir unmittelbar der Rechnung fiir den Profiltropfen, Abschnitt 2. Die 
Anwendung der Gleichungen (28) und (29) liefert dann den Anstellwinkel x) und die Ordinaten 
der Skelettlinie y(), beide von der Profilsehne aus gemessen. 

Die bisher vernachlassigte Geschwindigkeitskomponente w, berechnet sich nach F’. Riegels, 


Gleichung (33) der II. Mitteilung, zu 


1 
Ws; aa 1 i 1—2x 
Vv (x) = yi+(er az ¥Yl+2x 
U2 


2 
zh w : A Ba (37) 
L ey (7) X—X 


l 


Unter Beriicksichtigung des Fourier-Reihenansatzes der Gleichung (4) tibernehmen wir die L6- 
Dyas 
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sung unmittelbar von F’ Riegels, Gleichung (38) der I. Mitteilung: 
w, _ 1—cosp 38 
nat = > cos 1p (38) 


sing 


mit p,=—= > 00, und By = 2 & [Ar + 4 — IV Barty TEM bara - 
1 | 


Setzen wir jetzt Gleichung (38) in Gleichung (36) und diese wiederum in Gleichung (23) ein, dann 
erhalten wir nach Ausfiihrung der Integration 


vAat = — 2 (28, —Bi1—Bry1) fir 9+] 


und 
Aat = — 5 (2, — 2B, — fa) « (39) 


Wir differentiieren Gleichung (20) nach y und setzen die Werte der Gleichung (39) ein, dann er- 
halten wir nach einiger Zwischenrechnung 


: A - = 5 ig ie Pace B,, sing] 
= aces (vcosvp + Sing) : | @) 
2 sing = sing 
Die Integration iiber ¢ liefert 
Ay* = —otg $5 3b, sinvp =— arg $9. (41) 


Die letzte Beziehung folgt aus Gleichung (4). Wir kénnen also den EinfluB, der durch die An- 
stellung des Profiltropfens hervorgerufen wird in sehr einfacher Weise als eine zusatzlich Skelett- 
linienanderung deuten. 


Die entsprechende Anstellwinkelanderung ergibt sich nach Gleichung (25) zu 
eae ati (tg 2 ) eee (42) 


l L dp 


gon 
Den gesamten Anstellwinkel erhalten wir aus der Summe 
o = al) + Ao. 


Setzen wir Gleichung (42) ein, dann wird 


al’) AO re 
eae 5 te i) 
Bi 


Um die strémungsfesten Ordinaten der Gleichung (41) auf die profilfesten Ordinaten Ay) um- 
zurechnen, benutzen wir Gleichung (26) und erhalten fiir die Skelettliniendnderung 


Tabelle 5. pg, N=12. Cs = yOu Oe 
Ay x (1 — cosq) ane + ae (ou)yayc (44) 
An den Punkten n ergibt sich 
es dee VS aes 5 ROC TE NT e, Ey, 

1 | 0,98297 | 0,13165 | 0,03 407 AVS A (¢ ie ®) wise « (Pay? + my - (45) 
2 | 0,93302 | 0,26795 | 0,13397 ’ 

3 | 0,85 356 | 0.41421 | 0.29289 Die Formeln fiir die Konstanten Pn und q, sind in Abschnitt 6 
5 ee at tee tes 0.50000 angegeben. Sie sind fiir N—12 berechnet und in Tabelle 5 
6 0.50 ae Oe ae a mitgeteilt. Ist also der Anstellwinkel nach Gleichung (43) 
7 | 0,37059 | 1,30323 | 1,25882 berechnet, so kann auch sofort Ay) berechnet werden. 

8 | 0,25000 | 1,73205 | 1,50 000 Wir setzen nunmehr die drei Anteile y() (Abschnitt 2), 
= Hee oe ae se a ie y (Abschnitt3) und Ay) (Abschnitt 4) zusammen und er- 
il 0.01 704 7°59 575 196 593 halten die gesuchten Ordinaten wie folgt: 


¥y = y(4 +y() ol. Ay) : (46) 

AbschlieBend sei noch darauf hingewiesen, da die Ergebnisse aus den Integralformeln mit 
denen aus den Summenformeln identisch sind, wenn sich die Profiltropfenform durch eine sin- 
Reihe und die Skelettlinienform durch eine cos-Reihe bis zum Gliede N (= 12) darstellen lassen. 
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5. Beispiele. Um die Brauchbarkeit der in den Abschnitten 2, 3 und 4 abgeleiteten Formeln 
zu zeigen, haben wir zunachst fiir einige bekannte Falle die Rechnungen durchgefithrt. Dabei 
sind sowohl die Integral- als auch die Summenformeln gepriift. Im Anschlu8 daran wird fiir eine 
willkirlich gewahlte Geschwindigkeitsverteilung die Profilform ermittelt. 

a) Ellipsenprofil ohne Anstellwinkel. (Abschnitt 2.) Die Geschwindigkeitsverteilung 
ist gegeben durch 

Mg i aaah ae as (47) 
i 1 +a? ctg?p 
Mittels Gleichung (10) haben wir fiir a= 1/4 die Profilneigungen des gesuchten Tropfens berechnet. 


: : ¥ : : Ay, 1 I 
Bereits nach der dritten Naherung ergeben sich Zahlenwerte die von den Werten y, = = cos —~ 


8 N 
nur noch geringfiigig abweichen. Durch Integration finden wir 
y(t) (Peale 
-i5 — 2 sin , (48) 


was die Ordinaten eines Ellipsenprofils vom Dickenverhaltnis t/l = a—=1/4 sind. Da der Wurzel- 
ausdruck in Gleichung (47) den Faktor x nach Gleichung (7) darstellt, kénnen wir fiir die Geschwin- 
digkeitsfunktion nach Gleichung (6) schreiben f—a, was, in Gleichung (15) eingesetzt, ebenfalls 
das Ergebnis der Gleichung (48) ergibt. 

b) Ebene Platte mit Anstellwinkel. (Abschnitt 3.) Die Geschwindigkeitsverteilung ist 
gegeben durch 


a Q 
Vv == || re atg -> 6 (49) 


Mit g = atg f nach Gleichung (22) liefern sowohl die Integralformeln (24), (25) und (26) als auch die 


Summenformeln (28) und (29) das richtige Ergebnis 
ev=e ud y)=—0. (50) 
c) Parabel-(Kreisbogen)-Skelett ohne Anstellwinkel. (Abschnitt3.) Die Ge- 
schwindigkeitsverteilung ist gegeben durch 
= l+asing bzw. g=—asing. (51) 


Auch bei diesem Beispiel liefern Integral- und Summenformeln das richtige Ergebnis 


of) = 0 und = = + sin’. (52) 
Es stellt a/4 = f/l die Wélbungshéhe dar. 

d) Skelettprofil mit konstanter Geschwindigkeitsverteilung. (Abschnitt 3.) 
Wir wollen den Fall betrachten, da®B auf Ober- und Unterseite jeweils konstante Geschwindig- 
keiten herrschen. An der Hinterkante muf sich die Geschwindigkeitsverteilung unstetig andern, 
um die AbfluBbedingung zu erfiillen, die fordert, daB dort Anstrémgeschwindigkeit herrscht. 
Ks ist 

W 
Tithe 
Die Anwendung der Gleichungen (25) und (28) ergibt a) —0. Gleichung (26) in Verbindung mit 
Gleichung (24) liefert in bekannter Weise 


a “bawe. g=> ea. (53) 


— cos Pp 


O82 (14 cos gyn S22 (54) 


ieee £ (1 — cos ¢) In } 9 


l aot 


Die Werte nach der Summenformel (29) weichen gegeniiber den exakten etwas ab. Siehe Abb. 1 
und Tabelle 6. 
e) Symmetrisches Joukowski-Profil mit Anstellwinkel. (Abschnitt 2, 3 und 4.) 
Die Geschwindigkeitsverteilung ist gegeben durch 
w 1+ ¢(1—2cosq) + a(1— 2¢ cos} tet an 
ee 1 7 (eye eee: 


sin ~ 
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Wenn 


ae = = sin p (1 — cos 7) (56) 


die gesuchte Profilform ist, dann stellt der Wurzelausdruck in Gleichung (55) den Faktor % nach 
Gleichung (7) dar, womit man die Geschwindigkeitsfunktion nach Gleichung (34) zu 

f=e(1— 2 cos¢) (57) 
erhalt. Setzen wir fin Gleichung (15) ein, dann finden wir das angenommene Ergebnis der Glei- 
chung (56) bestatigt. Die Profilkontur entspricht dem symmetrischen Joukowski-Profil. Fiir 
Gleichung (35a) ergibt sich 


g=atg + (1—2¢cosg)=a (1+ 26) tg 4 — 2e sing ; (58) 
Hierfiir konnen wir das Ergebnis sofort aus den Beispielen b und c entnehmen. Und zwar ist 
a)—= (1+ 2¢)a und ae = — > sin’pa. (59) 
(s) 
02 Tabelle 6. Zum Beispiel 5d. 
‘exakt yt 
n QP E/l : 
exakt ummen 
fi © nach Summentorme/ | formel 
E | a 0 0 0 
7 8 | 22/3 | 1/4 | 0,1790a | 0,1838a 
0 02 o4 06 08 10 6 | 2/2 1/2 | 0,2206a | 0,2257a 


Abb. 1. Skelettprofil mit konstanter Geschwindigkeitsverteilung. 
(Beispiel 5 d.) 


Beriicksichtigen wir jetzt die Formeln des Abschnittes 4, dann wird fiir die Skelettlinienanderung 
nach Gleichung (44) 


Am =~ sintpa. (60) 


Addieren wir diesen Wert zu dem der Gleichung (59), so ergeben sich die Ordinaten der Skelett- 
linie zu null, d.h. wir haben es mit einem symmetrischen Profil zu tun. Fiir den Anstellwinkel 
finden wir nach Gleichung (43) 

C—@ (61) 


in Ubereinstimmung mit der exakten Lésung. 


f) Profil mit willkirlich gewahlter Geschwindigkeitsverteilung. (Abschnitt 2, 
3und 4). An Hand einer willkirlich gewahlten Geschwindigkeitsverteilung wollen wir jetzt den 
Gang einer vollstandigen Rechnung kurz erlautern. 


Die gewahlte Geschwindigkeitsverteilung ist in Abbildung2 dargestellt. Auf Ober- bzw. 
. . . . . . E a 
Unterseite hat die Geschwindigkeit im Bereich 0,1< wes 0,6 konstante Werte, und zwar W,/V=1,3 


c 
bzw. W,,/V=1,0. Im Bereich 0,6< ; <1,0 fallt die Geschwindigkeit linear auf den Wert an der 


Hinterkante W/V= 0,8 ab. Dem Verlauf im Bereich dene <0,1 liegt kein besonderer analytischer 


Ansatz zugrunde. In Tabelle 7 sind die Werte W/V fiir die Stellen m angegeben. 


Wir berechnen zunichst die Werte = ame und 1 Wo—We 


V 5 a . Die beiden Kurven sind in 


Abbildung 2 gestrichelt eingezeichnet. 
Aussehend von Gleichung (34) bestimmen wir gema8 Gleichung (13) die nullte Naherung fiir 
f, d.h. f. Die Anwendung der Gleichung (10) liefert dann die nullte Naherung fiir die Neigungen 


aes gesuchten Profiltropfens y‘. Wir bestimmen diese fiir die Stellen n=1, 3,...11. Aus 
y berechnen wir nach Gleichung (14) den Faktor #1), und damit nach Gleichung ‘12) die Ge- 
schwindigkeitsfunktion f(), Eine nochmalige Anwendung der Gleichung (10) ergibt y’@, die wir 
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fir n=2,4...10 berechnen. Nun bestimmen wir x) und f@) und berechnen y’@) (n=1,3... 11). 


Die weiteren Naherungsschritte werden entsprechend durchgefiihrt. Beim fiinften und auch letzten 
Schritt haben wir noch die Werte an der Nase yn und an der Hinterkante y, ermittelt. Wie aus 
der Zusammenstellung hervorgeht, konvergiert das Verfahren sehr schnell. Bereits nach der 
dritten Naherung ist nahezu der Endwert erreicht. Der Vollstandigkeit halber haben wir in der 
Tabelle7 noch die Werte fiir %, und f,, nach der letzten Naherung angegeben. Durch graphische 
Integration der Werte y’ iiber m und Auftragung iiber €/l erhalten wir die Ordinaten des Profil- 
tropfens y(), wie sie in der Tabelle und in Abbildung3 mitgeteilt sind. Die gréBte Dicke betragt 
t/l= 0,124 und befindet sich an der Stelle E,/I—0,4. Im Bereich der Hinterkante ergeben sich 
schwache negative Werte, was bedeutet, da die vorgegebene Geschwindigkeitsverteilung in 
diesem Bereich nicht erreichbar ist. Wir andern daher die Kontur an der Hinterkante etwas ab 
so daB nur positive Werte fiir y() auftreten. 


@ he 
i SS < Endwert 


Vie °Naherung . 


[ee \ 


x 
A skelettinie y! 


G04 


0 b2 OY G6 G8 S10 


a =104° 


Abb.2 ee 


Abb. 2. Geschwindigkeitsverteilung zum Beispiel 5f. Abb. 4. Profilkérper zum Beispiel 5f. 


Interessant ist noch das Ergebnis, das sich ergibt, wenn wir nur die nullte Naherung fiir f be- 
riicksichtigen. Mit Hilfe von Gleichung(i6) finden wir unmittelbar die Profilordinaten der 
nullten Naherung, wenn wir setzen f,, =f). Wie aus der Tabelle und Abbildung 3 ersichtlich, ist 
der Unterschied gegeniiber der letzten Naherung nicht sehr gro, so daf fiir rohe Uberschlags- 
rechnungen mit Gleichung (16) gerechnet werden kann. Wenn die vorgegebene Geschwindigkeits- 
verteilung nicht gar zu ungewohnlich ist, sind die Ordinaten der nullten Naherung immer etwas 
kleiner als diejenigen der endgiiltigen Rechnung. 

Ausgehend von der Gleichung (35a) berechnen wir g,. Den Faktor x, entnehmen wir der 
letzten Naherung fiir die Berechnung des Profiltropfens. Wir setzen g,,in Gleichung (29) ein und 
erhalten die Ordinaten der Skelettlinie y). Siehe Tabelle7 und Abb.3. Die gréBte Wélbungs- 
hohe betragt f/J—= 0,036 und befindet sich an der Stelle &-/l—0,5. Gleichung (28) liefert den An- 
stellwinkel zu o)=0,0205 (=) 1,17°. 


Den endgiiltigen Anstellwinkel erhalten wir nach Gleichung (43) zu x0= a208 


aoaq = 010182 (=)1,04°. 
Die zusatzliche Skelettlinienanderung Ay() bestimmen wir aus Gleichung (45). Wie aus Tabelle 7 
hervorgeht, ist diese verschwindend klein. Gema® Gleichung (46) haben wir die Ordinaten 
auf Ober- und Unterseite bestimmt. Abb.4 zeigt den gesuchten unsymmetrischen Profilkérper. 
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Tabelle 7. Zum Beispiel 5f. (Vgl. Abbildung 2 bis 4) 
n,m 0 1 2 3 4 5 6 7 8 9 10 11 12 
ra 6 Wh UE ee al ee ee 
WIV 0,8000 | 0,8213| 0,8837| 0,9831| 1,1125| 1,2632] 1,3000| 1,0300] 1,3000| 1,3000] 1,2870| 1,1050| 0,0000 
w,,|V | 0,8000| 0,8085| 0,8335 | 0,8728 | 0,9250} 0,9853| 1,0000| 1,0000| 1,0000| 1,0000) 0,9900 | 0,8500 | 0,0000 
fO 0,1851 | —0,1414 |—0,0721 | 0,0187| 0,1243] 0,1500| 0,1500| 0,1500| 0,1500] 0,1385 |—0,0225 | — 
ine —  |_0,1850 | —0,1366 |—0,0603 | 0,0352 | 0,1239} 0,1509| 0,1503| 0,1521 | 0,1578| 0,1723| 0,1085| — 
mm at 1,0020| 1,0056| 1,0127| 1,0162| 1,0077] 1,0008| 1,0003| 1,0018| 1,0068| 1,0297| 1,1340| — 
p=9 |= 010001 ENA! (010632)|)) = !||=0'0036i ieee == en =0:0439)1| ear=—= tine ==0)00 110 eae 
3: | 1 ae as 0,0259| 0,0791 | — — 020195) kee 00264, |e | 0106345 eee a 
Yn 2 SE 000231 ee 0.0567 — ial) SS ennneys|| 20 |Snoma|) = oo) — 
3 = aoe 0:02653|0 ee 0,0784| — 00198'|/ 25) ||—0;0263) |) =— Fl|—0;0614 a = 
4 | — |-0,0023) — | 0,0565| — | 0,0599| — |-0,0036| — |-0,0411} — _ |-0,0693; — 
5) [00m | = 0,0264 | 0,0784| 9 — 0.0198; — |-0,0262) —  |—0,0616| — | -0,0622 
i) 0,0000 | -0,0021 | 0,0010 | 0,0120 | 0,0306 | 0,0500| 0,0594| 0,0615 | 0,0574 | 0,0486 | 0,0352 | 0,0178 | 0,0000 
7) 0,0000 | -0,0031 | 90014 | 0,0082| 0,0254| 0,0450| 0,0560| 0,0581 | 0,0539 | 0,0444 | 0,0302 | 0,0124| 0,0000 
8m = 0,0064 | 0,0252| 0,0558| 0,0952| 0,1400| 0,1501] 0,1500] 0,1503| 0,1510| 0,1529| 0,1446} — 
~(s) 
- 0,0000 | 0,0013| 0,0055 | 0,0130| 0,0227} 0,0319| 0,0358| 0,0346| 0,0296) 0,0220| 0,0131| 0,0046| — 
(Ss 
47, 0,0000| 0,0001 | 0,0002| 0,0002| 0,0002} 0,0001 | 0,0001 | 0,0000 | —0,0001 | —0,0002 | —0,0003 | —0,0002 | 0,0000 
Vo 0,0000 |—0,0007 | 0,0066| 0,0252| 0,0535 | 0,0820 | 0,0953| 0,0961| 0,0869| 0,0698| 0,0480| 0,0222] 0,0000 
z 0,0000 | 0,0035 | 0,0046 | 0,0012 | —0,0077 | —0,0180 | —0,0235 | —0,0269 | —0,0279 | —0,0274 | —0,0224 | -0,0134 | 0,0000 
u 
| 
6. Formelzusammenstellung. In den Abschnitten 2 und 3 wurden der klareren Darstellung 
wegen die Formeln fiir die Quadraturkonstanten nicht angegeben. Mit Hilfe der Besselschen 


Formeln (vgl. Hiitte, 27. Aufl., S.208) lassen sich fiir eine endliche sin-Reihe die Fourier-Koeffi- 
zienten aus Gleichung (5) wie folgt berechnen: : 


Net 
2 : um. 
vb, = ene oe eis (v< N). (62) 
Entsprechend gilt fiir die endliche cos-Reihe aus Gleichung (21): 
N—1 
peat m (—1)’,. ; 
yay ve cos » | sin [> 8m + ay melts sing) (vy< N—1l), 
N-1 (63) 
1 m I Say yee : 
Nice een sin ar em ! Se pe ale 


“i 
Gleichung (11) gegeben. Das zweite Glied in Gleichung (63) tritt auf, wenn die Funktion g fiir 
a (Vorderkante) gegen unendlich geht. Werden die Gleichungen (62) und (63) sinngema® in 
die Beziehungen der Abschnitte2 und 3 eingesetzt, dann ergeben sich daraus nach einigen 


Zwischenrechnungen die Quadraturkonstanten in den Gleichungen (10), (16), (27), (28), (29) 
und (30) wie folgt: 


Dabei ist f,, und: f= AG = Wr) . Die Lage der Punkte in x-Richtung ist durch 


@mn| _ 1—(—1)"*"_. am[{_ a(m—n) x (m+n) 
cet LN Oe oy, (ct ONGhe Go oe iP 
Pu (ae aC 
Ban Can — ’ CNn aN tg oN . (64) 
= SOUP 
a! (m—n) a (m+n) 
ae = m 
mn aN > (cos » WV cos V i ). (65) 
- wztem 
OU ar fe On fn) x ah ) ( 1) ae 
d* > = m+n —— A) 
mn IN we - (cosy N + cos Vv N | N . 
in —il 
i Eyl be N-+n 
ped meas Ct) n , (—1) 
Nn = oN > 7 008 Ye | oN ‘ (66) 
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Mit Hilfe dieser Formeln berechnet sich weiterhin 


d,, = dn. a di, 0 (67) 
und 
1 1 
dnn = Inn —— (det + dino) + > cos ~* (diy —di) . (68) 
Fiir die in den Gleichungen (14) und (45) auftretenden Konstanten gilt : 
AD oe eure mun T% 
ay, = 2 sin Nig > Pn = tg oN , Qn = 1 —— COs aN 5 (69) 


Samtliche Konstanten wurden fiir N= 12 berechnet. Ihre Zahlenwerte sind in den Tabellen 1 
bis 5 mitgeteilt. 


7. Zusammenfassung. Es wird ein Verfahren zur Berechnung der Profilform (Profiltropfen, 
Skelettlinie und dazugehoriger Anstellwinkel) aus einer vorgegebenen Geschwindigkeitsverteilung 
an der Profiloberflache mitgeteilt. Das Problem wird zuriickgefiihrt auf die Auswertung einfacher 
Summenformeln, in denen jeweils nur die Produkte des Geschwindigkeitsbetrages an gewissen 


festen Stellen der Sehnenachse mit gewissen festen ein fiir allemal berechneten Konstanten 
auftreten. 


8. Nachtrag. Aus Abschnitt 4 folgt, daB man die zusatzlichen Strémungsverhiltnisse, die 
durch den angestellten Profiltropfen hervorgerufen werden, dadurch erfassen kann, daB man sich 
eine d4quivalente Anstellwinkelanderung, G1.(42), und eine aquivalente Skelettlinienanderung, 
G1. (44), berechnet. Von dieser Tatsache kann man auch recht vorteilhaft Gebrauch machen, 
wenn man fiir eine vorgegebene Profilform (Profiltropfen, Skelettklinie und Anstellwinkel) die 
Geschwindigkeitsverteilung zu bestimmen hat. 

Man ermittelt dann zuerst den rechnerischen Anstellwinkel 


a 2 dy) 
cal sini (x) } x 


und die rechnerischen Skelettordinaten 


be (t) . 
HO= y+ a (1—cosg)(2 + F(a) dew. 9 = y+ (Pm I+ Gy 9) 


sin p 


Bei Annahme kleiner Anstellwinkel wiirde sich die Geschwindigkeitsverteilung nach dem Riegel- 
schen Quadraturverfahren (vgl. hierzu E. Truckenbrodt1, S.326) dann folgendermafen berechnen: 


W 1 N—I1 N—1 
7 (Sn) = pase [lead + D>, Ann 2 yo + ie Cn 299) + 6 hs) 
Crean m=l1 \m=1 Hf 


Man kann zeigen, daB diese Formel der bisher angewandten identisch ist. 


(Eingegangen am 11. April 1951.) 
Anschrift des Verfassers: Dr.-Ing. E. Truckenbrodt, (20b) Braunschweig, KérnerstraBe 12. 


1 FuBnote2 auf S.1. 
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Tragwerksauslenkung unter bewegier Last. 
Von E., Pesiel. 


1. Aufgabe und Voraussetzungen. Uber das Thema dieser Arbeit liegen bereits viele Ver- 
éffentlichungen vor. Mit Riicksicht auf die zahlreichen Literaturhinweise in den Arbeiten 
von Steuding! und Odman? kann auf eine Erérterung friitherer Untersuchungen, die das gleiche 
Problem oder ahnliche Aufgaben behandeln, verzichtet werden, zumal in dem vorliegenden 
Aufsatz nicht an bekannte Darstellungsweisen angekniipft wird. Hier werden zunachst grund- 
sitzliche Fragen bei der Formulierung des Problems erértert und danach die Entwicklung eines 
Berechnungsverfahrens erlautert, das im wesentlichen eine Ausarbeitung des in der Dissertation® 
des Verfassers nur andeutungsweise erwahnten ,,zweiten“ Verfahrens darstellt. 

Wir untersuchen die lotrechte Auslenkung eines Tragwerks, iiber das ein Massenpunkt und 
eine Anzahl von lotrecht gerichteten Kraften, die auch zeitlich veradnderlich sein mégen, mit 
beliebiger Geschwindigkeit wandern. Dabei beschrinken wir unsere Betrachtung auf waage- 
recht gelagerte Balken-Tragwerke, deren elastisches Verhalten unter den Voraussetzungen der 
Navierschen Biegetheorie bestimmt werden kann. Nebeneinfliisse, wie z. B. rotatorische Trag- 
heit des Tragwerks, Wirkung der Schubkrafte, Reibungseinfliisse usw., werden vernachlassigt. 
Der Massenpunkt ist bei seiner Wanderung iiber das Tragwerk stets in unelastischer Berithrung 
mit der Fahrbahn des Tragwerks. 

Es ist grundsatzlich méglich, mit Hilfe der vorgelegten Methode auch Tragwerke von kompli- 
zierterem Aufbau (z.B. Hangebriicken) zu behandeln und die hier vernachlassigten Nebenein- 
fliisse zu beriicksichtigen. 


2. Formulierung des Problems. Der wandernde Massenpunkt (Masse M) habe die Stelle i 
(Laststellung in 7) erreicht. Die zu diesem Zeitpunkt t; vorhandene Auslenkung des Tragers 
driicken wir wie folgt aus: 


y= 2 dni (ti) Pol) (2.1) 


Zeitfunktion Formfunktion 


Damit sind bei Laststellung in i folgende Krafte am Trager wirksam: 
a) Die Tragheitskrafte der Tragermasse m(x) infolge Auslenkung des Tragers#: 


— fm S Gni Pni dx. 
B n=1 
b) Die Tragheitskrafte des wandernden Massenpunktes M infolge Auslenkung des Tragers: 
— M = (ni Pri a= 2 CQni Pri ale C Gri Pait Qni Pri %) x, : 


Hier wird auch die Formfunktion Pni nach t differenziert, da die Koordinate x; des wandernden 
Massenpunktes von t abhangt (dx,;/dt= ¢; = Lastgeschwindigkeit). 

c) Die Tragheitskrafte des wandernden Massenpunktes M infolge vorhandener Abwei- 
chungen p(x) der Fahrbahn von der Waagerechten (z. B. Unebenheiten, Uberhéhung usw.; y’ (x) 
wird als stetig in B einschlieBlich seiner Grenzen vorausgesetzt) : 


— M(ey" + vy’ X) ag : 


1 H. Steuding, Ing.-Arch. 5 (1934) S. 275. 


29S. eA: Odman, Differential Equation for Calculation of Vibrations Produced in Load-bearing 
Structures by Moving Loads. Preliminary Publications. International Association for Bridge and Struc- 
tural Engineering. 3rd Congress Liége 13.—18. 9. 1948. 

° E. Pestel, Kin strenges allgemeines Verfah 
den Lasten. Dissertation Hannover 1947. 


* Der Buchstabe ,,B‘‘ am Integralzeichen bedeutet, da® iiber den gesamten Bereich des Tragers zu 


Ee yeneen a Punkte bzw. Striche bedeuten partielle Ableitung nach der Zeit t bzw. nach der Orts- 


ren zur Berechnung von Tragerschwingungen unter wandern- 
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d) Die elastischen Riickstellkrafte : 
— f S dni (ES 4)!" dx. 
B n=1 


e) AuBere, auch zeitlich veranderliche Krifte (z.B. unausgeglichene Fliehkrafte von Fahr- 


zeugen usw. einschlieBlich des Gewichtes P des Massenpunktes M). Die Angriffspunkte Xp 


der wandernden Krafte stehen zu jedem Zeitpunkt in eindeutiger Beziehung zu der Stellung 
des Massenpunktes in x;: 


Le P(t). 


Wir variieren den Auslenkungszustand des Tragers bei Laststellung in i lings der Form- 
funktion @,;(x). Dabei leisten die vorgenannten Krafte folgende Arbeiten: 


= bani fm ni Zi Ini Pn dx ° 


— 6gni Pni (xi) MS (in Pri t+ 2C Ge Patt © Int Pas + Ini Und *) a, > 
— 04ni Pni (xi) M (?y" + y’X)s, » 
— O9ni [ei 3 (EJ pri)” dx, 

69ni [Pni(%i) P + = Pni (Xp) Pp(t)]- 


Bevor wir zur Aufstellung der Differentialgleichungen weiterschreiten, wollen wir aus der un- 
endlich groBen Menge zulassiger Formfunktionen 9, ; (x) zwei auswahlen, die fiir die Lésung des 
Problems bei dem gewahlten Ansatz (2.1) geeignet erscheinen. 


1. Eigenschwingungsformen des unbelasteten Tragers. Formfunktionen ¢,;, die 
den Eigenschwingungsformen des unbelasteten Tragers entsprechen, sind nicht mehr von der 
Laststellung i abhangig, so daB man unter Fortlassung des Index i einfach @, schreiben kann. 
Ferner bestehen folgende Orthogonalitatsbeziehungen : 


[mon,gidx=0, 
B 


[on (EJ pi)” dx = 0 | AUe hes 
B 


Weiterhin sollen die Formfunktionen wie folgt normiert sein: 
f Mt O, Ux == 1 
B 

Indem wir den Rayleighquotienten 


f gn (EJ pf)” dx 


Sm 2 dx 
B 


==) (HE, 


setzen, erhalten wir durch Gleichsetzung der virtuellen Arbeit der Tragheitskrafte und der 
elastischen Krafte einerseits und der auBeren Krafte andererseits folgendes System von unendlich 
vielen gekoppelten Differentialgleichungen : 


Gni + On Ini = Pn (%p) P, (t;) + @n (%i) X 
P 


co : (2.2) 

2 (pw (2y" LW Ba 3 Cini Got 2etni Gat & des Pot Gui Pa) a|| (M= 12,3, ----) 
n=1 

2. Eigenschwingungsformen des Systems bestehend aus Trager und wandern- 

dem Massenpunkt M. Wir wollen nun die Eigenschwingungsformen des Systems bestehend 

aus Trager und wanderndem Massenpunkt als Formfunktionen ¢,; wahlen. Dann sind die Form- 

funktionen nicht mehr von der Laststellung i unabhangig. Sie gehorchen folgenden Ortho- 
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gonalitatsbedingungen : 
Jn Pni Pi dX + M qi (%:) Pri (%i) = 0, | 
Jousseiy! dno | 
Wie zuvor seien die Formfunktionen normiert gemaB 


[m gis dx + Mop; (x) = 1. 
B 


ibe pp sell « 


Mit dem Rayleighquotienten 
fon (EJ pf)" dx 
fm g? dx + M y?, (xi) 
B 


= Men 


erhalten wir, ahnlich wie unter 1., folgendes System von gekoppelten Differentialgleichungen: 
dnt + 2s ai = 3 Pai (%p) Pp (8) + Gas (04) X | 
P 


x ee M (ety 4. i) te 2, (2 c ae Pri == ce Qni Pni Sig Qni Prix) «| (n= 1, 2, 3, ae -) s | 
n=1 
3. Erorterung der Gleichungssysteme (2.2) und (2.3). Zum Zwecke der Unterscheidung 
wollen wir den Beschleunigungsarten des wandernden Massenpunktes folgende Bezeichnungen 
zuordnen : 


Qni Pri (x;) > 2 ¢; Ini Pri (x;) > | c} Qni Pri (x;) : 


Zentripetalbeschleunigung 


(2.3) 


Vertikalbeschleunigung Coriolisbeschleunigung 


Die Wahl der Formfunktion in 2 2) hat zur Folge, daB® die Vertikalbeschleunigung des Massen- 
punktes, die er mit dem Tragerelement an der Laststelle i gemeinsam hat, gegeniiber der Coriolis- 
und der Zentripetalbeschleunigung eine ausgezeichnete Stellung einnimmt. Eine Entkopplung 
des Systems von Differentialgleichungen (2.3) wurde dadurch nicht erreicht. Diese scheitert, 
physikalisch gesehen, grundsatzlich daran, daB der wandernde Massenpunkt anderen Tragheits- 
kraften unterworfen ist als die nur in lotrechter Richtung ausgelenkte Tragermasse, wahrend 
vom mathematischen Standpunkt aus die tiefere Ursache darin zu suchen ist, da die das 
Problem beschreibende partielle Differentialgleichung eine Ubergangsbedingung an der Last- 
stelle besitzt', welche die Auflésung der partiellen Differentialgleichung in ein System von 
entkoppelten gewohnlichen Differentialgleichungen durch einen Bernoullischen Ansatz ent- 
sprechend (2.1) vereitelt. 

Die unterschiedliche Behandlung der Beschleunigungsarten des Massenpunktes infolge der 
speziellen Wahl der Formfunktion in 22) bringt es ferner mit sich, daB die nach Fortlassung 
der Kopplungsglieder verbleibende Restgleichung von (2.3) Koeffizienten und Stérungsglieder 
enthalt, die im Gegensatz zu denen in der entsprechenden Restgleichung von (2.2) aus Form- 
funktionen berechnet werden miissen, die mit der Laststellung i veradnderlich sind. Damit ist 
ein so erheblicher Rechenaufwand verbunden, daB die Lésung des vorgelegten Problems mit 
Hilfe von (2.3) nur dann ins Auge gefaft werden sollte, wenn der Einflu®B der Coriolis- und Zentri- 
petalbeschleunigungen gegeniiber dem der Vertikalbeschleunigung lediglich den Charakter einer 
Korrektur tragt, und man die Léisung bei n= 1 abbrechen kann. Ein derartiger Fall liegt 
z.B. vor, wenn die wandernde Lastmasse gro gegeniiber der Tragermasse ist, der Trager jedoch 
die im Brickenbau iibliche Biegesteifigkeit besitzt. 


4. Verfahren zur Lésung des Systems (2.2)2. An Stelle des Gliedes 


) (Fni Pn + 2 Ani Pn+ Ce Ini Pn ai Qni Pn XY ees (4.a) 


n=l + 


' Vel. Gleichung (6) der Arbeit Steuding. 
* Das Verfahren zur Lésung von (2.3) kann analog dem hier Vorgelegten entwickelt werden, so daB 
auf eine gesonderte Darstellung verzichtet werden kann. 


* Die mit Buchstaben bezeichneten Gleichungen werden lediglich bei der Entwicklung des Lésungs- 
verfahrens benutzt. 
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das in (2.2) erscheint, kénnen wir schreiben: 


2 ge 
dee (Ini Pn)xe: (4.b) 
n=1 
Der Zusammenhang dieser zwei gleichwertigen Schreibweisen lat sich auch physikalisch deuten. 
Bei (4.a) registriert man als Beobachter auf dem in beschleunigter Bewegung befindlichen 
Massenpunkt die verschiedenen auf diesen cinwirkenden Beschleunigungen, wahrend man bei (4.b) 
von einem ruhenden Standpunkt aus das Durchlaufen der Bahn des Massenpunktes (Lastweg) 
beobachtet und aus der Kriimmung des Lastweges (nicht zu verwechseln mit der Kriimmung 
des Tragers!) und der Horizontalgeschwindigkeit des Massenpunktes die auf ihn wirkende Be- 
schleunigung bestimmt. 
Der kiirzeren Schreibweise wegen fassen wir die aus der Fahrbahnabweichung von der Hori- 
zontalen resultierenden Triagheitskrafte 


M (cy! + py" %) x5 


zu einer bei bekannter Geschwindigkeit c; des Massenpunktes ecpeveres auBeren Kraft B; 
zusammen. Damit kénnen wir (2.2) wie folgt schreiben: 


4 3 oo d2 
dni +8 dai— >” gal Hp) Pp (tt) + Gn (a) [PB —M > Ee Gos Fae | (m= BPE ccain ea 
P n= 


Wir wollen das System (4.1) so lésen, da wir den wandernden Massenpunkt M schrittweise 
auf seinem Wege iiber den Trager verfolgen. Der Massenpunkt habe zur Zeit t;_, die Stelle i— 1 er- 
reicht und der Trager besitze an dieser Stelle die Vertikalauslenkung : 


Yi (Mi) = = Gn, i—1 Pn (%i_4) 
und die Vertikalauslenkungsgeschwindigkeit : 
Via (Xi) = 24 1-1 Pn (%:_4)- 


Damit finden wir, indem zunichst der Einflu8 der Massenkrafte des wandernden Massenpunktes 
unberiicksichtigt gelassen wird, die Auslenkung und Auslenkungsgeschwindigkeit des Tragers 
an der Stelle i zur Zeit t; auf elementarem Wege! zu 


sin @,, (t; — t;_1) ®,,(t;) | et (4.2) 


| 
On oO, (2) 


oo 


[yi (2), = 3 Pn (%) Moi COS Wp (t; — ti) + Qni—a 


n=1 


n 


[vi (*)], = S (xi) x 


(42x) a) (4.3) 
: : dt /t | dt /*; 
X | — Wn Yn, i—; 811 Wn (t; — 53) + Qn. i—1 COS Wn (t; aE t;_4) " 4 a 


(ti) = Sf Pp (tat) Gn (%p1+ 8) sin w, (44; — #) dt + 
p==0 


At; 
+ f Bt. +1) o (%1+ ¢#) sinw, (At, — t) dt, 
0 
At; 


Y(t) =) P Gn (%i_4 + €1) sin w, (At; — t) dt . 


Der Index 1 bei [y;(x;)], und [¥; (xi)], besagt, daB es sich um erste Naherungen handelt, da in UF At) 
der Einflu® der Massenkrafte des wandernden Massenpunktes unberiicksichtigt blieb. 

Wir wollen uns nun aufeinanderfolgende Punkte des Lastweges, gegeben durch ihre Ab- 
szissen x;_, bzw. x; und ihre Ordinaten y;_, bzw. y;, durch eine kubische Parabel verbunden 


1 K. Klotter, Einfiihrung in die Technische Schwingungslehre, 2. Aufl., Bd. 1, Salo tteeeeexiias 
Gottingen, Heidelberg 1951. 
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denken, deren Tangentenrichtung an der Stelle i—1 bzw. i mit der Richtung der Vektorsumme 
aus Horizontalgeschwindigkeit c;_, bzw. ¢; und der vertikalen Auslenkungsgeschwindigkeit 
Via (%i_4) baw. 9; (xi) zusammenfallt. Eine einfache Zwischenrechnung ergibt dann als erste 
Naherung der Lastwegsbeschleunigung in der Intervallmitte?: 
_ (ten? [yi(*)h Vi- (x;-1) | 
DiC), = aagel tc  qee en a 


as 3 i (%i) yg — Yi-1 Xi_1) 1 (ly: (%))1 Vi-y (Xi_1) 
+ @asea(3 B th ae rie | ‘ Ci_-4 }): 
2 


Dies vereinfacht sich, wenn wir im folgenden die Horizontalgeschwindigkeit des Massenpunktes 
als konstant betrachten, zu 


(¥; (x;)]; — [yi(x;)] a (%;-1) : (4.d) 
Mit diesem Wert fiir die Lastwegsbeschleunigung kénnen wir den Term fiir die Massenkrafte 
des wandernden Massenpunktes in (4.1) berechnen und somit zweite Naherungen fiir die Aus- 
lenkung und Auslenkungsgeschwindigkeit des Tragers an der Stelle 1 zur Zeit t; entsprechend 
den Gl. (4.2) und (4.3) ermitteln. 

Bevor wir den Iterationsproze8 einleiten, wollen wir in (4.d) in Hinblick auf (4.3) die Anteile 
abspalten, die bei der Iteration unverdndert bleiben. Dies sind zunachst die beiden ersten 
Glieder in (4.3) und y,;_, (x;_,), da sie nur GréBen enthalten, die durch die vorangegangenen 
Rechenschritte bis zur Stelle i—1 bereits endgiltig bestimmt worden sind; ferner gehért dazu 
der Term (d@,,/dt),,/w, in (4.3), der den Einflu8 von Kraften wiedergibt, die von der Tragerbe- 
schleunigung unabhangig sind. Wir fassen diese Anteile wie folgt zusammen: 


ed 5 5 in (X;_ \ d@,, 
2 Pn (i) i On Qn, 7-1 8in O, At; + Qn, ;-1 |cos w, At; — re | I ( di _|| 
oe ee (4.4) 


At; 
Fir das Restglied in (4.3) (d¥,/dt),,/@, schreiben wir nach Division durch At; die Abkiirzung 
[¥:(%)|, und kénnen nun fiir die Lastwegsbeschleunigung erster Naherung schreiben: 
[yi (mi) = Vit [yi (™) hr - 

Indem wir die mit Hilfe von (4.4) erhaltene Massenkraft erster Naherung mit dem Gewicht 
P= Mg des Massenpunktes zu einer ,,auBeren Kraft‘ (P;), erster Naherung 

fy Ve Le) 
P.), = P abu tri dL 
( i) it g g 


vereinigen und (P;); anstelle von P einsetzen, erhalten wir als Lastwegsbeschleunigung zweiter 
Naherung wahrend des i-ten Zeitintervalles: 


ae) ee Vie ee ee i ee cece ais 
[yi(%i) |, = Vi fe z g j bi (x) ], Lyi ils pa i ae 


Damit bestimmen wir, wie zuvor, eine ,,auBere Kraft‘ zweiter Naherung: 


o ch 2 oe 
P), = Py —Wi@h\ _ ply Wildl | idl Vi(y _ ied. 
(Pie \ ge J \ Suiie uke g ( g ) 
Durch fortgesetzte Iteration erhalten wir als ,,auBere Kraft‘* m-ter Naherung: 


r=0 r=1 


Fir m— o ergibt sich dann?: 


P (1 ze Z) 
einer Gk (4.5) 
yp = 


' Die Bestimmung der Lastwegsbeschleunigung als Funktion von t im Intervall At; bereitet grund- 
satzlich keine Schwierigkeiten; doch steht der geringe Gewinn an Rechengenauigkeit in keinem Verhiltnis 
zu der sich dann ergebenden Vergréerung des Rechenaufwandes. 

» Siehe FuBnote 3 von S. 378. 
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Indem wir P; anstelle von P in die Terme WY (t;) und (d¥,,/dt),, von (4.2) und (4.3) einsetzen, 
erhalten wir endgiiltige Werte fiir die Lastwegsordinate y;(x;) und die lotrechte Lastwegs- 
geschwindigkeit y;(x;) am Ende des i-ten Zeitintervalles und kénnen danach zur Berechnung 
der Stelle i-+-1 in der gleichen Weise fortschreiten. 


5. Anwendung. Fiir die Behandlung eines bestimmten Problems sind lediglich die Ermittlung 
der Formfunktion @p,, die Auswertung der Integrale ®, und YW, ihrer Ableitungen und die Auf- 
stellung des Rechenschemas unter Benutzung von (4.2), (4.3), (4.4) und (4.5) zu erledigen. Die 
Berechnung von Beispielen zeigte, da8 auBer in der Umgebung stirkerer Lastwegskriimmung 
die Schrittweite c; At; —= x; — x;_, ziemlich groB gewahlt werden kann. Einen Uberblick 
verschafft man sich am besten durch eine mit wenigen groBen Schritten durchgefiihrte Probe- 
rechnung, die nur die Grundeigenfrequenz und -schwingungsform (n= 1) beriicksichtigt. Die 
Untersuchung eines extremen Beispiels (weicher Trager, groBe Lastmasse im Verhaltnis “ur 
Tragermasse und hohe Lastgeschwindigkeit), die unter Benutzung einer jeweils verschiedenen 
Anzahl von Fouriergliedern durchgefiihrt wurde, erwies eine Uberschreitung von n=3 zur 
Erzielung ausreichender Genauigkeit als unnétig. Eigenschwingungsformen héherer Ordnung 
brauchen daher nur Beriicksichtigung zu finden, wenn pulsierende Krafte auf den Trager wirken, 
deren Frequenz in der naheren Umgebung einer der héheren Eigenfrequenzen des Tragers hegt. 
Die Anwendung des Verfahrens beantwortet gleichzeitig die fiir den Ingenieur wesentlichen 
Fragen nach der Tragerdurchbiegung und dem Querkraftsprung P; an der Laststelle. 


(Eingegangen am 20. April 1951.) 


Anschrift des Verfassers: Dozent Dr.-Ing. E. Pestel, Hannover, Nienburger Strafe 2. 
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Uber den Einflu8 der Absaugung auf die Lage der Umschlagstelle 
an Tragfliigelprofilen. 


Von H. Kriiger. 


1, Einleitung. Der Reibungswiderstand eines Tragfliigels 1aBt sich durch Laminarhalten der 
Grenzschicht betrachtlich vermindern. Ein Mittel hierzu ist die Verwendung von Profilen mit 
groBer Dickenriicklage. Bei diesen ,,Laminarprofilen‘‘ beginnt das den Umschlag laminar/turbu- 
lent begiinstigende Druckanstiegsgebiet weit stromabwarts. Eine andere Méglichkeit, die Grenz- 
schicht iiber groBe Laufstrecken laminar zu halten, bietet die Absaugung eines Teils der Grenz- 
schichtsubstanz durch die porése Kérperoberflache. Gleichzeitige Anwendung beider MaB- 
nahmen ermdglicht weitgehendes Laminarhalten bei geringen Absaugmengen und -leistungen. 
Wie Profilgestaltung und Absaugung die Umschlagspunktlage beeinflussen, soll im folgenden 
an Hand theoretischer Untersuchungsergebnisse aufgezeigt werden. 


2. Berechnungsgrundlagen, Ist fiir das an- t iy. 
gestrémte Profil (Abb. 1) der Verlauf der AuBen- lee 
geschwindigkeit U und der Sauggeschwindigkeit uy? 
— vy ttber der Laufstrecke s vorgegeben, so ]abt & 
sich die laminare Grenzschicht nach einem von 

H. Schlichting’ angegebenen Naherungsverfahren Z 
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Abb. 1. Erlauterungsskizze fiir die Grenzschicht Abb. 2. Die kritische Reynoldssche Zahl (U 6*/»))_;+ 
mit Absaugung. in Abhangigkeit vom Formwert 6*/3 = H. 


berechnen. Zu jedem Punkt der laminaren Laufstrecke erhalt man das ,,Geschwindigkeits- 
profil u(y) der Grenzschicht, damit Verdrangungsdicke 6* und Impulsverlustdicke 9 ?. Vom 
Quotienten 6*/} — H, dem ,,Formwert“‘, hangt nun im wesentlichen die kritische Reynoldssche 
Zahl der Reibungsschicht (U) 6*/v);,:, ab, bei deren Erreichen die laminare Grenzschicht 
instabil wird, also Umschlag eintreten kann. Abb. 2 zeigt diese Abhangigkeit, wie sie A. Ulrich 
fiir die exakten Geschwindigkeitsprofile der abgesaugten Staupunkts- und Plattenstrémung 
fand%, und wie sie sich nach einem geniigend genauen Naherungsverfahren? fiir nicht abgesaugte 
Kérper ergibt. Den nachstehend beschriebenen Untersuchungen wurde eine eindeutige Zuord- 


nung zwischen Formwert H und kritischer Re-Zahl (U-6*/v),,i, entsprechend der ausgezogenen 
Kurve zugrundegelest. 


3. Rechnungsgang: Ergebnisse fiir Profile iiblicher Form, Als Beispiel fiir Fliigelschnitte der 
tblichen Form betrachten wir ein symmetrisches Joukowsky-Profil von 15% Dicke und 25% 


* H. Schlichting: Ein Naherungsverfahren zur Berechnung der laminaren Reibungsschicht mit Absaugung. 
Ing.-Archiv, XVI. Bd., 3. und 4. Heft, 1948. 


2 Esist 6* U = | (U—u)dy (siehe Abb.1), &U2= i u(U —u) dy. 
e=0 ae 
* A. Ulrich: Theoretische Untersuchungen iiber die Widerstandsersparnis durch Laminarhaltung mit Ab- 
saugung. Schriften der Deutschen Akademie der Luftfahrtforschung, Band 8b, 1944, Heft 2. 


* H. Schlichting und A. Ulrich: Zur Berechnung des Umschlages laminar/turbulent. Bericht S 10 der Lilien- 
thal-Gesellschaft, 1940. 


XIX. Band 1951. 


Kriiger: Einflu8 der Absaugung auf die Lage der Umschlagstelle usw. 


385 


Dickenriicklage mit der potentialtheoretischen Geschwindigkeitsverteilung U(s) bei sym- 
metrischem Anstrom. Es liege homogene Absaugung vor, d.h. die Sauggeschwindigkeit — vy sei 


iiber der Lauflange konstant. 


—v/U,yU, lly liefert die Grenz- 
schichtrechnung die dimensionslose 
Verdrangungsdicke U-6§*/y itber der 
Lauflange s/l’. Die Auftragung in Abb. 3 
(diinn ausgezogene Kurven) fiir U,1/v 
= 10° und 108 zeigt, wie die Grenz- 
schicht mit steigender Absaugung diin- 
ner wird. Mit Hilfe der Rechenkurve 
von Abb. 2 kénnen weiterhin die kri- 
tischen Re-Zahlen der Grenzschicht 
(U-6*/v),,:: ttber der Laufstrecke auf- 
gezeichnet werden. Der Verlauf dieser 
,, Stabilitatsgrenzen‘*‘ — in Abb.3 stark 
ausgezogen — wird von der Absauge- 
starke wesentlich beeinfluBt. Bei kleinen 
Sauggeschwindigkeiten fallen die Sta- 
bilitatsgrenzen stetig ab; bei starker 
Absaugung bilden sie ein Minimum aus 
_und steigen zum Fliigelende hin an. 
Der Schnitt der Kurven fiir kri- 
tische Re-Zahl und Verdrangungsdicke 
gibt die Lage der Instabilitatspunkte 
oder theoretischen Umschlagspunkte 
auf der Profilkontur an. Der tatsich- 
liche Umschlag laminar/turbulent liegt 
gewohnlich weiter stromabwarts; spa- 
testens tritt er am Ablésungspunkt der 
laminaren Grenzschicht ein. Da die 
Verhaltnisse an dieser Stelle durch das 
benutzte Naherungsverfahren nur un- 
genau wiedergegeben werden, wurde 
fiir die Eintragung derAblésungspunkte 
in Abb. 3 ein exaktes Ablésungsprofil 
(nach Hartree) zugrunde- 
gelegt. (Die gestrichelte 
Fortsetzung des Grenz- 
schichtverlaufes iiber die 
Ablésungspunkte hinaus 
entspricht dem Nahe- 


Fiir verschiedene Re-Zahlen U_-l/y und Absaugestirken 
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Abb. 3. Ermittlung der Instabilitaitspunkte fiir ein symmetrisches 
Joukowsky-Profil mit homogener Absaugung. 
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Abb. 4. Lage der Instabilitatspunkte s/l/ beim symmetrischen Joukowsky-Profil 
in Abhingigkeit yon der Reynoldsschen Zahl Ux I/. 
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zur Profilnase. Bei groBen Absaugestarken treten neben dem Ablésungspunkt zwei Instabilitats- 


punkte auf, zwischen denen die Grenzschicht zum Umschlag neigt. 
Die Kurven der Abb. 5 zeigen die 


Ergebnisse der Abb. 4 in anderer Dar- 
stellung. Sie kénnen angenahert als 
ee optimale Absaugeverteilungen betrach- 
ae PSS tet werden, bei denen die Grenzschicht- 
s G20 id dicke knapp unter der Stabilitatsgrenze 
bleibt, die also mit kleinstem Mengen- 
aufwand den Umschlag verhindern. 
(DaB die Naherung gut ist, erwiesen 
genauere Ermittlungen, bei denen die 
Absaugung nicht jeweils konstant, 
sondern entsprechend der Optimalver- 
teilung zugrundegelegt war.) Kleine 
Re-Zahlen U,, l/y verlangen nach hinten 
zunehmende, groBe Re-Zahlen nach 
hinten fallende Absauggeschwindig- 
keiten. Bis zur Fliigelhinterkante lami- 
nare Grenzschichten erhalt man, wenn 
man der punktierten Verbindung der 
Ablésungspunkte gemaB nach hinten 
y 2 7) a a5 ie verstarkt absaugt. 
Abb. 5. Lage der Instabilitats- und Ablésungspunkte pee poms eae ee bret ene 
Joukowsky-Profil in Abhingigkeit von der Absaugestirke all Ue l/r. wit ee Moapeaheee oes an ie 
(Die Kurven entsprechen der ,,optimalen Absaugeverteilung‘.) Grenzschicht bis zur Hinterkante lami- 
: ; ed nar zu halten, bendtigt man beispiels- 
Nida bei U,lly= 10" ein cg = eg jv, l/y = 1,10 entsprechend c¢ = 3 5-104: Die Werte 
liegen in der GréBenordnung experimenteller Ergebnisse 2. e j 
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Abb 6. Lage der Instabilitaitspunkte s/l/ beim Laminarprofil in Abhangigkeit yon der Reynoldsschen Zahl U _ l/y 
oo l?: 


4, Ergebnisse fiir Laminarprofile. Fliigelschni i i 
, gelschnitte, bei denen sich der Druckabfall iiber ei 
weiten Bereich der Kontur erstreckt, erméglichen groBe laminare Laufstrecken. Wie widhistee 
1 


“1 Es ist of = | —%)/U VUx lv d(s/U’) = cg Vu Wr 
2 Holstein fand co = 1,1-10~4 bis 2,8 - 10-4, allerdings bei Punktabsaugung. 
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Absaugung bei diesen Profilen auf die Umschlagslage auswirkt, sei am Beispiel eines Laminar- 
profils mit extrem groBer Dickenriicklage gezeigt. Abb. 6 bezieht sich auf ein verallgemeinertes 
Joukowsky-Profil von 15% Dicke und 70% Dickenriicklage bei symmetrischem Anstrom und 
homogener Absaugung. Ein Vergleich mit Abb. 5 zeigt: die Stabilitatskurven des schwach 
abgesaugten Laminarprofils sind von gleichem Charakter wie die des itblichen Profils bei starker 
Absaugung. Auch hier entstehen in bestimmtem Re-Bereich drei Instabilitatspunkte, wobei 
zwischem vorderem und mittlerem und nach dem hinteren die Grenzschichtdicke den kritischen 
Wert ttberschreitet. In diesem Bereich wird der Umschlagspunkt sprungartig nach vorn verlegt. 
Verstarkte Absaugung verschiebt diesen Re-Zahl-Bereich nach oben, beeinfluBt jedoch die Lage 
des vorderen und hinteren Umschlagspunktes nur geringfiigig. Dementsprechend wird man sich 
bei diesen Tragfliigelprofilen auf die Beseitigung des vorderen Instabilititsbereichs beschranken, 
denn die geringe Verlagerung des hinteren Umschlagspunktes durch verstarktes Absaugen ist un- 
wirtschaftlich. Zur Laminarhaltung bis zum hinteren Umschlagspunkt kommt man mit kleinen 
Mengenbeiwerten der optimalen Absaugeverteilung aus: bei Re= 107 ist ein ch-Wert von 0,35 
entsprechend cg = 1,1-10~¢ erforderlich. 


5. Zusammenfassung. Die Untersuchungsergebnisse zeigen, wie man durch kontinuierliches 
Absaugen die im Druckanstiegsgebiet des Tragfliigels vorhandene Umschlagsgefahr beseitigen 
kann. Erméglicht werden dadurch weitgehend laminare Grenzschichten bei groBen Reynolds- 
schen Zahlen und, wie weitere Rechnungen! ergaben, auch bei hohen Auftriebsbeiwerten und 
Dickenverhaltnissen. Eine Dosierung der Sauggeschwindigkeit tiber der Profilkontur (gema8 der 
Optimalverteilung) bringt gegeniiber der homogenen Absaugung Mengen- und Leistungserspar- 
nis mit sich. Besonders geringe Absaugmengen benétigt man bei Profilen mit grofer Dicken- 
ricklage. : 


(Eingegangen am 30. April 1951.) 
Anschrift des Verfassers: Dipl.-Ing. H. Kriiger, Gottingen, Wiesenstr. 16, 


1 Braunschweiger Diplomarbeit des Verfassers. 
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Zur Veranschaulichung und Erweiterung der Theorie des Pendels 
mit oszillierendem Drehpunkt. 


Von W. Kucharski. 


1. Kurze Inhaltsangabe. Im Abschnitt 2 der folgenden Mitteilung werden nach Aufstellung 
der grundlegenden Gleichungen die Vorginge bei einem Pendel mit schnell oszillierendem Dreh- 
punkt bei ebener Bewegung der physikalischen Anschauung nahergebracht. Der physikalische 
Pffekt** der Oszillationen besteht hiernach darin, da® die Impulseinwirkungen auf die Pendel- 
masse in der einen Drehrichtung iiberwiegen; die Pendelmasse erhalt so im Ganzen einen Antrieb, 
der den spitzen Winkel zwischen der Pendelstange und der Oszillationsrichtung zu verkleinern 
sucht. Man kann auch hier von einer ,,Tendenz zum Parallelismus‘ sprechen. Diese Wirkung 
ist zu vergleichen mit derjenigen einer kontinuierlichen elastischen Riickstellkraft, welche dem 
Sinus des doppelten Winkels mit der Oszillationsrichtung proportional ist. 

Im Abschnitt 3 wird hieraus die Konsequenz gezogen, dafi die Pendelbewegung im allge- 
meinen aus einer langsamen Grundbewegung bestehen wird, die sich aus der genannten resul- 
tierenden Impulseinwirkung der Oszillationen und etwaigen eingepragten Kraften ergibt, und 
einer dariiber gelagerten Zusatzbewegung von hoher Frequenz und kleiner Amplitude. Im 
Rahmen der fiir diese Theorie ttblichen und gerechtfertigten Annaherungen wird dies verifiziert ; 
die Zusatzbewegung ergibt sich als eine Schwebung, deren Amplitude dem Sinus des Winkels 
zwischen Pendelachse und Oszillationsrichtung proportional ist. 


2. Die Grundgleichungen; einige Betrachtungen zur Veranschaulichung. a) Der Schwerpunkt A 
einer starren Masse m, ist laut Abb.1 in der x-Richtung reibungslos gefiihrt. Um A kann sich 
ein Pendelkérper ebenfalls reibungslos drehen, welcher der 
Kinfachheit wegen als eine ,,konzentrierte‘‘ Masse m mit einer 
masselos gedachten starren Stange von der Lange / dargestellt 
ist. Der Winkel zwischen der Pendelstange und der negativen 
x-Richtung sei y. Das System hat zwei Freiheitsgrade; offen- 
bar sind x und  ausreichende und unabhangige Koordinaten. 
Die Einheitsvektoren in der x-Richtung bzw. in der Richtung 
der bzw. senkrecht zur Pendelstange seien ¢, bzw. e, bzw. €,3; ex 
ist auch der Richtung nach konstant; e, und e, drehen sich 
mit der Winkelgeschwindigkeit m der Pendelstange, sind also 
mit der Zeit t veranderlich. 

Die kinetische Energie ist 


] : = he 
Abb. 1. Anordnung mit Bezeichnungen. E= 2 m, x2 + = m (xe, + ly 6)e5 (1) 
hieraus erhalt man mit e, +e, = sin @ leicht die Bewegungsgleichungen zu 
% dae 
(m, + m) x-+ ml tt (y sin gy) = P, (2) 
und 
ml (x sing + lp) = P,, (3) 


worin P, bzw. P, die eingepragten KraftgréBen fiir die Koordinaten x bzw. y bedeuten; P,, ist 
naturgema4B ein Moment. re 

Die Bewegung von A wird in den wichtigsten hier interessierenden Fallen als zwangslaufig 
angesehen werden. Gleichung (3) enthalt dann x als gegebene Funktion von t und ist bei ge- 
gebenem P,, das auch Null sein kann, fiir sich zu integrieren; Gleichung (2) dient jetzt zur 
Berechnung von P,, das durch die Idealisierung der Bewegung von m, zur Reaktionskraft 
geworden ist. 

Das zu behandelnde System ist einfach, hat aber merkwiirdige, durchaus nicht trivial er- 
scheinende Eigenschaften, deren Herausarbeitung zum Zweck einer physikalischen Anschauung 
der Vorgange beim oszillierenden Pendel die Aufgabe dieses Abschnitts sein soll. Vom Pendel 
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mit oszillierendem Drehpunkt oder gelegentlich auch kiirzer vom oszillierenden Pendel soll 
gesprochen werden, wenn die zwangslaufige Bewegung in einer periodischen Bewegung von sehr 
hoher Frequenz und kleiner Amplitude besteht, etwa nach dem Gesetz x— acoswt, wobei 
dann w z,B. gegeniitber der Frequenz des entsprechenden Schwerependels eine sehr gro®e und a 
(der halbe Hub von A) gegeniiber | eine sehr kleine Zahl jst. 


b) Es sei gestattet, mit einer reizvollen und naheliegenden Abschweifung zu beginnen: mit 
75-0 | Sind Po==0 (4) 

werde die ,,kraftefreie‘‘ Bewegung des Systems kurz untersucht, die bei passenden Anfangs- 
bedingungen eintritt, wenn z.B. die Zeichenebene von Abb. 1 senkrecht zum Vektor der Schwer- 


kraft liegt. Die kinetische Energie und der Gesamtimpuls in der x-Richtung sind jetzt offenbar 
konstant; mit einfacher Rechnung erhalt man 


es C, ne : ekg 70 
Peace) Em. ) 
Die Integration liefert 
Cy (t— t) = E(g, FY, (6) 


wo E das elliptische Integral zweiter Gattung (nach Legendre) zum Modul k? bedeutet, fiir welches 
dieser Fall eines der einfachsten physikalischen Beispiele ist. 

Man findet leicht, und es entspricht wohl der etwas geschulten Anschauung, da8 ¢ im allge- 
meinen nicht Null sein kann; m, erfahrt also von m her Reaktionskrafte in Richtung der Pendel- 
stange; x kann daher nicht konstant bleiben. Dagegen hat die Geschwindigkeit des jeweiligen 
Schwerpunktes von m, und m eine konstante x-Komponente x*, und es wird x= x*-+ x, 
wobei Xx, periodisch verlauft. Die weiteren Einzelheiten sind in bekannter Weise zu berechnen. 

Hervorzuheben ist noch, daf die DrehimpulsgréBbe mi? des Pendelkérpers um den Dreh- 
punkt A trotz P,—0 nicht konstant bleibt. Es existiert kein einfacher Drallsatz ohne Zusatz- 
glied fiir die Koordinate ~, da der Bezugspunkt A mit nicht konstanter Geschwindigkeit, die 
nicht standig der Schwerpunkts-Geschwindigkeit parallel ist, fortschreitet. Aber naturgemaf ist 
der allgemeine Momentensatz fiir A als Bezugspunkt giiltig; aus iam und dem Schwerpunktssatz 
kénnte man die Bewegungsgleichungen ebenfalls gewinnen (und auch aus den Gesetzen fiir die 
Relativbewegung). 

Die erwahnte Besonderheit der Koordinate pm auBbert sich im Ausdruck (1) fiir die kinetische 
Energie dadurch, daB sie hierin explicite durch den Einheitsvektor e, vorkommt, der mit der 
festen x-Achse den Winkel y-- 2/2 einschlieBt. 


c) Von jetzt an sei x als bekannte Funktion von t gegeben und A durch passende 
Einrichtungen zwangslaufig in entsprechender Weise bewegt. 

Eingepragte Krafte fiir die Koordinate ym mégen zunachst nicht vorhanden sein. Damit 
werden die Bewegungen des .,kraftefreien’* Pendelkérpers untersucht, etwa fiir den Fall, daB 
die Bewegungsebene horizontal liegt wie in Abschnitt b. Mit P,—0 wird jetzt 


p= —ésing (: = 7} (7) 


mit bekanntem E(t) die maBgebende Gleichung. 

Sie hat eine anschauliche Bedeutung. Die Pendelstange ist eine Gelenkstange, kann daher 
auf m nur eine durch A gehende Zentralkraft von der Richtung -+-e, iitbertragen. Bei Abwesenheit 
sonstiger Krafte, wie sie hier vorausgesetzt ist, liegt also die absolute Beschleunigung dv/dt 
von m standig in der Richtung der Pendelstange. Nun ist 


LO eet le 
dt ; 


hiervon liegt 21 bereits in der Stangenrichtung; die zu dieser senkrechte Komponente von x 
muf daher entgegengesetzt gleich pl sein, wie (7) es ausspricht. Wenn es daher in einer Literatur- 
stelle hei®t, daB die Bewegung von m in die Stangenrichtung fallt, so ist dies nicht richtig; die 
Geschwindigkeit von m wird im allgemeinen auch Komponenten senkrecht zur Stange haben, 
und es wird sich zeigen, daf diese fiir die bekannten besonderen ,,Effekte’* an diesem System 


maBgebende Bedeutung haben. 
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Es werden nun einige Fille behandelt, die zur Veranschaulichung dieser Erscheinungen bei- 
tragen; es hat sich gezeigt, daB auch bei diesem System die Betrachtung von StoBvorgangen 
instruktiv ist. 

1. Zur Zeit t=t, sei das Systemin Ruhe, etwa mit p=q. In dem darautffolgenden 
Zeitintervall At werde A durch geeignete Krafte auf die Geschwindigkeit 4x gebracht. In 
bekannter Weise werde nun At— 0 angenommen, die Krafte auf A aber so gesteigert, daB Ax 
den in Aussicht genommenen Wert erhalt. Bei geniigend kleinem At wird dann die gleichzeitige 
Lageninderung von m, und auch diejenige von m beliebig klein ausfallen. Mit anderen Worten: 
fiir einen Sto®vorgang dieser Art kann (7) bei konstant bleibendem @ iiber t wihrend At— 0 
integriert werden und liefert 


Ap= — Aé sin (8) 


fiir einen StoBvorgang. 
Diese GesetzmaBigkeit besteht auch, wenn die GeschwindigkeitsgréBen vor dem StoB von 


Null verschieden sind. 

Der Drehkérper erhalt also durch den Sto8 (Hammerschlag od. dgl.) auf den Kérper mj, eine 
plotzliche Anderung der Winkelgeschwindigkeit, obwohl die verbindende Stange nur Krafte in 
der radialen Richtung iibertragen kann. Besonders paradox erscheint dies bei y= 2/2; man 
stellt leicht fest, daB hierbei die in der Stange auftretende Kraft be- 
schrankt bleibt, als Zugkraft vom Betrage m(Aq)*l, entsprechend 
der sofort auftretenden Centripetalbeschleunigung. 

Wird nach dem StoB x konstant ge- 

/ x(t) halten, so bleibt auch @ konstant; der 

Pendelkérper dreht sich dann gleich- 
formig um A. 

Man kann sich in den Vorgang etwas 


— 


fi 
x Golde, wen x 
iy 


ee eee a mehr ,,hineinfiihlen‘‘, wenn man die Ab- 

7 solutbewegung betrachtet. Bei gy = 7/2 

wv | und vollkommener Ruhe vor dem StoB 

ae bewegt sich m auf einer einfachen 
Abb. 2. Eine kriftefreie Abb. 3. StoBbewegung des Zykloide > in Abb. 2 ist dies skizziert ; 
Bewegung. Drehpunktes. drei in gleichen Intervallen aufeinander- 


folgende Lagen 00, 11, 22 der Stange 
sind gezeichnet, die vom freien Endpunkt beschriebene Zykloide angedeutet. 

Man beachte, daB bei AE>0 stets Ap <0 fiir 0<gy<a. Bei yy >a2/2 wird die StoBkraft 
in der Stange eine Druckkraft; dies ist ebenfalls wesentlich. Bei Ersatz der Stange durch einen 
Faden ware die Giiltigkeit von (8) einzuschranken durch die Zusatzbedingung: Reaktionskraft >0. 

Die einfache Tatsache, daB die rechte Seite von (8) den Faktor sing hat, erweist sich spater 
2 ame wichtig; |4@| ist um so gréBer, je gréBer |sing| ist, entsprechend fiir kleinere 

sin g|. 

Nach dieser Vorbereitung wird nun ein Fall behandelt, der bereits die wesentlichen Ziige des 
oszillierenden Pendels aufweist. 


2, A werde nach dem aus Abb.3 ersichtlichen Gesetz hin- und herbewegt: || = a5 == compe, 
aber zu den Zeiten tj, t,,... im Abstand At/2 werde durch einen passenden StoBimpuls die 
Geschwindigkeitsrichtung jeweils bei |x|—a umgekehrt, so daB hier jeweils ein 
Geschwindigkeitssprung vom Betrage 2 vp erfolgt; — 2 v9 bei t,, +2, bei t,. 

Fir t<t, sei p=0, und 0< 9 <a/2. 


Dann ist nach (8) unmittelbar nach dem ersten Sto® auf m, 

“ele Zia A 

P=M1= 7 sing fiir ee Gi Ae (9) 
Mit dieser konstanten Winkelgeschwindigkeit dreht sich der Pendelkérper wahrend des Zeit- 


intervalles At/2; der zweite Sto® trifft ihn auf dem Winkel et . At Date 
(da dt =+2v fir t—t,) P= Por V1 5 Daher is 


er: hy . Zoe oo A 
P2= Pi — 2 Ff sin Y= Gy — = oe (7 +P 2 fiir ly <t<ty. (10) 


Offenbar ist p>; die Anderung von @ beim zweiten StoB fiir tt, ist daher dem Betrage 
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nach groBer als diejenige beim ersten Sto fiir i=t,. Diese war positiv, die zweite ist negativ; 
durch die zwei Sté8e hat also der DrehkérpereinenZuwachs an Winkelgeschwindig- 
keit in negativer Richtung erhalten, namlich (da py=—0 fiir t=0 vorausgesetzt war) 
nach Einsetzen von (9) in (10) 


Ap = Gy = —**Jein (+44) — sin Po . (11) 


Nun mégen ‘die StéBe auf m,in so kurzem Zeitabstand erfolgen, da® der Winkel gp, At/2 sehr 
klein ist; werden dann GréBen zweiter und héherer Ordnung gegeniiber 1 gestrichen, so wird 
nach Kinsetzen von @, aus (9) 
. 2 2 

Ap = — 2 Fy sin py cos py At= —"F sin 2 gy At. (12) 
Die zwangslaufige Bewegung von m, werde nun periodisch fortgesetzt; At ist dann die Zeit- 
periode. Es hat jetzt einen leichtverstandlichen Sinn, (12) durch At zu dividieren und das 
Ergebnis folgendermaBen auszusprechen: 


Der Drehkérper erfahrt wahrend einer Periode die durchschnittliche Winkel- 

beschleunigung 
: , : 
at = — 2% sin gy cos yy = — "8 sin 2 Q. (13) 

Die hiermit angedeutete Vorstellungsweise wird um so besser zutreffen, je gréBer die Frequenz 
des Antriebes von A ist, je mehr StéBe also ausgeiibt werden, ohne daB sich die Lage gy im 
Ganzen wesentlich dndert. Der maB- 
gebende ,,Effekt** der Anordnung be- 
steht darin, daB der Drehkérper durch 
zahlreiche StéBe ein hieraus sich 
summierendes Drehmoment erfahrt, 
das seine mafgebende Mittellinie auf 
dem kiirzesten Wege parallel zur x-Rich- 
tung zu bringen sucht und dem Sinus 
des doppelten Winkels proportional 
ist, ganz dhnlich wie in der kinetischen 
Gastheorie der Gasdruck durch Sum- 
mation der sehr haufigen StoBimpulse 
der Teilchen errechnet und erklart wird. 
Diese Formulierung stimmt offenbar 
auch, wenn in einem anderen Qua- 
dranten liegt. In jedem Falle kann 
man von einer Tendenz zum Par- 
allelismus sprechen. 

In diesem ,,makroskopischen“ Sinn 
hat demnach der kraftefreie Pendel- 
kérper bei oszillatorischem Antrieb des 
Drehpunktes drei Gleichgewichtslagen : 
Po=90, Y= A und Yo=z/2. Es ist 
leicht einzusehen, da davon die ersten 
beiden stabil die dritte labilist. Natur- 
gemaB kann bei dieser summarischen 
Betrachtungsweise nichts dariiber aus- Abb. 4. Schwerpunktsbahn bei StoBerregung. 
gesagt werden, ob die hierbei nicht be- 
achteten Zusatzbewegungen von hoher Frequenz mit t anwachsen, abnehmen oder konstant 
bleiben. 

In Abb.4 ist eine Bewegung der behandelten Art aufgezeichnet. A bewegt sich zwischen A, 
und A, nach dem Gesetz von Abb.3. Die Masse m nimmt mit ~—2/4 nacheinander die Lagen 
0,1, 2,3... 10an; die von A, ausgehenden Lagen der Pendelstange sind ausgezogen, die von A, 
ausgehenden gestrichelt eingetragen. Die anwachsende Winkelgeschwindigkeit nach der Oszil- 
lationsachse hin tritt deutlich hervor. Die Figur ist insofern nicht genau, als die Bahnen 01, 
12, ... tatsichlich zykloidenartig schwach gekriimmt sind, womit sich auch die Lagen der 
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jeweiligen Endpunkte geringfiigig andern. Der Zweck dieser anschaulichen Skizze wird dadurch 
nicht beriihrt. 
Der Anfangswert von g im Punkte 0 ist hierbei mit — “2. sin g derartig gewahlt, dah 


sofort eine durchschnittliche Bewegung nach der Oszillationsachse hin eintritt. Wiirde man 
z.B. Gy 0 wahlen, so witrde mim ersten Teil der Bewegung (mit abnehmender », Winkelgeschwin- 
digkeit“‘) noch in den von den beiden Kreisen um A, und A, gebildeten Zwickel eindringen, aber 
durch die erlauterten Impulswirkungen in dieser Richtung verzégert und zur Umkehr gezwungen 
werden, worauf dann eine Bewegung nach Abb.4 einsetzt. Wird hierbei im ersten Teil gy >2/2 
so tritt keine Umkehr ein, sondern eine nach p= hinstrebende Bewegung. Derartige weitere 
Beispiele sind bei Ersatz der Zykloidenstiicke durch Gerade einfach und schnell zu zeichnen ; 
dies ist um so weniger ungenau, je kleiner a/] angenommen wird, das fiir Abb.4 im Interesse der 
Deutlichkeit nicht gerade sehr klein gewahlt wurde. 


3. In praktisch wichtigen Fallen wird der Antrieb von m, etwa nach dem Gesetz 
x= acoswt (14) 


erfolgen. Die bekannte Theorie liefert hierfiir den Wert 


: 2 2 

a — — + © a sin Po * COS Po = — 45 wt sin 2 y. (15) 
Der Aufbau dieser Beziehung ist genau derjenige von (13). Hierbei ist die Geschwindigkeits- 
anderung wahrend einer Halbperiode Ax —=2aw, im Falle von (13) dagegen 2 v9; danach 
betragt der auf das Quadrat dieser Impulsinderung bezogene Zahlenfaktor in (14) ein Viertel 
von demjenigen in (13). Dies hat seinen Grund darin, da bei sinusartiger Erregung auch » (t) 
einem Sinusgesetz von gleicher Frequenz folgt; die Vorgange spielen sich stetig ab; die mit der 
Anderung des Winkels verbundenen Differenzwirkungen der Impulse verteilen sich gleichmaBig 
und werden dadurch in ihrer Wirkung abgeschwacht. 


4. Bekanntlich hat Herr Meifner die vollstandige Theorie der linearen Gleichung auf der 
Basis durchgefiihrt, da in (14) die Sinuskurve durch passende Parabelbégen ersetzt 
wird. In (8) tritt dann als Zeitfunktion eine jeweils fiir eine Halbperiode konstante Beschleu- 
nigung von wechselndem Vorzeichen auf. Demgegeniiber stellt der Antrieb nach Abb.3 eine 
noch weitergehende Idealisierung dar. In der entsprechenden Hillschen Differentialgleichung ist 
jetzt die periodische Funktion im allgemeinen Null und fiir successive Zeitwerte At/2; At; 3/2 At..: 
unendlich mit wechselndem Vorzeichen und von solcher Machtigkeit, daB eine jedesmal gleiche 
beschrankte Impulsanderung resultiert. 


3, Zur Erweiterung der bisherigen Liésungen, a) In den bisherigen Lésungen hat man sich 
darauf beschrankt, Quasi-Gleichgewichtslagen zu untersuchen, in denen im Durchschnitt z. B. das 
ih Gk , 
= ms Ww? sin Pp COS Po 
die Waage halt. Nach den hier gegebenen Erlauterungen erscheint es naheliegend, eine Erwei- 
terung der Lésungen in folgender Weise zu versuchen. 


Die Gesamtlisung y(t) wird zerlegt, etwa 


VE 8s Heche (16) 
hierin soll nun aber # nicht wie bisher im Mittel als konstant angesehen, sondern dadurch definiert 
werden, daf es eine ,,Grundbewegung“ darstellt, wie sie sich unter dem Einflu8B elngepragter 
Krafte zusammen mit jenem Impulsglied allein ergeben wiirde, wobei dieses ebenfalls als ein- 
gepragte Kraftgréfe behandelt wird. Dabei braucht @ nicht klein zu sein; lediglich # und # 
werden gegeniiber den entsprechenden von den Oszillationen herriihrenden GréBen als geniigend 
klein angenommen. 


Es wird also # definiert durch die Differentialgleichung 


Moment der Schwerkraft der unter 2. veranschaulichten Impulswirkung — - 


- leat ae 
i= — or 7 w sin } cos 0 + k(d) (17) 


unter der Voraussetzung 


d<w, b<a%, (18) 
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worin k(?) elngepragte KraftgréBen bezeichnet, die im allgemeinen von gy und ? abhangen 
werden. Durch die Einschrankung (18) werden u. a. extreme Anfangsbedingungen, wie starke 
StdBe u.dgl. ausgeschlossen, die z.B. den Pendelkérper in stindige Rotation mit groBer durch- 
schnittlicher Winkelgeschwindigkeit versetzen wiirden, Falle, die in diesem Zusammenhang nicht 
von Interesse sind. 

Uber diese Giundbewegung lagert sich x als Schwingungsbewegung von sehr hoher Frequenz 
und kleiner Amplitude. 

Wie im folgenden kurz dargelegt wird, fiihrt dicser naheliegende Gedanke zum Erfolg in 
dem Sinne, daB die noch verbleibenden Restglieder klein von héherer O1 dnung ausfallen und 
mit 0/@-—>0, a/l—0 stirker nach Null gehen als die ».Lésung“ selbst, worauf auch die bisherige 
Theorie basiert. 

Mit dieser Genauigkeit ergibt sich y als eine Schwebung von der Kreisfrequenz w mit einer 
zu sin ? proportionalen Amplitude. 

b) Als erster Fall werde 

k= 0 (19) 


behandelt; die Oszillationsbewegung von A sei dargestellt durch 
X= acoswt; 


mit «= a/l entsteht dann aus (7) die zu lésende Differentialgleichung 


P= a w* cos wtsin Gp. (20) 
Nun sei 
p=O+y. 7 <1; (21) 
und # definiert durch 
— atu? sin } cos } (22) 
mit der Einschrankung 
3b a 
a ee 23 
A Ga ae (23) 


Gleichung (22) ist allgemein mit elliptischen Funktionen und fiir geniigend kleine 3, die immer 
noch groB gegen y sein kénnen, durch sin und cos ohne weiteres zu lisen, und es ist klar, dab 
durch die Einschrankung (23) die praktisch vor allem wichtigen Falle nicht ausgeschaltet werden; 
dies 1a8t sich auch von Fall zu Fall priifen. 
Mit ¥ <1 wird mit der itblichen Genauigkeit 
sin y = sin 0-+ x cos @. (24) 
Fir y werde nun angesetzt 


4=— asin Pcoswt+y, (25) 


wobei durch das zunachst unbekannte y die verbleibenden Restglieder usw. erfaBt werden. 


Bei der Bildung von ¥ aus (25) erscheinen Glieder mit 8 und @; diese sind aber nach (23) zu 
streichen [Gleichung (24) beriicksichtigt auch nur die erste Ordnung], und es bleibt 


%¥ = aw sin Icoswt+y. 


— 


. 1 ; ns 
Setzt man alles in (20) ein und beriicksichtigt cos? m t= —> =e os cos 2wt, so verbleibt fiir y 
a 


die Differentialgleichung 


ip — 0 02 cos H cos wt-y = — 5 aw sin Pcos P cos 2 wt. (26) 


Dies ware eine gewohnliche inhomogene Mathieusche Gleichung, wenn # konstant ware. Das 
ist nicht der Fall. Entwickelt man aber # von einem gegebenen t aus nach Potenzen von t, was 


zweifellos erlaubt ist, so erscheinen wieder Glieder mit @ und @, die nach (23) im Rahmen der 
gesamten Genauigkeit zu streichen sind, so daB (26) in diesem Rahmen als Mathieusche Gleichung 


zu diskutieren ist. 
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Fir die itbliche partikulare Lésung ist die rechte Seite von (26) maBgebend. Diese ist, da 
auch « als Kleine Zahl zu gelten hat, von einer Ordnung kleiner als diejenige von (20). Der 
entsprechende Anteil von y ist daher um eine Ordnung kleiner als y nach (25) und daher (ab- 
gesehen von Resonanzfallen od.dgl.) im Rahmen der allgemeinen Genauigkeit zu streichen. 

Es verbleibt der homogene Anteil von y, maBgebend fiir die Stabilitat ,,im kleinen‘*. Diese 
wird zum SchluB in einem gemeinsamen Abschnitt behandelt; unter den iiblichen Verhaltnissen 
ist die homogene Lésung von (26) stabil. 

Hiernach setzt sich die kraftefreie Bewegung des Pendelkérpers im Rahmen der iiblichen 
Genauigkeit zusammen aus einer langsamen Schwingung um die im ,,Grofen“ stabilen Gleich- 
gewichtslagen 0— 0 und )=z, deren Frequenz durch (22) und bei gréBerem ? durch die Anfangs- 
bedingungen gegeben ist, und aus einer Schwebung nach (25) mit der Frequenz w und der 
Amplitude « sin #. 

Fiir kleine Werte von ? mit sin? ~ # und cos ? ~ | 1aBt sich alles leicht ausrechnen; 
man erhalt dann z. B. 


ios se aaael 
ee fiir sin Hy ~ Jo. (27) 
y= — 0%, sin xwt cos Wt | 


Das Verhaltnis der Grund- zur Schwebungsfrequenz ist dann \z a. In Abb.5 und 6 ist dies fiir 


VX 


Abb. 5. Hauptbewegung und Schwebung Abb. 6. Schwerpunktsbahn bei kraftefreiem Pendel. 
bei kraftefreiem Pendel. 


a = 1/10 skizziert. Abb.5 zeigt @ und y als Funktion von t; in Abb.6 ist die Bahnkurve von m 
angedeutet, wobei der Kleinheit von %) entsprechend die schwach gekriimmten Kreise um A 
als gerade Linien angenommen wurden. Die Gleichung der Kurve ist durch die angegebene 
Lésung ohne weiteres in t als Parameter gegeben, samtliche Einzelheiten sind leicht zu disku- 
tieren. (In der Skizze Abb.6 ist zur Vereinfachung 10 y2 ~ 15 benutzt worden.) 

Hiernach kann ein kraftefreier Drehkérper, etwa das Pendel 1 in Abb.7, 
als Anzeigeinstrument fiir die Richtung der Oszillation des Dreh- 
punktes benutzt werden. Ein gleiches Pendel 2, zu dem ersten unter 180° an- 
geordnet, tut dasselbe und hat bei gleichen Anfangsbedingungen dieselbe Grund- 
schwingung um die Oszillationsrichtung. Man kann aber die beiden Pendel 
nicht zu einem starren Kérper verbinden und dem Instrument dadurch etwa die 
doppelte Richtkraft und Unabhangigkeit z.B. von der Schwere bei anderer Lage 
der Bewegungsebene geben. Denn durch eine starre Verbindung entsteht einfach 
ein ausbalanzierbar starrer Kérper, der sich um eine Schwerpunktsachse drehen 
kann, und die Bewegungsgleichung in y wird unabhangig von derjenigen in x. 
Die beiden Pendel hatten jetzt stets dasselbe y, @ usw., und die Richtung von 
x ist standig fiir das eine giinstig, fiir das andere ungiinstig, und umgekehrt im 
Sinne der Tendenz zum Parallelismus; die Wirkungen auf das eine neutralisieren die gleich- 
zeitigen auf das andere. 


Abb. 7. Eine 
Doppelanordnung. 


c) Nach Abb.8 werde nun die Schwerkraft bei horizontaler Oszillationsrichtung 
beriicksichtigt. 


Mit w= g/l, wobei wy die Kreisfrequenz des entsprechenden Schwerependels bei kleinen 
Schwingungen um die Vertikale bedeutet, wird 


k= a} cos y, 
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und die zu lésende Differentialgleichung lautet 
PY = aw* cos wt sin p + Ww? cos —/. (28) 
Wieder sei 
g=0-+y. sing=sin d+ x cost; 
aber & wird jetzt definiert durch 
a Faw? sin } cos U-+ we cos U. (29) 


Das hat nur einen Sinn, wenn ws von gleicher Ordnung wie a?w?, also gegeniiber w? von zweiter 
Ordnung klein ist; dies wird vorausgesetzt. 
Fir x ergibt sich die Differentialgleichung 


c— = ooo? sind cos? + a w* sin? coswt+ aw cos? coswt.y— wasind «x, (30) 

die wieder durch den Ansatz 
Y=—asinicoswt+yp, Yra-w'sind coswt-+ wp (31) 
gelést wird, wobei in ¥ die Glieder mit # und @ korrekterweise zu streichen waren. Fiir p verbleibt 
y+ (aj sin Bd — aw? cos # coswt) p= w2 sin? 9 cosmt — E a?w* sin? cos2wt. (32) 


Die rechte Seite hiervon ist wieder um eine Ordnung héher als die 
entsprechenden Glieder in der Gleichung (30) fiir y, die partikulare 
Lésung von (32) ist also zu streichen. 

Wie im vorhergehenden Abschnitt gezeigt wurde, ist # in der 
verbleibenden homogenen Gleichung als ,,sehr langsam verander- 
licher Parameter zu betrachten, die Gleichung daher im Rahmen 
der gesamten Genauigkeit als Mathieusche zu diskutieren, in welcher 
jetzt auch ein konstanter Faktor von wy auftritt. Die hiermit ver- 


RY) 


bundene Untersuchung der Stabilitaét im kleinen erfolgt zusammen- Abb. 8. Schwerewirkung, 


Oszillation horizontal. 
fassend zum SchluB. 


Fir @ ergibt (29) wieder die Méglichkeit von Gleichgewichtszustanden, etwa fir ? = 0); 
namlich fiir cos J) = 0. Jj) = +27/2; (also fiir die beiden vertikalen Lagen) und fiir 


; 2m? 

S10 y= rice (33) 
Die rechte Seite ist stets positiv; damit aber sin %, reell ausfallt, mul offenbar 

1 ‘ 

= oF w* > wf (34) 


“ 


sein. Wenn dies der Fall ist, dann gibt (33) die bekannten Gleichgewichtslagen, namlich bei 
5 oto? > we zwei zur Vertikalen symmetrisch liegende, und fir = a? w? = w? die vertikale 
Lage } = + 27/2. 

Durch eine einfache Betrachtung der bei Abweichungen aus diesen Gleichgewichtslagen 
auftretenden Anderungen von i a2w* sin Pcos @ und wi cos @ stellt man elementar fest, dab 


im ,,GroBen“ 


stabil | 1 
die Lage J) = + x/2 indifferent ; ist fiir > a? jw? = we, 
|labil (35) 
die Lage %) = — 2/2 stets labil, 
; 2 
die Lagen sin J) = — <1 stabil 


sind. 

Bei geniigend groBem w? wird # nach (33) beliebig klein, dann ist also, wie bekannt, die im 
GroBen stabile Pendellage mit beliebiger Annaherung horizontal. 

Der endgiiltige Vergleich mit den bisher gewonnenen Ergebnissen wird bei Diskussion der 


homogenen Gleichung durchgefiihrt. 
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d) Bei vertikaler Oszillationsrichtung im Schwerefeld (Abb. 9) tritt 
~ = 4? cos wt sin p — We sin P (36) 
an die Stelle von (28); # soll jetzt die Gleichung 
j= — - a2w? sin 3 cos # — w? sin (37) 
erfiillen; fiir y ergibt sich an Stelle von (30) 
Ve Fate sin 9 cos 0+ a w* sin 3 cos wt — w2 cos Oy + aw? cos P cos wt: x, (38) 
und mit dem gleichen Ansatz (31) fiir y verbleibt fiir y 
y+ (w2 cos 0 — aw* cos & cos wt) y = 
= aw? sin ? cos Pcos wt — 2 a?w* sin ? cos P cos 2wt, (39) 
wovon wieder die partikulare Lisung zu streichen ist. 
Die Gleichgewichtslagen fiir # sind jetzt offenbar 
2w2 
Uo 0) SU) ol Und ey abccos ia al ; (40) 


a2 m2}? 


bei = a? uw)? > 2 w? gibt es jetzt also auBer den beiden vertikalen Lagen noch 


je zwei symmetrische in der oberen Halbebene, in Abb.9 strichpunktiert 
angedeutet. Kine elementare Betrachtung gibt wieder Auskunft iiber die 
Abb. 9. Schwerewirkung, ee ; we 
Oszillation vertikal. Stabilitat ,,im GroBen‘: 


Gleichgewichtslage %)— 0 stabil fiir alle w’, 


stabil seme drat Sate 
Gleichgewichtslage 3, — z | indifferent fiir 9 ee S Wo> (41) 
| labil 
Gleichgewichtslage 2/2< &j< 7 labil. 
y e) Die noch zu behandelnden homogenen Glei- 
chungen fiir y werden mit wt—2z auf die haufig 

y 20 benutzte Normalform 
Uy dy 
Ys qq + (A+ keos2 t) y= 0 (42) 
Y gebracht. Das hier benutzte k entspricht 2h? bei Strutt}, 
Yy y Klotter und Kotowski? benutzen A, y, x anstatt A, k, 27 


in (42). Die hier in Frage kommenden Werte von k und 
A sind klein; fiir die Beurteilung der Stabilitat kommt 
es in einer Darstellung k(/) auf die Umgebung des Punk- 
tes A= 0 an. In Abb.10 sind die hier interessierenden 
Gebiete dargestellt; in dem schraffierten besteht Stabili- 
tat, auf den Grenzkurven und in dem freigelassenen 
Gebiet Labilitat, soweit es auf die Lésung der homogenen 
Gleichung ankommt. 

Da es sich hier um angendherte Lésungen handelt, 
bei denen nur von erster Ordnung kleine Glieder beriick- 


G2 


Abb. 10. Gebiete der Stabilitat und der Labilitat 


»im kleinen‘, sichtigt werden, geniigen auch fiir die Grenzkurven die ein- 
fachsten Abschatzungen, die in k bis zu k? gehen und in 

A die jeweils kleinsten Glieder beibehalten, namlich m4 
|k| = 1,80(1—0,91)A und ~ {k| = y8(—A) , (43) 


wovon die erste Gleichung die obere Grenzkurve dieses Bereiches fiir / = 0 darstellt, die zweite 
die untere fiir A< 0. Auf der oberen ist k nicht mehr klein; trotzdem ist die Ubereinstimmung 


és M. J. O. Strutt, Lamésche, Mathieusche und verwandte Funktionen in Physik und Technik. Berlin 
1932. 


* K. Klotter u. G. Kotowski, Z. angew. Math, Mech. 19 (1939) S, 289. 
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des angegebenen Wertes z.B. mit der bekannten sog. Struttschen Karte bereits vollauf geniigend ; 
die Beriicksichtigung des Gliedes mit k* liefert Korrekturen, die hier nicht von Bedeutung sind. 

Hieraus ergibt sich die zu dieser Mitteilung gehérende ,,Tabelle fiir die Stabilitat nach der 
homogenen Gleichung“*. Das Hauptergebnis ist: Bei geniigend groBem w und passendem 
(weder zu groBem noch zu kleinen ) a ist in allen Fallen Stabilitat der Bewe- 
gung nach der homogenen Gleichung zu erzielen. 

Diese Bedingungen sind nun zu vergleichen mit denjenigen, die sich bei Existenz stabiler 
Gleichgewichtslagen 9 = 0, ,,im Gro®en“ fiir diese ergeben haben: 

Im Falle a (kraftefreies Pendel) ist 9,— 0 stabil, ohne besondere Bedingung hierfiir. In diesem 
Fall wird also durch die homogene Gleichung laut Tabelle 1 « nach oben begrenzt; das ist 
aber unwesentlich, da a< 0,45 stets erfiillt sein muB, wenn die ganze Annaherungsrechnung 
einen Sinn haben soll; sie verlangt a <]. 


Tabelle 1. Fiir die Stabilitét nach der jeweiligen homogenen Gleichung. 


N Abschnitt Nr. der Kennzeichnung 7 Stabilititsbedingung 
T.  |dieser Arbeit Gleichung des Falls [K| (fiir die homogone Gleichung) 
a |») EG ee 0 4 |cos 8| | 0 <alcos 8|< 0,45 
EE ee a 
Schwerkraft 4 aii O 
1 c) (32) 4 0.-Richtung @? 40 |cos B| | 0<a|cos B|<0,45 (1—0,91 A) 
| | horizontal sin ? >0 | 
[Sy sae 
| oe cin b 0) /2|sind O| <0,45 
2 Ce eG) o wo? 4a |cos 8| | = // 2lsin Al<4 leos 4] <0, 
| | sin O< 0 (1 — 0,91 A) 
Schwerkraft | 405 eaatp 
1 d) | (39) | 0.-Richtung | © 4.0 |cos 3| | 0 Sa|cos d|< 0,45 (1—0,91 A) 
vertikal | cos0>0 
4 wo? Wo RT a 
ae 008 0 Wee ri. 2|sin B\<a | cos #|< 0,45 
cos 0< 0 (1 — 0,91 A) 


Im Falle b 1 (Schwerkraft, 0.-Richtung horizontal, untere Halbebene) gibt es nach (33) 


sin Hy 
2 


und (35) stabile Lagen, wenn —?— a|/ ist; diese Bedingung kommt also gegebenenfalls 


zu derjenigen fiir « nach der Tabelle hinzu. Hier bildet #—2/2 eine Ausnahme; diese Lage ist 
nur stabil, wenn 


ee gf 

@ > o | Ded 
Im Falle b 2 (wie vorstehend, aber obere Halbebene) gibt es keine stabilen Lagen fiir o. 
Im Falle c 1 (Schwerkraft, 0.-Richtung vertikal, untere Halbebene) ist %)—0 stabil im 


GroBen, ohne Bedingung hierfiir; Einschrinkung von a lediglich durch die Bedingung der 
Tabelle. 


; o iL 
Im Falle c 2 (wie vorstehend, aber obere Halbebene)ist nach (41) 09a stabil, wenn rf < al/> ; 


dies ist in diesem Falle die gleiche Bedingung wie diejenige nach der Tabelle fiir die Zusatz- 
aes folgt: Der stabile Bereich fiir die Lésungen der homogenen Gleichungen 
enthalt alle stabilen Gleichgewichtslagen (im » GroBen") und auSerdem weit- 
gehende Méglichkeiten fir stabile ,,Grundbeweguugen”, wie sie jeweils Corals }(t) 
mit #<w und §<w? dargestellt werden und iiber welche sich die durch 
—=—asin ?coswt dargestellten Schwebungen lagern. oat 
f) Eine Verfeinerung der Theorie ist z.B. in fofgender Weise miglich : O(t), wofitr hier 
nur die periodischen Bewegungen beriicksichtigt werden sollen, kann immer als Fouriersche 
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Reihe angesetzt werden; in vielen Fallen werden wenige Glieder hiervon vollkommen ausreichen. 
Hiermit und mit dem Ansatz fiir y nach (25) wird dann in die Differentialgleichung fiir y hinein- 
gegangen, in welcher nun aber alles beibehalten werden mu, was dem gewiinschten Genauig- 
keitsgrad entspricht. Es hat also wenig Sinn, Glieder hoherer Ordnung z. B. bei einer partikularen 
Lésung fiir y zu diskutieren, nachdem man sie durch sin y= sin &-+ y cos @ an anderer Stelle 
des Gedankenganges gestrichen hat. Bei konsequentem Vorgehen erhalt man so die Méglichkeit, 
auch z.B. in ¥ das bisher konsequenterweise gestrichene # und # sinngemaB zu beriicksichtigen. 
Zum SchluB entsteht so eine Differentialgleichung fiir y, die wie bisher linear ist, und deren 
Faktor von y ebenfalls periodische Funktionen von t enthalt, mit dem wesentlichen Unterschied, 
daf hierin jetzt zwei Perioden auftreten, namlich aufer 2 2/@ noch diejenige von 0. Sehr haufig 
werden diese nicht in einem rationalen Verhaltnis stehen; man hat es dann weder mit einer 
Mathieuschen noch mit einer allgemeinen Hillschen Differentialgleichung zu tun. Fiir die An- 
wendung wird man aber haufig geniigend genau und sicher vorgehen, wenn man in solchen Fallen 
ein benachbartes rationales Periodenverhdltnis fiir die rechnerische Abschatzung zugrunde legt 
(also z. B. 14 oder 15 anstatt 10 72), womit dann wieder eine Grundperiode erzielt wird und die 
bekannten SchluBweisen fiir die Hillsche Gleichung zur Verfiigung stehen. In dieser und anderer 
Richtung gibt Hamels! Arbeit zahlreiche Ausfiihrungen und Anregungen nebst der notwendigen 
Kritik. 

Der Verfasser hat sich iiberzeugt, daB dieses Verfahren zum Ziele fiihrt; zur Veréffentlichung 
sind die naturgemafi etwas umstandlichen Rechnungen im Ganzen wenig geeignet. 

g) Dagegen erscheint es lohnend, auf eine hierbei auftretende Teilaufgabe in kurzen Ziigen 
einzugehen. Die zu lésende Differentialgleichung sei 


2 
= (A + ksin - cos 2 1] c= 0 (m = ganze Zahl). (44) 


Gegeniiber der gewéhnlichen Mathieuschen Gleichung 
d*x 
qa t+ (A+ koos2t) x= 0 (45) 


ist also dic Amplitude des periodischen Faktors mit cos 27 nicht konstant, sondern auch noch 
periodisch veranderlich, und zwar mit einer anderen Periode. Diese Aufgabe tritt auf, wenn man 
entsprechend den Andeutungen des vorigen Abschnittes auch die von héherer Ordnung kleinen 
Glieder beriicksichtigen will, und zwar ist dann das in (44) allein herausgegriffene Glied eines 
unter mehreren anderen, denen gegeniiber k sehr klein ist. Entsprechend der Aufgabenstellung 
dieser Mitteilung wird m eine verhaltnismafig groBe Zahl sein, doch laBt sich (44) auch fiir 
kleinere m ohne weiteres lésen, wenn es nur als ganze Zah] angenommen wird. Dann ist namlich 
der Faktor von x wieder periodisch (mit der Grundperiode mz); (44) ist eine Hillsche Differential- 
gleichung, und die bekannten Methoden zur Lésung und zur Beurteilung der Stabilitat stehen 
zur Verfiigung. Das Ergebnis einer im wesentlichen in bekannter Weise durchgefiithrten Rech- 
nung fiix die Grenzkurve bei A <0, die hier hauptsachlich interessiert (entsprechend der zweiten 
Gleichung (43), ist folgendes. 


Die Grenzkurve fiir 2 <0 hat die Gleichung 
|k| = Ye(—A); | 


@== 32 fir m— 1. 


hierin ist 
6 G8 n (46) 
=a Gey Hae FT SSS AN | 


c—>16 fiir m—> o. 


Die folgende Tabelle 2 enthalt einige Werte von cin Abhangigkeit von m. Fiir m> 3 wird 
also das fiir die Stabilitat erforderliche k vergréBert, wenn diese AmplitudengréBe nicht 
konstant ist, wie in (45), sondern entsprechend (44) periodisch um Null verdnderlich. Dies ist 
plausibel, da nach (44) im Durchschnitt ein kleinerer Amplitudenwert fiir cos 2 1 wirksam ist. 
Fir die kleineren Werte von m verliert diese summarische Anschauung ihren Sinn; die Gesamt- 


1 G. Hamel, Math. Ann. (1912). 
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bewegung laft sich dann nicht mehr als eine Modifikation von (45) mit langsam veranderlichem k 
auffassen. 


Tabelle 2. 


c = 8 fiir Gleichung (45) 


Der Wert fiir m=1 fiigt sich nicht in diese GesetzmaBigkeit. Das zeigt sich auch darin, dab 
: ; : 2 2 ee Cee 
bei m= 2 die gleichen Werte fiir c erhalten werden, wenn man in (44) cos = anstatt sin — ein- 


setzt, wahrend hierbei der Faktor von x mit m—1den Wert A = k+ = k cos — erhalt, also 
etwas wesentlich Verschiedenes. 


Ks lage nahe, in (44) cos a anstatt sin ee zu schreiben und dann mit m—> co den Ubergang 


zu konstanter Amplitude von cos 27 zu versuchen. Dies fiihrt aber nicht zum Ziel; denn zur 
Aufstellung der Stabilitatsbedingung, die zu Gleichung (46) fiihrt, wird jedesmal die gesamte 
Periode mz herangezogen, die mit m-—> oo ebenfalls nach Unendlich geht. 


(Eingegangen am 10. Mai 1951.) 
Anschrift des Verfassers: Professor W. Kucharski, Berlin-Charlottenburg, Hardenbergstrafe 34. 
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Bemerkung 
zu meinem Aufsatz in Band XVIII, S. 363 des Ingenieur-Archivs: 


»Die Biegeschwingungen eines Stabes unter pulsierender Achsialkraft 
bei beliebigen Randbedingungen“ 


Von K. Klotter. 


Von den Herren W. Kucharski und I. Szab6 in Berlin erhielt ich eine ausfithrliche, von Herrn 
S. Woinowsky-Krieger in Quebec eine kurze Zuschrift mit Einwanden gegen die Behandlung 
der Stabilitatsfrage in dem genannten Aufsatz. Ich habe mich von der Richtigkeit der Einwande 
itberzeugt und muB feststellen, daB die Schliisse, die ich hinsichtlich der Stabilitatsbedingungen 


des Stabes mit nicht gelenkig gelagerten Enden zog, unhaltbar sind. 


(Eingegangen am 23, Juli 1951.) 
Anschrift des Verfassers: Professor Dr.-Ing. K. Klotter, Stanford (Calif.), 621 Alvarado Row. 


Berichtigung 
zu meinem Aufsatz in Band XVIII, $.344 des Ingenieur-Archivs: 


»Systematische Sechskomponentenmessungen an Pfeilfliigeln“, I. Mitteilung. 
Von W, Jacobs, 


In der genannten Arbeit sind irrtiimlicherweise in Tabelle2 (Seite 353) die Werte yy nicht auf 
b/2, wie angegeben, sondern auf 6’/2 bezogen worden, auBer fiir Absatz 3 (Rechnung mit kon- 
stanter Cy s,,-Verteilung) und Absatz 4 (Rechnung mit Mittelwert von konstanter ¢,5,-Ver- 
teilung und elliptischer Auftriebsverteilung). Die angegebenen Zahlenwerte fiir yy/(b/2) sind 
also mit b’/b (siehe Tabelle auf Seite 347) zu multiplizieren, um die richtigen Werte zu erhalten. 
AuBerdem sind fiir Absatz 4 (Rechnung mit Mittelwert von konstanter Car Verteilung und 
elliptischer Auftriebsverteilung) folgende neuen Werte fiir yy/(b/2) und xjy/l,, einzusetzen: 


(Eingegangen am 14, November 1951.) 


Anschrift des Verfassers: Dr. habil. W. Jacobs, Stockholm (Schweden), Flygtekniska 
forséksanstalten, Ulvsunda. 
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Die Edelstihle. Von Professor Dr.-Ing. Franz Rapatz, Stahlwerk Gebr. Buhler & Co. - 


A.-G. Kapfenberg (Steiermark).. Vierte, verbesserte und erweiterte Auflage. Unter 
Mitwirkung von Dr.-Ing. Helmut Krainer und- Dipl.-Ing. Josef Frehser, Stahlwerk Gebr. 
Bohler & Co. A.-G., Kapfenberg (Steiermark). Mit 338 Abbildungen und 121 Zahlen- 
tafeln. V, 730 Seiten. 1951.  ~ 5 f Ganzleinen DM 49.50 


Aus den Besprechungen: 

In der Einleitung prizisiert der bekannte und erfahrene Stahlfachmann den Begriff ,,Edelstahl .... 

In ansprechender Weise unterrichtet der folgende Abschnitt tiber die Eigenschaften, Behandlungsvorschriften, Hig- 
nung und Verwendung der Stihle nach Legierungselementen geordnet. Konkrete Angaben wie Vorschriften dienen 
nicht nur zum Kennenlernen eines Stahles, sondern eignen sich ebenso fiir die praktische Anwendung, sei es bei der 
Konstruktionsberechnung, bei der Warmebehandlung oder beim Einkauf oder Priifen entsprechender Werkstoffe. 


Im folgenden, mit ,,Stahle nach Verwendungsgebieten geordnet* betitelten Abschnitt geht der Autor mehr als nur 
auf die allgemeinen Werkstoffeigenschaften, wie Zugfestigkeit, Streckgrenze, Elastizititsmodul etec., ein, indem der 
Kerbwirkung, Dauerfestigkeit, Belastbarkeit, Gestaltfestigkeit und wie die den neueren Erkenntnissen entsprechend 
benannten Beanspruchungen heifen, besondere Aufmerksamkeit géschenkt ist. Der Leser gewinnt so ein zuverlassiges 
Bild.um Méglichkeiten und Grenzen in der Anwendung von warm- und dauerstandfesten, chemisch bestindigen, 
unmagnetisierbaren Stahlen, solchen fiir tiefe Temperaturen, fiir Ventile, Turbinenschaufeln, Zerspanungswerkzeugen, 
den Heizleiterlegierungen und wie die ,,Edelstihle‘ letztlich in der Praxis ihre Einteilung finden. SchlieBich ver- 
dienen die miteingeschlossenen Erzeugnisse nach Sonderyerfahren (StahlguG, gesinterte Stahle, Hartmetalle und 
Hartlegierungen, Hartverchromungen etc.) ebenso gebiihrende Beachtung wie die auf die’ Praxis zugeschnittenen 
Ausfiihrungen tiber bestimmte Gebrauchseigenschaften, wie Zerspanharkeit, VerschleiBwiderstand, Wiarmeausdeh- 
nung, Metallurgie des Schweifens etc. In einem besonderen Abschnitt werden schlieBlich die infolge unrichtiger 
Behandlung entstehenden Stahlfehler aufgedeckt und hierzu fachmannisch Stellung genommen. 


Zusammenfassend ist zu'sagen, da8 es sich hier um ein ausgezeichnetes Werk mit bekannt vorziiglicher Springer- 
Verlag-Betreuung handelt. Insbesondere gilt es auf die zahlreichen Mikrophotos, die graphischen Darstellungen, auf 
den jedem Abschnitt eingeflochtenen einschlagigen Literaturnachweis und das 20 Seiten umfassende Sachverzeichnis 
hinzuweisen. Wer sich je im engeren oder weiteren Sinne mit Edelstahlen zu befassen hat, wird nach einem Versuch 
mit besonderer Vorliebe wieder zu diesem vorziiglichen Buch greifen, ; >» 1 echnische Rundschau‘‘ 


% 


a e 
Hochwertiges GuBeisen, (GrauguB) seine Eigenschaften und die physikalische 


Metallurgie seiner Herstellung. Von Dr.-Ing. habil. Eugen Piwowarsky, ord. Professor der 
Eisenhiittenkunde, Direktor des Instituts fiir allgemeine Metallkunde und das gesamte 
GieBereiwesen an der Rheinisch-Westfalischen Technischen Hochschule Aachen, Zweite, 
verbesserte Auflage. Mit 1063 Abbildungen. XII, 1070 Seiten. 1951. 

' Ganzleinen DM 135.— 


Inhalt at bersicht: I. Einleitung. — II. Die konstitutionellen Grundlagen der Eisen-Kohlenstoff-Legierungen. - 
— III. Ubermolekulare Higenheiten kyhlenstoffhaltiger Lésungen. — IV. Der EinfluB des Siliziums auf die Gleich- 
ge wichts- und Graphitisierungsvorgange. — V. Der Mechanismus der Graphitisierung siliziumhaltigen GuBeisens. — 
Vi. Die strukturelle Beherrschung der metallischen Grundmasse. — VII. Die Primarkristallisation des GuBeisens. — 
VIII. Der Einflu& der standigen Eisenbegleiter. — IX. Der EinfluB der Gase. & X. Technologische Eigenschaften 
flussiger und erstarrender Eisen-Kohlenstoff-Legierungen. — XI. Die mechanischen und elastischen Eigenschaften 
des GuBeisens. — XII. Die physikalischen Eigenschaften des CuSeisens. — XIII. Verhalten des Gufeisens bei hohen 
und tiefen Temperaturen. — XIV. Technologische Eigenschaften des festen GuBeisens,— XV. Die chemischen Higen- 
schaften des Gufeisens. — XVI. Der zusatzliche Oberflachenschutz yon GuSeisen. — XVII. Festigkeitseigenschaften 
yon GrauguS und Temperguf nach Vorkorrosion. — XVIII. Der Hinflu8 thermischer Nachbehandlungen auf me Ge- 
fiigeinderungen und die Eigenschaften perlitischer Grundmassen (Hartungs- und Glithprozesse). — XIX. Legiertes 
GuBeisen. — XX. Die spanlose Verformung des Gufeisens. — XXI. Das SchweiBen von GuSeisen. — XXII. Fehbler- 
quellen bei der Gefiigeuntersuchung yon Gufeisen. — XXII. Das Schmelzen von Gufeisen im Kupolofen, — XXIV- 
Andere Schmelzéfen. — XXV. Einige besonders wichtige Anwendungsgebiete fiir GuBeisen. — XXVI. GieBen und 
SchweiBen. — XXVII. Anhang. — Sachverzeichnis. ; 
Das einschlagige in- und auslindische Schrifttum, dessen Studium dem Fachmann ein abgerundetes Bild vom hye 
der Technik hinsichtlich Herstellung, Eigenschaften uni Hinsatzméglichkeiten des heutigen hochwertigen GuBeisens 
zu vermitteln vermag, ist so tiberaus umfangreich geworden, da nur wenige Fachlewte Zeit und MuBefinden diirften, 
davon in vollem Umfange Gebrauch zu machen. Hier setzt der Zweck des vorliegenden Werkes ein. Hp Role sg Fach- 
mann in geschlossener Buchform einen Uberblick bringen tiber die Metallurgie des GuBeisens, eg s Beziehungen 
zwischen dem Gefiigebau und den Eigenschaften desselben sowie tiber den zweckmafPigen Einsatz dieses hikabilis im 
und seine kinftigen Entwicklungsméglichkeiten. Die Darstellung uud kritische Sabet des bhai ea aiad er- 
folgte in einer Form, ‘die den Leser immer wieder auf offene Probleme hinweist und ihn immer yards za eigenem 
produktiven Denken anregt, um auf diese Weise der weiteren Entwicklung dieses Werkstoffes zu dienen. 


Uber die erste, 1942 erschienene Auflage schrieb die Zeitschrift des Vereins deutscher Ingenieure: 
Es ist ein Verdienst Piwowarskys, da er mit dem vorliegenden Werk eine schr umfassende, auf den letzten Stand 
‘cebrachts kritisch wiirdigende und zur weiteren Mitarbeit anregende Zusammenfassung des gesamten Wissens- 
; - of “eé 
gebietes gibt, das sich bislang in zahllosen Hinzelveréffentlichungen zerstreut fand. 
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